Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,287,719	0	T	A	48.0%	74.4 / 58.5	77	intergenic (‑348/‑45)	acs/nrfA	acetyl‑CoA synthetase/nitrite reductase, formate‑dependent, cytochrome
Reads supporting (aligned to +/- strand):  ref base T (34/6);  new base A (0/37);  total (34/43)
Fisher's exact test for biased strand distribution p-value = 7.97e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACC  >  NC_000913/4287578‑4287867                                                                                                                                               |                                                                                                                                                     aaGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGc                                                                                                                                             >  2:826883/1‑151 (MQ=255)   gATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc                                                                                                                                            >  1:725242/1‑150 (MQ=255)   gATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGcc                                                                                                                                            >  1:577073/1‑150 (MQ=255)      aCTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTg                                                                                                                                        >  1:1072461/1‑151 (MQ=255)         gattAACTCTAAGGTGTTATTTTACATGCCTTTCCAATTGTTTAAAGACAACATTCCCGGGGTGTTTTTTTGTA‑ACAAATAGGGGCAAAAAATTTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGca                                                                                                                                      <  2:486648/147‑1 (MQ=255)             aaCTCTAAAGTGGTATTTTACATGCACTTAAAATTGATTAAAGACAACATTCACAGGGTGGTTATTTGTT‑ACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGcaaca                                                                                                                                   <  2:3661/149‑1 (MQ=255)              aCTCTAAAGTGTTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGGGGTTTTTTGTT‑ACACATGGGGGCGAGCAATGTCATGACAGTGTGGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATg                                                                                                                                <  1:412524/151‑1 (MQ=255)              aCTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGAATTGTCATGACAGGGTGGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATg                                                                                                                                <  1:134426/151‑1 (MQ=255)                 ctAAAGCGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGgcg                                                                                                                             >  2:292745/1‑151 (MQ=255)                   aaaGTGGTTTTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTTTCGTAAAAAACAATAAAAACCGCCATTGCAACAATGgc                                                                                                                              <  2:593369/148‑1 (MQ=255)                        ggTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                      >  2:727726/1‑151 (MQ=255)                        ggTATTTTACATGCACTTACAATTGATAAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGAAAGGTCATGACAGTGTAGGTGGGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCg                                                                                                                      <  1:149439/151‑1 (MQ=255)                              ttACATGCACTTAAATTTGATTAAAGAAAACATTCACGGTGTGGTTATTTGTT‑ACACATAGGGGCAAGCAATGTCAGGACGGTGTAGGTAGGTTTACGTTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGa                                                                                                                 <  2:98419/150‑1 (MQ=255)                                     cACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCt                                                                                                          >  1:177284/1‑150 (MQ=255)                                        ttACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTAt                                                                                                        <  1:249932/149‑1 (MQ=255)                                         tACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGAc                                                                                                     >  2:684409/1‑151 (MQ=255)                                          aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGa                                                                                                      >  2:301858/1‑149 (MQ=255)                                          aCAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGa                                                                                                      <  1:301914/149‑1 (MQ=255)                                           cAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                                    >  1:305766/1‑150 (MQ=255)                                            aaTTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACa                                                                                                    >  2:919738/1‑149 (MQ=255)                                            aaTTGATAAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGAGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGa                                                                                                      <  2:863480/147‑1 (MQ=255)                                                 aTTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATa                                                                                              >  2:359910/1‑150 (MQ=255)                                                   tAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                            >  2:922085/1‑150 (MQ=255)                                                   tAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaa                                                                                            >  1:1111138/1‑150 (MQ=255)                                                    aaaaACAACATTCACATTGTGTTTATTTGTT‑ACACATAGGGGCGACAAATGTCATGACAGTGTAGGTGCGGTAACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaa                                                                                             <  2:244892/148‑1 (MQ=255)                                                    aaaGACAACATTCACAGTGTGTTTATTTGTT‑ACACATGGGGGCGAGAAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATaaaaa                                                                                          <  1:680161/151‑1 (MQ=255)                                                          aacaTTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACGGTGTAGGTGCGGTTACTGTCGTGAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATaaa                                                                                      <  1:849978/149‑1 (MQ=255)                                                             aTTCACAGTTTGTTTATTTGTT‑ACACATGGGGGCGAGAAATGTCATGACAGTGTAGGTGCGTTTCCTGTCGTAAAAAAAAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgca                                                                                  <  1:589214/150‑1 (MQ=255)                                                             aTTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcac                                                                                 >  1:148688/1‑151 (MQ=255)                                                               tCACAGTGTGTTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAacgcacg                                                                                <  1:128352/150‑1 (MQ=255)                                                                 acgggtggtttATTTGTT‑ACACAAAGGGGCGGGCAGTGTCAAGGCAGTGTGGGTGCGGTTACTGTCGTGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg                                                                             <  2:411910/142‑1 (MQ=255)                                                                  caGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCg                                                                             >  2:328640/1‑150 (MQ=255)                                                                  caGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                            >  2:266493/1‑151 (MQ=255)                                                                   aGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGt                                                                            >  1:807931/1‑150 (MQ=255)                                                                          ttATTTTTT‑CCACAGAGGGGCGAAAAAGGACAAGACAGTTTAGGGGCGTTTACTGTCGAAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCtt                                                                       <  2:684696/148‑1 (MQ=255)                                                                           tttttGTTAACAAAAAGGGGCGAGAATTGCAATAAGGGTGGGGGGGCGGTTGCTGGAGAGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCa                                                                     <  2:301021/148‑1 (MQ=255)                                                                           tATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGc                                                                   >  2:653499/1‑151 (MQ=255)                                                                               ttttaacaaaagggggCGAGAATTGTCAGGACGGTGTGGTTGGGTTTCCTGTCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCtt                                                                 <  1:1084636/138‑1 (MQ=255)                                                                                 tt‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTg                                                             >  2:91890/1‑151 (MQ=255)                                                                                    acacATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAt                                                           >  2:70821/1‑151 (MQ=255)                                                                                    acacATAGGGGCGAGAAATGTAATACCAGGGTAGGTGGTTTTACTGTCGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTAtt                                                              <  2:959566/148‑1 (MQ=255)                                                                                     cacaGAGGGGCGAGAATGGTCAGGACAGTGAAGGGGAGGTTATTGTTGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                          <  1:585150/151‑1 (MQ=255)                                                                                      ataggggggggAGAAATGACAGGAGAGTGGAGGTGCTGTTAGTGTCGGGAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                          <  2:194806/145‑1 (MQ=255)                                                                                      acaTGGGGGCGAGAATTGAAATGACAGTGAGGGTGGGTTTACTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                          <  1:844724/150‑1 (MQ=255)                                                                                        aTAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCt                                                       >  2:424931/1‑151 (MQ=255)                                                                                        aTAGGGGCGAGAAATGTAAGAACAGTGTGGGTGCGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTcc                                                        <  2:1080033/150‑1 (MQ=255)                                                                                           ggggggAGAAATGTCATGAAGGGGTGGGTGGGTTTATTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCttt                                                     <  1:459933/150‑1 (MQ=255)                                                                                           gggcggaaaaaTGTAACGACAGTGTAGGTGGGGTTACTGCCGGAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGAtt                                                          <  2:134418/136‑1 (MQ=255)                                                                                             ggggagaAATGTCATGACAGTGTAGGTGGTGTTACTGTGGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                  <  2:104304/148‑1 (MQ=255)                                                                                             ggCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                  >  2:210223/1‑151 (MQ=255)                                                                                             ggCGAGAAAGTTCATGACAGGGTAGGTGGGTTAACTGCGGAAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCtttttt                                                  <  2:745188/151‑1 (MQ=255)                                                                                                   tgctttcccgattttgtgGGGGCGTTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                <  1:1053207/133‑1 (MQ=255)                                                                                                    aatttaataacggtgtAGGTGCGTTTACTTCCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttcttttt                                            <  2:486817/138‑1 (MQ=255)                                                                                                    aaTGTCATGACGGTGTGGGGGGGTTTCTTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttt                                             <  1:664893/149‑1 (MQ=255)                                                                                                    aaTGTCATGACAGTGTGGGTGGGTTTCCTGTCGTAAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttttc                                           <  1:580185/151‑1 (MQ=255)                                                                                                    aaTGTCATGACAGTGTAGGTGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttctttt                                             <  2:992204/149‑1 (MQ=255)                                                                                                      tGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGATTATTGATTCCTTTTTTCTTTTTCAc                                         >  1:1048103/1‑151 (MQ=255)                                                                                                       gTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACt                                        >  2:103912/1‑151 (MQ=255)                                                                                                         cATGACAGTGTGGGTGCGTTTCCTGCCGTAAAAAACAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACtt                                       <  1:392332/150‑1 (MQ=255)                                                                                                          aTGACAGTGTGGGTGGGTTTACTGTCGTAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTc                                      <  1:827095/150‑1 (MQ=255)                                                                                                             aCAGTGTAGGTGGGGTTACTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGt                                   <  2:145173/150‑1 (MQ=255)                                                                                                              cAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTc                                 >  1:339548/1‑151 (MQ=255)                                                                                                                gtgtGGGGGCGTTTACTGTCGGAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTg                                    <  1:507273/146‑1 (MQ=255)                                                                                                                 tgtAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTtttttct                                                >  2:1052931/1‑133 (MQ=255)                                                                                                                 tgggggTGCGGTTTCTGTCGTGAAAAAAAAAAAAACCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCa                                <  1:691782/145‑1 (MQ=255)                                                                                                                   taggggggtttACTGTCGAAAAAAACAAAAAAAACCGCCATTGCAACAATGGCGCAATTCTGTTGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCAc                               <  2:261512/139‑1 (MQ=255)                                                                                                                     ggTGGGGTTACTGTCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGt                                   <  2:105901/142‑1 (MQ=255)                                                                                                                            ttACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg                   >  2:811247/1‑151 (MQ=255)                                                                                                                            ttACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg                   >  1:775731/1‑151 (MQ=255)                                                                                                                            ttACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg                   >  1:60207/1‑151 (MQ=255)                                                                                                                            ttACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg                   >  1:1065768/1‑151 (MQ=255)                                                                                                                              aCTGTCGAAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACg                   <  2:520184/149‑1 (MQ=255)                                                                                                                                tGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCCGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctg               >  2:342213/1‑151 (MQ=255)                                                                                                                                  tCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgct             >  2:948555/1‑151 (MQ=255)                                                                                                                                  tCGTAAAAAAAAAAAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct                <  1:999730/148‑1 (MQ=255)                                                                                                                                     taaaaaaCAAAAAAAACCGCCATTGAAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgctgc              <  1:743619/145‑1 (MQ=255)                                                                                                                                       aaaaaaaaaaaaaaCCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGgct                <  1:727907/143‑1 (MQ=255)                                                                                                                                             aaTAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAAcc  >  1:332048/1‑151 (MQ=255)                                                                                                                                               |                                                                                                                                                     AAGATACTGACTAACTCTAAAGTGGTATTTTACATGCACTTACAATTGATTAAAGACAACATTCACAGTGTGGTTATTTGTT‑ACACATAGGGGCGAGCAATGTCATGACAGTGTAGGTGCGGTTACTGTCGTGAAAAACAATAAAAACCGCCATTGCAACAATGGCGCAATTCGGATGAAGCCCCTATGACAAGGATAAAAATAAACGCACGCCGTATCTTCAGCTTATTGATTCCTTTTTTCTTTTTCACTTCTGTTCACGCTGAACAAACGGCTGCTCCCGCAAAACC  >  NC_000913/4287578‑4287867

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