Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 246,747 G→T S12S (TCG→TCT)  yafL → NlpC/P60 family protein YafL

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913246,7470GT97.8% 141.5 / NA 45S12S (TCG→TCT) yafLNlpC/P60 family protein YafL
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base T (27/17);  total (27/18)
Fisher's exact test for biased strand distribution p-value = 4.00e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGGGGTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAATGGTTTCCTGCTCCGCTC  >  NC_000913/246503‑246990
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tGGGGTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGAtt                                                                                                                                                                                                                                                 >  2:320466/1‑250 (MQ=255)
   ggTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTg                                                                                                                                                                                                                                             >  2:686517/1‑251 (MQ=255)
        ccccTGGGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATTACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGAt                                                                                                                                                                                                                                        <  2:40079/249‑1 (MQ=255)
          cctttgCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGAttt                                                                                                                                                                                                                                      <  2:144684/247‑1 (MQ=255)
                       tttattCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGCGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTAGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGttttt                                                                                                                                                                                                                        >  2:688988/1‑250 (MQ=255)
                            tCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCa                                                                                                                                                                                                                    <  1:403766/250‑1 (MQ=255)
                                cTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGGGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTCTTGACTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGc                                                                                                                                                                                                               >  2:17658/1‑251 (MQ=255)
                                  tGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGcc                                                                                                                                                                                                              >  1:3159/1‑250 (MQ=255)
                                     aTTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAgcg                                                                                                                                                                                                          >  1:297127/1‑251 (MQ=255)
                                                    aGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACata                                                                                                                                                                                           <  1:599497/251‑1 (MQ=255)
                                                            cGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTc                                                                                                                                                                                   <  1:594521/251‑1 (MQ=255)
                                                             gACAATCAATGAGTAATTTAATAAACGGATGTTTTTACATAATAAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAAGGATCCTGGATACAC‑TGCAACATTAAAGCATGACCAGCAATTATTAAGTATGTCC‑TTGCCGTCAATTCCATACTTTGTATTGTCTGGCTTACTGTTTATT‑GTATTACCAATTTCTTAATTTGCCAGCGCCACCACATCATATAtttctttct                                                                                                                                                                                  >  2:273077/1‑250 (MQ=255)
                                                                                          tGTTTTTACATAATTAACTGTAATATCTATATTAATTTTATACTAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCAACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAAc                                                                                                                                                     >  1:582367/1‑251 (MQ=255)
                                                                                           gTTTTTATATAATTAATTTAAATATATATTTTAATTTTATAA‑AATGTGTTTATTAATTTTGTTCAGAAACTACCCGGGATACAC‑TGCAAAAGTAAAGCATGCACAGCCTTTTTTGGGTATGTCT‑TTGCCGTCGATTCCATCCTTTGTATTGTCGGGTTTTCTTTTGATTTGT‑TTGCCTTTTTTTTCTTTTGCCAGCGCCGCCACACCCCATATCTCTTTCAGCTACACCTCCCGCCAGCGGAGGCAAAAcc                                                                                                                                                    <  2:402533/251‑1 (MQ=255)
                                                                                           gTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACcaca                                                                                                                                                                                                  >  2:536780/1‑205 (MQ=255)
                                                                                           gTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACcaca                                                                                                                                                                                                  <  1:536780/205‑1 (MQ=255)
                                                                                                  cATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGcgcg                                                                                                                                             >  1:653565/1‑251 (MQ=255)
                                                                                                    aaaTTAACTGAAATATATATTTTAATTTTATAATAATG‑GTTTATTAAGTTTGTACAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGtt                                                                                                                                           <  1:17658/250‑1 (MQ=255)
                                                                                                            taaaataaaatttttAAATTTATAAAAATG‑GTTTATTAATTTTGTACAAAAAAGA‑CCGGGAAACCCATTCACCATAAAACCATGCCCAGCCATTATTGAGTATTTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTaaa                                                                                                                                    <  2:3159/237‑1 (MQ=255)
                                                                                                                             ttttATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCt                                                                                                                  >  1:403905/1‑251 (MQ=255)
                                                                                                                                      aatG‑GTTTATTAAGTTTGTTCAGAACCGATCCGGGATACAC‑TGCAACATTAATGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCCCCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCa                                                                                                         <  1:555346/251‑1 (MQ=255)
                                                                                                                                              tattaAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACCTATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAACCAGTACCAAACTCATCTGAAAAAGCAGGCCAGc                                                                                                  >  2:639274/1‑251 (MQ=255)
                                                                                                                                                  aaGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTACGCGTTTATTAAAACAGTACCAAACACATCTGAAAAAGCAGGCCAGCtat                                                                                               >  2:82225/1‑250 (MQ=255)
                                                                                                                                                  aaGTTTGTTCAGAAACGACCCGGGATCCAC‑TGCAACATTAAAGCATGACCAGCCTTTATGTAGTATTTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCTTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCtat                                                                                               <  1:639274/250‑1 (MQ=255)
                                                                                                                                                               aaCGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACt                                                                                                                       >  2:11392/1‑213 (MQ=255)
                                                                                                                                                               aaCGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACt                                                                                                                       <  1:11392/213‑1 (MQ=255)
                                                                                                                                                                 cGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATg                                                                               >  1:644943/1‑251 (MQ=255)
                                                                                                                                                                       gggATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGAAATGCCGaa                                                                          >  2:488059/1‑250 (MQ=255)
                                                                                                                                                                             cac‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCCCCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCaaaa                                                                   <  1:686517/251‑1 (MQ=255)
                                                                                                                                                                              ac‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCaaa                                                                    >  1:100508/1‑249 (MQ=255)
                                                                                                                                                                               c‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAgg                                                                 >  2:552534/1‑251 (MQ=255)
                                                                                                                                                                                 tGCAACATTATAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTAAGCCGCCCGCCAGCGGCTGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCACTGCCGAAAGCAAAAggg                                                                >  2:395755/1‑251 (MQ=255)
                                                                                                                                                                                   cAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGgcg                                                              >  1:491375/1‑251 (MQ=255)
                                                                                                                                                                                   cAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGgcg                                                              >  1:212429/1‑251 (MQ=255)
                                                                                                                                                                                           aaccagccccAGCCATTAGTGAGTA‑GTCCCTTGCCGTCGTTTCCTTCCTTTGGATTGTCTGGATTACTTTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCCGCGCCACCACATCACATATCTCTTTCAGCTACCCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCa                                                      <  2:491375/243‑1 (MQ=255)
                                                                                                                                                                                                ggACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGcaca                                                 <  1:552534/250‑1 (MQ=255)
                                                                                                                                                                                                     aGCCATTATTTAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCgg                                            <  1:560309/251‑1 (MQ=255)
                                                                                                                                                                                                                 gTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACACCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGAAAAAGGgcg                                                              >  2:111258/1‑221 (MQ=255)
                                                                                                                                                                                                                 gTATGCCC‑TTGCCGTCGTTTCCATCCTTTGTATTGTCTGGATTCCTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCCCATATCTCTTTCAGCTACGCCCCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGgcg                                                              <  1:111258/221‑1 (MQ=255)
                                                                                                                                                                                                                   aTGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGCTGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTCTATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGAGaaaaa                               >  2:686456/1‑249 (MQ=255)
                                                                                                                                                                                                                                  gATTCCATCCTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAATGGt                >  1:120266/1‑251 (MQ=255)
                                                                                                                                                                                                                                           cttttttttGTCTGGATTACTGGTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATAACATACCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCa                           <  2:704834/224‑1 (MQ=255)
                                                                                                                                                                                                                                           cTTTGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCa                           >  1:704834/1‑231 (MQ=255)
                                                                                                                                                                                                                                              tGTATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTAAGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAATGGTTTCCTGCTCCGc    >  2:287745/1‑251 (MQ=255)
                                                                                                                                                                                                                                                tATTGTCTGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATCCAAAACCGTGCGCGTTTATTAACACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATTCCGAAAGAAAAAGAGCGCTACGCCAGCACAACCGGGCGCAGATAAAACAGCATCCAGTTTGGTTTCCTGATCCGCTc  >  2:45528/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TGGGGTATCTCCTGTGCTGTATATTTATTCAGCTTGAATTGTAGCTCAAATGAGCAACAGCGACAATCCATGAGTAATTTAATAAACGGATGTTTTTACATAATTAACTGAAATATATATGTTAATTTTATAATAATG‑GTTTATTAAGTTTGTTCAGAAACGATCCGGGATACAC‑TGCAACATTAAAGCATGACCAGCCATTATTGAGTATGTCC‑TTGCCGTCGATTCCATCCTTTGTATTGTCGGGATTACTGTTGATTTGT‑TTGCCGTTTTCTTCATTTGCCAGCGCCACCACATCACATATCTCTTTCAGCTACGCCGCCCGCCAGCGGATGCAAAACCGTGCGCGTTTATTAAAACAGTACCAAACTCATCTGAAAAAGCAGGCCAGCTATATTGTGGAAGGCAATGCCGAAAGCAAAAGGGCGCTACGCCAGCACAACCGGGAGCAGATAAAACAGCATCCAGAATGGTTTCCTGCTCCGCTC  >  NC_000913/246503‑246990

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