Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,222,233	0	G	A	28.9%	64.3 / 14.2	46	intergenic (+3032/+72)	ymgF/ymgD	inner membrane protein YmgF/PF16456 family protein YmgD
Reads supporting (aligned to +/- strand):  ref base G (27/5);  new base A (0/13);  total (27/19)
Fisher's exact test for biased strand distribution p-value = 1.17e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACTCAGGAGAAAAACAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGC  >  NC_000913/1221983‑1222472                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 aCTCAGGAGAAAAACAGCGGTCCGCACGGTGACTACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑Tg                                                                                                                                                                                                                                                  <  1:931850/251‑1 (MQ=255)            aaaCAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGACA‑TATAAGtt                                                                                                                                                                                                                                      >  1:513922/1‑251 (MQ=255)               cAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTg                                                                                                                                                                                                                                   <  1:246038/251‑1 (MQ=255)                         cACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGAAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAA‑TGATGGTGGGCC‑TGAAA‑TTTAAGTTCTGATCCCGCCGt                                                                                                                                                                                                                         >  2:359360/1‑251 (MQ=255)                                  gacTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑GGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTaaa                                                                                                                                                                                                                 >  2:708331/1‑250 (MQ=255)                                     tttGAACCTTTTGTTGAAATGAACTGGCTCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑CCCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGGTCCCGCCGTTAGCTaaaaacc                                                                                                                                                                                                             >  2:377389/1‑249 (MQ=255)                                               ttgttgAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAGTTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTAt                                                                                                                                                                                                   >  1:948468/1‑251 (MQ=255)                                                    aaaacgAACTGTAACCATAACAATAAAGTTTTTGCTGTCAAAA‑GTAGGGGTGAGAAAGTCGAAAAAGATGGGGCTAGTAATTTGGGGG‑AAA‑T‑AATGTTGGGCGTAAAGGA‑CATC‑CTGAATCCGGCGGCCAGCGCCTGGGATATGTGGGCGTGCAGCTGGGTGTTAATGACTACAATGAC‑CCCGCAG‑TGATGGTGGGCA‑TCAAA‑TATAAGTTTGGCGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATc                                                                                                                                                                                              <  1:218188/246‑1 (MQ=255)                                                     aaaTGAACTGGATCCATAAAAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCg                                                                                                                                                                                             >  2:698078/1‑251 (MQ=255)                                                       aTGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                                                                                                           >  1:726507/1‑251 (MQ=255)                                                        agaatgggATCCAAAACAGAAAAGTTTTTG‑TGGCCCAAATGAACGTTGTAAAAGCAGAAAAAGATGGGGCGAGTATTTTGGGAA‑AAA‑TTAAATTGGG‑CG‑AAATGGACAAC‑CTGAACCGAGGGGCCAGCGTCTGGGAATTGTGGGGGTGCAGCGGGGTAATAAGGGCTAAAATGAC‑AC‑ACAGGTGAGGGGGGGCC‑TAAAA‑TAAAAGTTCTGTTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAcg                                                                                                                                                                                          <  1:264641/245‑1 (MQ=255)                                                            atGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAATAAACCGCGTCGTATTCATCGACGCggg                                                                                                                                                                                      >  2:667997/2‑250 (MQ=255)                                                             tGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGt                                                                                                                                                                                      >  1:859934/1‑250 (MQ=255)                                                                  ccATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACatat                                                                                                                                                                                >  1:30599/1‑251 (MQ=255)                                                                         aGTAAAGTTTTTGCTGTCTAAA‑TGCAAGGTGTGAAAGTCAAAAAAGATGGGGCGATTAATTGGGGGG‑AAA‑TTAATCTTGG‑CGTAAAGGG‑CACC‑TTGAGCCCAGGGGCCGGTGTCTGGGGAATGGGGGCTTGCAGGGGGTTGATAATGGCTAAAATGAC‑CCCGCGG‑TGA‑GGTGGGGC‑TGAAAATATAATTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAatat                                                                                                                                                                         <  1:359360/251‑1 (MQ=255)                                                                          gTAAAGATTGTGCTGTCTCAA‑TGAAGGGGGTGAAAGACGAAAAAAATGGGGGGAG‑AATTTGGGGGTAAA‑TTAAACTTGG‑GTTAAAGGG‑CAAC‑CTGAATCGAGGGGCAAGCGCGGGGGGAATGGGGGTGTGGAGTGGGGTGAAAATGAC‑AAAATTACCACCGGAG‑TGTTGGGGGGCT‑TAAAA‑TATAATTTCGGCTCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAAtatt                                                                                                                                                                        <  1:775756/251‑1 (MQ=255)                                                                                 tgtggctgccaaaa‑TGAACGGGGTGAAAGGCAAACAAGAGGGGGCGATTAATTGGGGGA‑AAA‑TAAAACTTGG‑CGAAAATGG‑CACC‑CTGAACCCGGCGGCCAGCGTCTGGGGAAGGGGGTCGGGCAGCGGGGTAATAATGGCAACAATCAC‑CCCGAGG‑TGATGGGGGTCA‑TAAAA‑TATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTAttttt                                                                                                                                                                  <  1:492220/238‑1 (MQ=255)                                                                                    tGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATACAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGcc                                                                                                                                                              >  1:971247/1‑251 (MQ=255)                                                                                        gtccaaa‑TGAAGGGTGAGAAAGTCAAAAAAGATGGGGGAAGTAATTGGGGGG‑AAA‑TTAAACTGGG‑CGTAAAGGA‑CACC‑CTGAACCCAGCGGCCAGCGTCTGGGGAATGGGGGCGTGCAGCGGGGTAATAATGGCTAAAATGAC‑CCCGCAG‑TGATGGTGGGCC‑TAAAA‑TATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGAt                                                                                                                                                           <  1:25030/245‑1 (MQ=255)                                                                                         tatcAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGACTGTGGGCGTGAAGCTGGGAGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGCCGCGGTACATCTGAAATATTATTTTTTGCCGAt                                                                                                                                                           >  2:618040/3‑249 (MQ=255)                                                                                             aa‑TGAACGGGTTAAAAGTCAAAAAAGATGGGGCGA‑TTATTTGGGGG‑AAAATTAAATTGGG‑CAAAAAGGG‑CAAC‑CTGAATCCCGCGGCCAGCGTCTGGGAAGTGGGGGCGTGCAGCGGGGTGATAATGGCTAAAAGGAC‑CCCGCAG‑TGAGGGGGGTCC‑TAAAA‑TATAAGTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACg                                                                                                                                                     <  1:981297/251‑1 (MQ=255)                                                                                              a‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGc                                                                                                                                                    >  1:79037/1‑251 (MQ=255)                                                                                              a‑TGAACGGTGGAAAAGTCGAACAAGAGGGGGGGAGTAATTGGGGGG‑AAA‑TTAAACTTGG‑CGAAAAGGG‑CAAC‑CTGAACCCAGCGGCCAGCGTCTGGGGAATGTGGGCGGGAAGCGGGGTGATAAGGGCTAAAATGAC‑CCCGCAG‑TGAGGGGGGGCC‑TAAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGAc                                                                                                                                                                                           <  1:737878/212‑1 (MQ=255)                                                                                                     ggtggaaaaggcaaaaaaGAGGGGGGGAGTAATGTGGGGA‑AAA‑TTAAGTGGGG‑CGTAAAGGG‑CAAC‑CTAACCCCGGCGGCTGGCGTCTGGGAAATGTGGGGGTGAAGTGGGGTGAAAATGACAAAAATGAC‑CCCGCTG‑TGATGGGGGGCC‑AGAAA‑TATAATTTCTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGgtg                                                                                                                                              <  1:189804/235‑1 (MQ=255)                                                                                                            aaggcGAACAAGTTGGGGGGAGTATTTTGGGGA‑AAA‑TTAAACTTGG‑CG‑AAATTG‑GACCACTGAATCCGGCGGCCAGCGGCTGGGAAATGTGGGCGTGCAGCGGGGTGATAATGGCTAAAATGCC‑CCCCGAT‑TGAGGGGTGGCT‑TGAAA‑AATAATTTTTGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCaa                                                                                                                                       <  1:218183/247‑1 (MQ=255)                                                                                                                cGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTAGGTGATAATGGCTACAATGAC‑CCCGCAG‑TGATGGGGGGCC‑TGAAA‑TATAAGTTCTGGTCCCGCCGTTAGCTCAAAAACCGCGTCGTATTCATCGACGCGGTCCATATGAAATATTATTTTTTGCCGATAGCACGCATGGGGTCATCAATTGc                                                                                                                                   >  2:218477/1‑251 (MQ=255)                                                                                                                    cAAGATGAGGAGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGAAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTg                                                                                                                               >  2:912165/1‑251 (MQ=255)                                                                                                                     aaGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGGGTCATCAATTGCCGTg                                                                                                                               >  1:238756/1‑250 (MQ=255)                                                                                                                              cGAGTAATTTGGGGG‑CAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCCTGAAGCTGGGTGATAATGGCGACAATGAC‑ACCGCAG‑TGATGGAGAGCC‑TGAAA‑TATAAGTTCTGATACCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTACCGTGATCAACAGc                                                                                                                      >  1:978735/1‑250 (MQ=255)                                                                                                                                 gTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAttt                                                                                                                  >  1:103083/1‑251 (MQ=255)                                                                                                                                         gggg‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGGATCCAGCGGCCAGCGTCTGGGGAATGCGGGCGTGCAGCTGGGTGAAAATGGCTACAATGAC‑ACCGCAG‑TGAAGGTGGGCC‑AGAAA‑TATAAGATCTGATCCCGCCGTTAGCTAAAACACCGCGTCGTATTCATCAACGCGGTACATATGAAAAATTATTTTTTGCCGATAGCACGCATGGTGTCATCAACTGCCGTGATCAACAACATTTGCGGGTCt                                                                                                          >  1:902110/1‑251 (MQ=255)                                                                                                                                              aaaaatAAACTTGG‑CAAAAAGGG‑CAAC‑CGTAATCGGGCGGCCAGGGGCTGGGGAATGGGGGCGTGCCGGTGGGTGATAATGCCAACAAGAAC‑CCCGAAG‑GGATGGGGGGCC‑TAAAA‑TATAAGTTCTGTGCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTa                                                                                                       <  1:377389/246‑1 (MQ=255)                                                                                                                                                    aaaCTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTATGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGACATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCa                                                                                                  >  1:894343/1‑250 (MQ=255)                                                                                                                                                        gtGG‑CTTAAAGGG‑CACC‑CTGAAGCCAGCGGCCCGCTGCTGGGGAATGTGGGGGTGCAGGTGGTTGAAAATGGCTAAAATGAC‑CCCGCAG‑TGATGG‑GGGCCATAAAA‑TATAAGTTCGGACCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCt                                                                                             <  1:308877/250‑1 (MQ=255)                                                                                                                                                              atAAAGGG‑CAAC‑CTGACTCCGGCGGCCAGCGTCTGGGGAAGGTGGGCGTGCGGCTGGGTGATAAGGGAAACAATGAC‑CCCGAAG‑TGATGGGGGGCT‑TAAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTc                                                                                        <  1:698078/250‑1 (MQ=255)                                                                                                                                                                       ccacccagaacccAGCGGCCGGCGTTTGGGAAAGGTGGGCGTGCAGTTGGGTAAAAATGGCAAAAATGCC‑CCCGCGG‑TGATGGGGGCCC‑TAAAA‑TAAAAGTTCTGTCCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGttttt                                                                                 <  1:897975/240‑1 (MQ=255)                                                                                                                                                                        aaC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCCG‑TGATGGTGGGCC‑TGAAA‑TATAAATTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTCTTCCTCGACGCGGTACATATGAACTATTATTTTTTGCCGATAGc                                                                                                                                                        >  1:928164/1‑178 (MQ=255)                                                                                                                                                                        aaC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCt                                                                               >  1:904706/1‑251 (MQ=255)                                                                                                                                                                        aaC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCt                                                                               >  1:223847/1‑251 (MQ=255)                                                                                                                                                                           ccttaaTCCAGCGGCCAGGGTGTGGGGAATGGGGGCGTGGAGCGGGGTGAAAAAGGTTAAAAAGCC‑CCCGCTG‑GGGTGGGGGGCC‑TAAAA‑ATTAATTTTTCTTCCCGCCGTTAGCTTAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTAc                                                                             <  1:708331/247‑1 (MQ=255)                                                                                                                                                                             tGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTAcc                                                                            >  1:33322/1‑250 (MQ=255)                                                                                                                                                                                             cGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTc                                                           >  1:234393/1‑251 (MQ=255)                                                                                                                                                                                              gtcgggggATTGGGGGCGTGCAGCGGGGGGAGAAGGGCTACAATGAC‑CCCGGAG‑TGTTGGGGGGCC‑TGAAA‑TATAATTGCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCAGTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTcc                                                          <  1:412391/247‑1 (MQ=255)                                                                                                                                                                                                  ggggaatgggtgcttgcgtcggggTGATAATGGCTCAAATGAC‑ACCGTAG‑TGGTGGGGGGAC‑AAAAA‑AATATGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGgag                                                      <  1:539307/230‑1 (MQ=255)                                                                                                                                                                                                                     cTGGGGGATAATGGCTACAATGAC‑ACAGCAG‑TGATGGGGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGCAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCCATTGCGGGTCTTTAGCGCAAACCTGCTTCAGTTTTTCTACCACTTTGGCGCTAAGTTCCGGAGAATGCTCAATTTTTAAAt                                    >  2:546255/1‑250 (MQ=255)                                                                                                                                                                                                                          taaTAAGGGCTAAAATGAC‑CCCGCGG‑TGATGGTGGGCA‑AAAAA‑TATAAGTTCGGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACgg                               <  1:259261/248‑1 (MQ=255)                                                                                                                                                                                                                                                        aTGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCACAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCACTACCCATTACCGc  >  2:721164/1‑251 (MQ=255)                                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                 ACTCAGGAGAAAAACAGCGGTCCGCACGGTGACGACTTTGAACCTTTTGTTGAAATGAACTGGATCCATAACAGTAAAGATTTTGCTGTCTCAA‑TGAACGGTGTGAAAGTCGAACAAGATGGGGCGAGTAATTTGGGGG‑AAA‑TTAAACTTGG‑CGTAAATGG‑CAAC‑CTGAATCCAGCGGCCAGCGTCTGGGGAATGTGGGCGTGCAGCTGGGTGATAATGGCTACAATGAC‑ACCGCAG‑TGATGGTGGGCC‑TGAAA‑TATAAGTTCTGATCCCGCCGTTAGCTAAAAAACCGCGTCGTATTCATCGACGCGGTACATATGAAATATTATTTTTTGCCGATAGCACGCATGGTGTCATCAATTGCCGTGATCAACAGCATTTGCGGGTCTTTAGCGCAAACCTGATTCAGTTTTTCTACCACTTTGGCGCTCAGTTCCGGAGATTGCTCAATTTTTAAATCACGGCTGGCAACGCTGGCATTACCCATTACCGC  >  NC_000913/1221983‑1222472

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