Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,576,080	0	G	A	41.8%	104.7 / 24.6	67	L514F (CTT→TTT)	yddB	putative TonB‑dependent outer membrane receptor
Reads supporting (aligned to +/- strand):  ref base G (33/6);  new base A (0/28);  total (33/34)
Fisher's exact test for biased strand distribution p-value = 2.70e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGAC  >  NC_000913/1575933‑1576228                                                                                                                                                    |                                                                                                                                                     cAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAat                                                                                                                                                   >  1:121247/1‑151 (MQ=255) cAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAat                                                                                                                                                   <  1:77655/151‑1 (MQ=255)          tCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGc                                                                                                                                           <  2:670291/150‑1 (MQ=255)             tACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCAt                                                                                                                                        >  1:971241/1‑150 (MQ=255)                cttttttCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCTTTTTAAAATCCTGATAACGCGTCAGAGTTTTATTCCCTGATACCGATTCCGTCCAGCTATGGCGGATATCACGTAACCCCATATAAAAAATATTCCCGCCAt                                                                                                                                        <  2:967575/144‑1 (MQ=255)                cGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGt                                                                                                                                     >  2:820263/1‑150 (MQ=255)                   ttttCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACg                                                                                                                                 >  1:797787/1‑151 (MQ=255)                      tCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGtt                                                                                                                              >  2:939916/1‑151 (MQ=255)                          ctgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTtata                                                                                                                           >  1:583754/1‑150 (MQ=255)                           tgcAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAAc                                                                                                                         >  1:836673/1‑151 (MQ=255)                              aaccccATTGCGAGTTCATCGTTATAAGGCTTTTTAAAACCTTAATACCGTGTCAGTTTTTTATTCCCTAACCCCAATCCCGCCCACCTATTCCGGATATCATTTAACCCCAAATCAAAAATTTCCCCGCCATAGTAACGGTTATaa                                                                                                                          <  2:422405/144‑1 (MQ=255)                                 cccattcccattTCATCTTAATAAGTCTTTTTAAAATCCTAATACCGCGTCAGATTTTTATTCCCTGATCCCTTTTCCGTCCAGCTATTGCGGATACCCCGTAACCCAATATCAAAAATATTCCCGCCATAGTAACGGTTATAACct                                                                                                                       <  1:431245/147‑1 (MQ=255)                                 cccATTGCCAGTTCTTCTTTATAAGGCGTTTTCAATTCCTGATAACGCGTCAGATTTTTATTCCCTGATCCCGATTCCGCCCAGCTATTGCGGATATCACGAAACCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACctgc                                                                                                                     <  1:434079/149‑1 (MQ=255)                                   cATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGt                                                                                                                  >  1:633729/1‑150 (MQ=255)                                        cctgTTCATCTTTATAAGGCTTTTTAAAATCCTGATAACGCGTCAGATTTTTTTTCCCTGATACCGTTCCCGCCCAGTTTTTGCGAATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCa                                                                                                             <  1:806287/147‑1 (MQ=255)                                           tttcttcTTTAAAGGGCTTTTTAAAATCTTAAAAACGCTTAAGATTTTTTTTCCCTAATCCCATTTCCTCCCACCTTTTGGGGATATCATAAAACCCCAAATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTg                                                                                                          <  2:669105/149‑1 (MQ=255)                                            ttctttgtaataCGGTGTTTTCAACTCCTAATAAGGCTTCAGATTTTTATTACCTGATCCCGATCCCTCCCACCTATTGCGGATATCACGTACTCCCATATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTg                                                                                                          <  2:721332/140‑1 (MQ=255)                                                tttttaaagggCTTTTTAAAACCCTAAAAACGCGTAAGATTTTTTTTCCCTGATACCGATTCCGTCCAGTTTTTGCGGATATCACTTAACCCCATAAAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGtt                                                                                                     <  1:922808/141‑1 (MQ=255)                                                tttttATAGGGTTTTTTAAACTCTTAATACCGCGTAAGTTTTTTATTCCCTGATCCCGATCCCCCCCACTTTTTGGGAATATCTCTTAACCCAATATAAAAATTATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAgg                                                                                                       <  2:132111/145‑1 (MQ=255)                                                   ttATAGGGCTTTTAAAATTCTTAATAACGCGTCAGTTTTTTATTCCCTGATACCGATCCCGCCCATCTTTGGGGGATATCACATAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTg                                                                                                   <  2:655433/149‑1 (MQ=255)                                                    tataAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                                >  2:248833/1‑151 (MQ=255)                                                     ataGGGCTTTTTAAAATCGTGAAACCGCGTAAGAGTTTTATTCCCGAATCCCGATTCCGCCCAGCTATGGCGGATACCACTAAATCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGa                                                                                                  <  1:242099/148‑1 (MQ=255)                                                     ataAGGCTTTTTAAAATCCTGATACCGCGTCAGATTTTTATTCCCTGATCCCGATTCCGCCCAGTTATTGCGGATATCACGTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTa                                                                                               <  1:22949/151‑1 (MQ=255)                                                      taAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                                >  2:243873/1‑149 (MQ=255)                                                      taAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGAtt                                                                                                >  2:30912/1‑149 (MQ=255)                                                        aGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGc                                                                                             >  2:723776/1‑150 (MQ=255)                                                          gctttttAAATTCTTAATAACGCGTCATTTTTTTTTTCCCTGATCCCGTTTCCGTCCAGCTATGGCGAATCTCACGTAACCCAATATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAg                                                                                              <  2:665068/144‑1 (MQ=255)                                                           cGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaa                                                                                          >  2:600882/1‑150 (MQ=255)                                                           cGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaa                                                                                          >  2:600758/1‑150 (MQ=255)                                                            tttttAAATTCCTAATACCGCGTAAGTTTTTTATTCCCTAATCCCGTTCCCCCCCAGTTTTTGGGAATATCATTTAACCCAATATAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                         <  2:470998/149‑1 (MQ=255)                                                            gTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCaaaa                                                                                         >  1:620220/1‑150 (MQ=255)                                                             ttttCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAat                                                                                       >  2:650614/1‑151 (MQ=255)                                                               ttCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAatat                                                                                     >  1:828292/1‑151 (MQ=255)                                                                         gATACCGCGTCAGATTTTTATTCCCTGATCCACTCCCCGCCCAGCTATTGCGAATATCACGTAATCCCAAATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCg                                                                            <  1:290725/150‑1 (MQ=255)                                                                          aaagcgctTCAGTTTTTTTTTCCCTAATCCCGATTCCGCCCAGTTTTTGCGAATATCACGTAACCCCATATAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGt                                                                           <  1:49515/146‑1 (MQ=255)                                                                          aaaccgctTCAGATTTTTTTTCCCTGATACCGATCCCGCCCAGCTTTGGCGAATATCACATAACCCAATAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCAtt                                                                               <  2:722905/142‑1 (MQ=255)                                                                          aTACCGCGTCATATTTTTTTTCCCTGATCCCGATCCCGTCCAGCTATTGCGGATATCACTTACCCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGt                                                                           <  2:598749/150‑1 (MQ=255)                                                                          aTAACGCGTCAGATTTTTTTTACCTGATCCCGATCCCGCCCAGCTTTGGCGGATACCACGTAACCCAAAAAAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCg                                                                            <  1:908674/149‑1 (MQ=255)                                                                          aTAACGCGTCAGATTTTTATTCCCTGATCCCGTTCCCGCCCACTTTTGGCGGATATCACGTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTcc                                                                             <  1:305033/148‑1 (MQ=255)                                                                                        ttttttTCCCTAATACCGTTCCCGTCCAGCTATTGGGGATATCACTTAATCCCATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCggg                                                              <  1:262876/149‑1 (MQ=255)                                                                                        ttttATTCCCTGATCCCGTTCCCGTCCAGCTATTGCGGATATCACTTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgg                                                               <  1:799975/148‑1 (MQ=255)                                                                                          ttattaCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCgggg                                                             >  1:210451/1‑148 (MQ=255)                                                                                           tattcCCTAATCCCAATCCCGCCCAGTTATTGCGATTATAACGAAACCCAATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGgaga                                                          <  2:357294/150‑1 (MQ=255)                                                                                           tattaCCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAACCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGgaga                                                          <  2:973239/150‑1 (MQ=255)                                                                                                ccTGATCCCTATCCCGTCCAGCTATTGCGAATATCCCGTAATCCAATATAAAGAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATAtt                                                      <  1:786012/149‑1 (MQ=255)                                                                                                ccTGATCCCGATCCCGCCCAGCTATTGCGAATATCAATTAACCCAATAAAAAAAATATCCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGatat                                                       <  2:1023003/148‑1 (MQ=255)                                                                                                    aTACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGtt                                                >  2:1013968/1‑151 (MQ=255)                                                                                                     tACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGttt                                               >  1:738116/1‑151 (MQ=255)                                                                                                       ccgttcccGTCCAGCTTTTGCGGATACCCCTTAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTg                                              <  1:1001063/144‑1 (MQ=255)                                                                                                        cGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGa                                             >  2:511714/1‑150 (MQ=255)                                                                                                         gATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGAc                                            >  1:206169/1‑150 (MQ=255)                                                                                                             ccGTCCAGCTATTGGGGATACCACTTAACCCCATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATa                                       <  2:443749/151‑1 (MQ=255)                                                                                                              cGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGAt                                        >  2:316189/1‑149 (MQ=255)                                                                                                                tCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGtt                                    >  2:739998/1‑151 (MQ=255)                                                                                                                tCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGtt                                    >  1:310048/1‑151 (MQ=255)                                                                                                                    gttATTGCGGATATCACGTAACCCAATATAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgt                                  <  1:293199/147‑1 (MQ=255)                                                                                                                        tggggaatatCACATAACCCAAAAAAAAAAATTTTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatag                             <  2:298256/144‑1 (MQ=255)                                                                                                                        tggCGAATCTCACTAAACCCAATATAAAAAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcata                              <  1:37378/147‑1 (MQ=255)                                                                                                                            ggATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTgtcatagtcat                         >  2:64043/1‑150 (MQ=255)                                                                                                                                  cagttAACCCCATAAAAAAATTATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCAc                  <  1:939940/147‑1 (MQ=255)                                                                                                                                     gTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGcc               >  1:172525/1‑151 (MQ=255)                                                                                                                                       aaTCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCgg             >  1:250206/1‑151 (MQ=255)                                                                                                                                        aTCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCggg            >  1:294655/1‑151 (MQ=255)                                                                                                                                             atatCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCAtt        >  1:949643/1‑150 (MQ=255)                                                                                                                                                  aaGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGac  >  1:831005/1‑151 (MQ=255)                                                                                                                                                  aaGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGac  >  2:816757/1‑151 (MQ=255)                                                                                                                                                  aaGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCAGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGac  >  2:377289/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CAATAACGTTCTTACCGATTTTCTGCTGCAATCCCATTGCCAGTTCATCGTTATAAGGCGTTTTCAAATCCTGATAACGCGTCAGAGTTTTATTACCTGATACCGATTCCGTCCAGCTATTGCGGATATCACGTAATCCCATATCAAGAATATTCCCGCCATAGTAACGGTTATAACCTGCTGTAATCATTGAGGTTTGATTAGCAAAAATATCCCATTCCGTCATAAAGCGCGGGGAGATATTGTGGTTTGACAGATAGTTGTCATAGTCATACCGCACGCCGGGCATTAATGAC  >  NC_000913/1575933‑1576228

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑