Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,233	0	G	A	51.5%	38.5 / 91.0	68	pseudogene (60/75 nt)	ycgH	pseudogene; putative ATP‑binding component of a transport system
Reads supporting (aligned to +/- strand):  ref base G (20/13);  new base A (0/35);  total (20/48)
Fisher's exact test for biased strand distribution p-value = 6.98e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.79e-10

TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGA  >  NC_000913/1222086‑1222383                                                                                                                                                       |                                                                                                                                                       tgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGaaa                                                                                                                                                     <  1:716107/151‑1 (MQ=255) tgAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGaaa                                                                                                                                                     <  2:382699/151‑1 (MQ=255)   aaaGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAata                                                                                                                                                   <  1:461836/151‑1 (MQ=255)       tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATaa                                                                                                                                                <  1:633555/150‑1 (MQ=255)       tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATaa                                                                                                                                                <  2:160019/150‑1 (MQ=255)       tCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAg                                                                                                                                               <  1:1274443/151‑1 (MQ=255)         gAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGtt                                                                                                                                             >  2:702256/1‑151 (MQ=255)          aaCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGtt                                                                                                                                             >  2:1014388/1‑150 (MQ=255)          aaCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTc                                                                                                                                            <  1:826695/151‑1 (MQ=255)          aaCAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTc                                                                                                                                            <  2:860990/151‑1 (MQ=255)            cAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGtt                                                                                                                                             >  2:658346/1‑148 (MQ=255)             aaGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGa                                                                                                                                         >  1:1460974/1‑151 (MQ=255)                aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATc                                                                                                                                       <  2:559447/150‑1 (MQ=255)                aTGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATc                                                                                                                                       <  1:1505539/150‑1 (MQ=255)                    ggCGAGTAGTGTGGGGGAAATTAACCTTGGCGTAAATGG‑AAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGTGTGCAGCTGGGTGATAATGGCACCAAGGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGAGCccgc                                                                                                                                   <  1:844916/150‑1 (MQ=255)                                   gaaaaTTAAACTGGGCTAAAAGGC‑CACCCTAAACCCGGCGGCCAGGGTCGGGGGAA‑TGTGGGCGTGCGGCGGGGTAAAAATGCCTCCAACCCCCCCGCGGTA‑ATGGTGGCCTTAAAATAAAAGTTCTAACCCCGCCGTTAGCTAAAAAAc                                                                                                                    <  1:988235/148‑1 (MQ=255)                                           aaCTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTa                                                                                                           >  1:781370/1‑151 (MQ=255)                                              tgggCGTAAATGG‑CACCCGGAATCCAGCGGCCAGCGTGTGGGGAG‑TGTGGGGGTGCAGCTGGTTGATAATGGCTACAAGGACCCCGAGGTG‑AGGGGGGGCCTAAAATAAAATTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                        <  1:1709479/149‑1 (MQ=255)                                               tGGGGTAAAAGGCCACCCTAAACCCAGCGGCCGGCGTGGGGGAAA‑GGGGGGGGGGCAGCTGGGTAATAAGGGCTAAAATGACCCCGAAGTG‑AGGGGGGGCTTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTc                                                                                                        <  2:91942/151‑1 (MQ=255)                                                  ggaaaaTGG‑AAACAGGAACCCAGGGGCCGGCGTGGGGGGAA‑TGGGGGGGTGCAGTTGGGTGATAATGGCTAAAAGGCCCCCGCAGGG‑AGGGGGGGCCTAAAATAAAAGTTCTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                     <  2:1306256/147‑1 (MQ=255)                                                  cttAAAGGG‑AAACCTAAACCCAGGGGCCAGGGTCTGGGAAA‑TGTGGGGGGGCGGCTGGGAAATATTGGCTACAATGAC‑CCGGAGGGAATGGTGGGCCTAAAATATAGTTTCGACCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                     <  1:1298714/148‑1 (MQ=255)                                                  cgaaaaTGG‑CACCCAGAACCCAGCGGCCAGCGTGGGGGAAA‑TGTGGGCGTGCAGGTGGGTAAAAATGGCAACAAGGCCCCCGCAGTG‑AGGGTGGCCCTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                     <  1:1013861/147‑1 (MQ=255)                                                      aaTGC‑CAACCGGAACCCGGCGGCCGGGGTCGGGGGAA‑TGGGGGGGGGCAGCTGGGTAATAAGGCCTAAAATGACCCCGAAGTG‑ATGGGGGGCTTAAAAAATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcgc                                                                                                <  1:1310974/151‑1 (MQ=255)                                                       aTGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGActcg                                                                                               >  1:779800/1‑151 (MQ=255)                                                            cacctgcaacccAGCGGCCGGCGTGGGGGAAG‑GGGGGGCGGGCATCTGGTTTATAAGGGCTAAAATGCCCCGGAGGTG‑TTGGGGGGCCTAAAATAAAAGTCCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                           <  1:1706515/141‑1 (MQ=255)                                                            cacccTGAACCCGGCGGCCAGCGTCGGGGGAA‑TGGGGGGGTGCAGCTGGGTAATAATGGCTAAAAGGACCCCGCGGTG‑ATGGGGGGCCTAAAAAAAAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTAc                                                                                           <  2:1155468/148‑1 (MQ=255)                                                               ccgaaacccAGGGGCCAGCGTGGGGGGAA‑GGGGGGGGGGCAGGGGGGTAAAAATGGCTAAAATAACCCCGCAGTG‑AGGGGGGCCCTAAAATATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACata                                                                                        <  2:1613869/144‑1 (MQ=255)                                                                cTGAAACCGGCGGCCGGTGGTGGGGGAA‑TGGGGGGGTGGGGGTGGGTGATAATGGCTAAAATGACCCCGGAGTG‑AGGGGGGGCAAAAAAAAAAGTTTCTACCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                       <  2:1129132/151‑1 (MQ=255)                                                                 taaacccggcgcCCAGCGTGGGGAGAA‑GGTGGGGGGGCGGCTGGGTGATAAGGGCTAAAAGGACACCGCAGTG‑TGGGGGGGCTAAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                      <  2:1136874/143‑1 (MQ=255)                                                                 aaaacccgggggCCAGCGGAGGGGAAA‑GGGGGGGGGGCAGGGGGGGGAAAATGGCTAAAAGGACCCCGCAGGG‑AGGGGGGGCCAAAAATAAAATTTCTAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                      <  2:826678/141‑1 (MQ=255)                                                                  gAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGGTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGa                                                                                     >  1:663641/1‑151 (MQ=255)                                                                  aacgccAGCGGCAAGCGTCGGGGGAA‑TGGGGGCGTGCAGCTGGGTAATAAGGGCAAAAATGACCCCGAAGTG‑ATGGGGGGCTTAAAATAAATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                      <  1:482341/146‑1 (MQ=255)                                                                     ccccgcgcCCGGCGTCTGGGAAA‑GGGGGGGGTGCGGGGGGGTAAAAAGGGCTACAAGGACCCCGAGGGG‑AGGGTGGGCTAAAAATAAAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAat                                                                                  <  1:1069286/147‑1 (MQ=255)                                                                             ccggggTCGGGGGAG‑GGTGGGGGGGCAGGTGGGTAATAGTGGCAAAGATGCCACCGGGGTG‑AGGGGGGGCCAAAAATAGAGTTACTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                           <  1:877595/145‑1 (MQ=255)                                                                             cccagggtgggggAATTGGGGGGGTGCGGTTGGGTAAAAATGCCTAAAATGACCCCGAAGTG‑ATGGTGGGCCTAAAATATAAGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                           <  1:658058/142‑1 (MQ=255)                                                                             ccAGCGTGTGGGGAT‑TGGGGGCGAGAAGGGGGATGATACGGACTACAAGCACCCCGTAGTG‑TTGGTGCTCATAAAAAAAAATTTTGAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttt                                                                           <  1:1032768/150‑1 (MQ=255)                                                                                 cGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTg                                                                       >  1:1693521/1‑150 (MQ=255)                                                                                   tCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCg                                                                    >  1:1039599/1‑151 (MQ=255)                                                                                         gAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCa                                                              >  1:1389168/1‑151 (MQ=255)                                                                                           a‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                            >  1:849056/1‑151 (MQ=255)                                                                                                gggtgggCAGCGGGGTAAAAAGGCCAACAAGCACCCCGGAGTG‑TGGGGGGGCTTAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCAt                                                         <  1:942744/144‑1 (MQ=255)                                                                                                  gggTGCAGGGGGGTGATAAGGCCTAAAATGACCCCGCAGTG‑ATGGTGGGCCTAAAATAAAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                       <  2:955330/148‑1 (MQ=255)                                                                                                  gCGTGCAGTTGGGTAATATTGCTTCAAATGCCCCCTCAGTG‑GTGGTGGGCCTAAAATATATGTCCTGACCCCGCCGTTAGCTAAAAACCCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATgg                                                       <  1:716237/150‑1 (MQ=255)                                                                                                    gTGCAGCGGGGTAATAAGGGCTAAAATGCCACCGAGGTG‑ATGGGGGGCCAAAAATATAAGTTCGAATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                     <  1:186418/150‑1 (MQ=255)                                                                                                          gtGGGTTATAATGGCTACAATGACCCCGCGGTG‑GTGGTGGGCTTAAATTATAAGTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatc                                                <  2:45739/148‑1 (MQ=255)                                                                                                          cgggttataaagggCAACAAGGACCCCGGATTG‑ATGGGGGGCAGAAAATATAAGTTCGGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGtcatca                                               <  1:261247/138‑1 (MQ=255)                                                                                                           tGGGTGATAATGGCTAAAATGCCCCCGGAGTG‑ATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAAt                                             <  2:556852/151‑1 (MQ=255)                                                                                                              gtaaTAATGGCTAAAATGACCCCGCAGTG‑AGGGTGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTg                                           <  1:485318/147‑1 (MQ=255)                                                                                                               taaaaaTGGCTAAAATGACCCCCCGGTG‑ATGGTGGCCTTAAAATATAAGTTTTGCCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGc                                          <  2:1734178/146‑1 (MQ=255)                                                                                                                aaaaaTGGCTACAATAACCCCGAAGTG‑ATGGGGGGCCTAAAAAATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGcc                                         <  2:1441381/147‑1 (MQ=255)                                                                                                                 aTAATGGCTACAATCCCCCAGAGGTG‑AGGGGGGGCTTGATATAAATGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCg                                        <  1:670208/150‑1 (MQ=255)                                                                                                                      ggCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                                  >  1:1392031/1‑151 (MQ=255)                                                                                                                      ggCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCACGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCa                                  >  1:1207684/1‑151 (MQ=255)                                                                                                                        caaaaaTGACCCCGAAGGG‑AGGGGGGGCCAAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                                <  2:977672/147‑1 (MQ=255)                                                                                                                        aaaCAATGACCCCGAAGTG‑ATGGGGGGCCTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAAc                                <  1:1077251/149‑1 (MQ=255)                                                                                                                           cAATGCCCCCGAGGTG‑ATGGGGGGCCGAAAATATAAGTTTTAACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGc                             <  2:1631717/151‑1 (MQ=255)                                                                                                                                cccccGAAGTG‑AGGGGGGGCCTAAAATATAATTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTg                        <  2:52143/148‑1 (MQ=255)                                                                                                                                   ccGCAGTG‑ATGGTGGGCCGGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCTCGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCgg                     >  1:713469/1‑151 (MQ=255)                                                                                                                                      cactG‑ATGGAGGGCCTGGAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                   >  1:1607358/4‑150 (MQ=255)                                                                                                                                      cAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGt                   >  1:580526/1‑150 (MQ=255)                                                                                                                                        gggggtgggggGCCTGAATTATAAGTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTc                  <  2:1590682/145‑1 (MQ=255)                                                                                                                                        gTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATg                                                                                      >  2:1084750/1‑81 (MQ=255)                                                                                                                                             tggtggGCCTAAAATAAAAGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAg             <  1:1591262/150‑1 (MQ=255)                                                                                                                                                tggGCTTAAAATAAAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAgcgc          <  2:224073/150‑1 (MQ=255)                                                                                                                                                   gCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCTCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaaa       >  1:1163713/1‑150 (MQ=255)                                                                                                                                                   gCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCaaa       >  2:1232590/1‑150 (MQ=255)                                                                                                                                                    ccTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAc      >  2:157271/1‑150 (MQ=255)                                                                                                                                                    ccTAAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAAcc     <  2:1243130/151‑1 (MQ=255)                                                                                                                                                       gAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGa  >  1:255017/1‑151 (MQ=255)                                                                                                                                                       |                                                                                                                                                       TGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGGAAATTAAACTTGGCGTAAATGG‑CAACCTGAATCCAGCGGCCAGCGTCTGGGGAA‑TGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGACACCGCAGTG‑ATGGTGGGCCTGAAATATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGA  >  NC_000913/1222086‑1222383

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