Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,224,0210AG37.1% 43.8 / 11.6 36intergenic (+572/+258)ymgJ/minEuncharacterized protein YmgJ/Z‑ring positioning protein MinE
Reads supporting (aligned to +/- strand):  ref base A (4/18);  new base G (13/0);  total (18/18)
Fisher's exact test for biased strand distribution p-value = 1.61e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ACGTGGAATATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAC  >  NC_000913/1223772‑1224271
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
aCGTGGAATATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGgag                                                                                                                                                                                                                                                                <  1:659430/251‑1 (MQ=255)
           aCCAGGCAGTTATTAAAGCCTCATTCAGCAGTTCGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGaat                                                                                                                                                                                                                                                      <  1:732703/250‑1 (MQ=255)
                      aTTAAAGACTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGTCGGTGTGGAATTCCCGTGGAACGCGGTGGCTGTTGTGAAATGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGTAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGAGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTg                                                                                                                                                                                                                                           <  2:1038849/249‑1 (MQ=255)
                                   tCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑GAACCTCATCCGATTGGTGGACAATGAATAATAAAA‑TTTTGGGG‑TTattaattt                                                                                                                                                                                                                           >  1:419464/1‑245 (MQ=37)
                                   tCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAAACTCATCCCACGGGGGGACAAAGAATAATTAAA‑ATGCGGGG‑GTatcattcc                                                                                                                                                                                                                           >  2:943688/1‑248 (MQ=37)
                                        agTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAAACTCATCCTATTGGGGAAAAAAGAAAAATGAAA‑ATGCGGGG‑GTATCATCTcccccg                                                                                                                                                                                                                      >  2:221614/1‑248 (MQ=37)
                                             gTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATc                                                                                                                                                                                                                 <  1:779555/251‑1 (MQ=255)
                                               ttAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGGGAAACC‑AATCCCATGTGGGAACAATAAAAAATGAAA‑ATGGGGGG‑GTattattctcccctatatta                                                                                                                                                                                                                >  2:1227009/1‑233 (MQ=25)
                                                 aGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCAt                                                                                                                                                                                                              <  1:943688/250‑1 (MQ=255)
                                                   cGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑GAAACTCATCCTATTGGGGAAAAATGAATAAAGAAA‑TTGTGGGG‑GTGTCATCTCCCAGTATAACCATACt                                                                                                                                                                                                           >  1:724112/1‑251 (MQ=37)
                                                             ggggtggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAAACCCCTCCTATGGGGGAAAAAAGAATAATGAAA‑TTGCGGGG‑GTATCATCTCCCCGTatatt                                                                                                                                                                                                                 >  1:466008/1‑234 (MQ=37)
                                                             ggggtggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACGGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATc                                                                                                                                                                                                                 <  2:466008/235‑1 (MQ=255)
                                                                  ggggtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑GAACCTCATCCTATGTGGGGACAATGACTAATGAAAATTCGGGGG‑TTTTC‑TCTCCCACTATATTCACTCTTACAATAAGGTTAtt                                                                                                                                                                                            >  1:48287/1‑251 (MQ=25)
                                                                     gtATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGgaga                                                                                                                                                                                                                                                               <  1:395477/183‑1 (MQ=255)
                                                                             aTGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACGGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑g                                                                                                                                                                                <  2:419464/251‑1 (MQ=255)
                                                                               gggtccGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGTTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCTCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATcc                                                                                                                                                                                                                <  2:790715/214‑1 (MQ=255)
                                                                                  gCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGGGAAACCCCCTCTTTGGGGGGAAAAAT‑AATAATGAAA‑TTTCGGGG‑GTTTTTTTTCCCCGAATATCCATACTAAAAAAAAGGTTAttttttctaacaagctc                                                                                                                                                                             >  2:213408/1‑236 (MQ=25)
                                                                                   ccgttgtgGAATCCCCGTTGAACGCGGTGGCTGGTGAGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTGCTGCGCATTTTTCGTTTACAGGATGTGGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCTCGCCACTTACGTCGGG‑AAACCTCATCCTATGTGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCATTATATCCATACTAACAATAAGGTTATTTACTCATCCA‑GGCATaa                                                                                                                                                                          <  2:1024493/247‑1 (MQ=255)
                                                                                               ccccGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCAGCCTATGGGAGAACAATGAATAATGaaa                                                                                                                                                                                                                                               <  1:286170/173‑1 (MQ=255)
                                                                                                                         aaCTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGTGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAAACTCCTCCTATGGGGGAACAAAGAATAATTAAA‑ATGCGGGG‑GTATCATCTCCCCCGATACCCCAAATAAACATAAAGTTTTTTACTCAACCC‑GGCCAAAAC‑ATTTTTTTTTGTTCGGGGGAACAGCCTT‑AAGgtggg                                                                                                                                  >  1:1220676/1‑249 (MQ=255)
                                                                                                                              gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑GAAACACATCCTCTGGGGGAACAAAGAATAATGAAA‑TTTCGGGG‑GTATTATTTCCCCATATATCCATAATAACAATAAAGTTATTTTTTCAACCC‑AGCACAAAC‑ATTTTTTTTTTTGTGGGGGTACAGCCTT‑TAGGGGGAAAggg                                                                                                                             >  2:661134/1‑251 (MQ=255)
                                                                                                                              gTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCACCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAggg                                                                                                                             <  1:795376/251‑1 (MQ=255)
                                                                                                                              gTCGTTCTGACAATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAACAAAA‑GCCCGCCAGTTACGGCGGG‑CAACCTCCCCCTTTGGGGAAAAAAAGAAATATTAAA‑ATGCGGGGGTTTTC‑TCTCCCCCTATCTCCATACTAACAATAAGGGTTTTTACTCACACA‑AGCCTAAACCATTTT‑TTTTTTGCGTGGGGACAATCTT‑AAGGGGTAAAggg                                                                                                                             >  2:795376/1‑251 (MQ=255)
                                                                                                                                cGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAAACCCATCCTATGGGGGAACAATGAATAATGAAA‑TTGGGGGG‑TTATCATTTCCCAGAATATCCATACTAAAAATAAGGTTATTTACTTAACCA‑AGCATAAAC‑ATTTTTTTTTTTGCGGGGGAACAGCCTT‑AAGGGGGAAAggggg                                                                                                                           >  1:252274/1‑251 (MQ=255)
                                                                                                                                   tCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAAAGCCCGCCAGTTACGGCGGG‑AAAACCCATACTATAGGAGAAAAATGAATAATTAAA‑TTGCGGGG‑GTATCATTTCCCAGTAAATCCATAATAACAATAAGGGTGTTTAATCCACCAAGGC‑TAAAA‑AATTTTTTTTTTGCGGGGGGAAAGCCTT‑TAGGGGTAAAggggg                                                                                                                           >  1:263501/1‑249 (MQ=255)
                                                                                                                                    cTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGGGAAAAAAGAAAAATGAAA‑TTGCGGGG‑GTATCATCTCCCAGCATATCCATACTAAAAATAAAGTTATTTTTTCCACCA‑CGCACAAAA‑ATTTTTTTTTTTTGGTGGGGACAGACTT‑TAGGGGGAAAggggggggg                                                                                                                       >  2:911806/1‑247 (MQ=255)
                                                                                                                                                      aCTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCa                                                                                                                                                                                  <  1:1175895/177‑1 (MQ=255)
                                                                                                                                                                      gTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACt                                                                                     <  1:221614/251‑1 (MQ=255)
                                                                                                                                                                           cGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTcgcc                                                                                 <  1:661134/250‑1 (MQ=255)
                                                                                                                                                                                   ggATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGGGAACAAAGAAAAATGAAA‑TTGCGGGG‑GTATCATCTCCCAGTATATCCATACTAACAATAAGGGTATTTACTCAACCA‑AGCATAAAC‑ATTTTTTTTTTTGCGTGGGAACAGCCTT‑AAGGGGTAAAGGGGGAGGGGGGAATAGCAATGAGGAGTATCAGCAAGAATACACGCCGCTTTAcc                                                                         >  1:1106467/1‑250 (MQ=255)
                                                                                                                                                                                     atggtGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACcaca                                                                      <  2:48287/247‑1 (MQ=255)
                                                                                                                                                                                                            gCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTTTCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGaaaaa                                               <  1:213408/251‑1 (MQ=255)
                                                                                                                                                                                                                      gCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCa                                      <  1:1149227/251‑1 (MQ=255)
                                                                                                                                                                                                                         cGCCAGTTACGGCGGG‑AAACCTCATCCTATTGGAGAACAATGAATAATGAAA‑TTACGGGG‑TTATCATCTCCCAGTATATTCATAATAACAATAAGGGTATTTACTCAACCA‑AGCATAAAC‑ATTTTTTTTTTTGCGTGGGAACAGCCTTTAAGG‑GGAAAGGGGGGGGGGGGAATAGCAATGAGGGGTATCAGCAAGAATAATCGCCGCTTTTCCACAACGGGGGTGAGAGGGGGGAAAAACTCaaaggag                                     >  1:14988/1‑244 (MQ=255)
                                                                                                                                                                                                                                       gg‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTg                     >  1:172620/1‑251 (MQ=255)
                                                                                                                                                                                                                                                cATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATaa              <  2:252274/250‑1 (MQ=255)
                                                                                                                                                                                                                                                   ccTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATgc          >  2:188314/1‑251 (MQ=255)
                                                                                                                                                                                                                                                      aTGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCt       <  2:1106467/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           agaACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAc  <  1:982058/251‑1 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                               
ACGTGGAATATACCAGGCAGCTATTAAAGCCTCATTCAGCAGTACGTTTAGCGGACATCTCGGGGTGGGGTATAGCCATGGTGCCGGTGTGGAATCCCCGTGGAACGCGGTGGCTGGTGTGAACTGGTCGTTCTGACCATCAACGATTAAACTGCGCTTCGGCGCAGTTTTCGTTTACAGGATGTTGAAAGGGAAAATTCTGGGGCAAAAAAA‑GCCCGCCAGTTACGGCGGG‑AAACCTCATCCTATGGGAGAACAATGAATAATGAAA‑TTGCGGGG‑TTATCATCTCCCAGTATATCCATACTAACAATAAGGTTATTTACTCAACCA‑GGCATAAAC‑ATTTTGTTTTGTGCGTGGGAACAGCCTT‑AAGGTGTAAAGGGGGAGGTGGAAATAGCAATGAGGAGTATCAGCAAGAATACTCGCCGCTTTACCACAACGTGGATGAGAGGGATGAAAAACTCAAGGCAGAGATAACTCTGCCTTGAAGATAAATGCGCTTTTAC  >  NC_000913/1223772‑1224271

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