Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,718,257 A→C 14.6% S90R (AGT→AGG)  yfiF ← putative methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,718,2570AC14.6% 153.4 / 4.0 65S90R (AGT→AGG) yfiFputative methyltransferase
Reads supporting (aligned to +/- strand):  ref base A (30/23);  new base C (5/7);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 5.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCGTGCCTGACGCTTTTGTCAATTCCGCTTCATCCACCACATGGTACGCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTG  >  NC_000913/2718020‑2718491
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
ttCCTTGCCTTGCGCTTTTTTCAATTCCGCTTCATCCCCCACATGGTACGCTTTGCTGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCAc                                                                                                                                                                                                                                              <  2:933443/239‑1 (MQ=255)
     tGCCTGACGCTTTTGTCAATTCCGCTTCATCCACCACATGGTACGCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCCGGTTTCTTCCGCACGCTGACGACAGCAA‑ACCTCCCGCTCACAAAAAACTTTACCAccact                                                                                                                                                                                                                                          >  2:967747/1‑234 (MQ=255)
      gCCTGACGCTTTTGTCAATTCCGCTTCATCCACCACATGGTACGCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATg                                                                                                                                                                                                                                         <  1:678619/238‑1 (MQ=255)
          gACGCTTTTGTCAATTCCGCTTCATCCACCACATGGTACGCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACCCTGACGCC‑CCCACACTTCCGGTCCAATAAAACATTTACCACGGAAGcgcac                                                                                                                                                                                                                                    >  2:684615/1‑235 (MQ=255)
               ttttGTCCATTCCCCTTCATCCACCACATGGTACGGTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACCCTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGa                                                                                                                                                                                                                                <  1:751242/238‑1 (MQ=255)
               ttttGTCAATTCCGCTTCATCCACCACATGGTACGGTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAACCGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGAt                                                                                                                                                                                                                               <  1:926463/239‑1 (MQ=255)
                       aTTCCGCTTCATCCACCACATGGTACGCTTTGCGGTCTGCTCCCATCCAGCGCCACGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCATGCACGGACAATCTCTTCCGGACGGCTCTGGCACAGTGCCTGACACGCCTATTCGCCGCAGACGAGGGTTTCTTCTGCACGCTGACGAA‑GCAACACTTCCGGGTCAATATAACTTTTTCACCCGAAGCCACCggttatggacttt                                                                                                                                                                                                                       >  2:1070462/1‑227 (MQ=255)
                                            ggTACGCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑CCAACCCTTCCGGATCAATAAAACTTTTACCCCTGATGCCACCCGGGTCGGCCTTTTCCGGCgcaccactcaccgg                                                                                                                                                                                                   >  2:129133/1‑225 (MQ=255)
                                              tACGCTTTGCGGTTTGCTGCCATCCAGCGAAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑gcg                                                                                                                                                                                                <  1:927699/238‑1 (MQ=255)
                                                 gCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCCCGCTGACGCC‑GCAACACTTCCGGATCAATAAAACTTTTACAACTCATGCCACCGTGAACGGCCTTTTCCGGCGTCCCATCACCCG‑Ccgcgcg                                                                                                                                                                                            >  1:1065765/1‑239 (MQ=255)
                                                 gCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGCCCTTTTCCGGCGTCTCATCACCCG‑Gcgcgcg                                                                                                                                                                                            >  1:49771/1‑239 (MQ=255)
                                                    ttGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGCCATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGAACCCTGACGCC‑CAAACACTTCGGTATCAATAAAACTTTTACCCCCGTAGCCACCGTGACAGCCCTTTTCCGGCTCGTCACAACCCG‑ACGCTCgggc                                                                                                                                                                                         >  2:648912/1‑238 (MQ=255)
                                                    ttGCGGTTTGCGGCCATCCAGCGCAAAGCTACTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGCCATGCATTTTAGCCGTAGACGCGGGGTTCTTCCGCACGCTGACGCC‑GCAACACGTCAGGGTCAATAAAACGTTTACCCCTGATGCCACCGTTACCCGCCTTTTGCGGCTACTCCTTACCCC‑GGGCACGGG‑a                                                                                                                                                                                        >  2:1066815/1‑239 (MQ=255)
                                                     tGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCgggga                                                                                                                                                                                        >  1:880769/1‑237 (MQ=255)
                                                             gctgcCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCCCGCTGACGCC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCCCCTCACCCG‑GCCCGCGGG‑AAACCTTgg                                                                                                                                                                                >  2:760000/1‑237 (MQ=255)
                                                              ctgcCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTACGCCCGCTGACGCC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGCTCCGCCTTTTCCGGCGCCCCATCACACC‑CCCCCCGGG‑AAACCGGgcg                                                                                                                                                                               >  2:702948/1‑238 (MQ=255)
                                                               tgcCATCCAGCCGAAAGGTTATTTAAAGCGCGGCGTTACACCCTGGATAAACCAGGCGCGAACAATCGCGTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATGTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTgcgc                                                                                                                                                                              <  1:120991/238‑1 (MQ=255)
                                                                                 ttcttTAAAACGCGGGGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTTCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACt                                                                                                                                                           <  1:702948/239‑1 (MQ=255)
                                                                                     ttAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCg                                                                                                                                                       <  1:681705/239‑1 (MQ=255)
                                                                                         aaCGCGGCTTTACCCCTGGGATAACCCAGGCCGGAAACATCGCTTCCGGACGGCTTTGGAAAAGTCCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCCCGCTGACGAC‑GCAACACTGCCGGATCAATAAACCTTTTACCACTGATGCCGCCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCa                                                                                                                                                   <  2:880769/239‑1 (MQ=255)
                                                                                         aaCGCGGCGTTACCCTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCa                                                                                                                                                   <  1:614710/239‑1 (MQ=255)
                                                                                          aCGGGGCGTTCCCCTTTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCAc                                                                                                                                                  <  1:252680/239‑1 (MQ=255)
                                                                                          aCGCGGCGTTACCCTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCAc                                                                                                                                                  <  1:967747/239‑1 (MQ=255)
                                                                                              ggcgTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCCCGCTGACGAC‑GCACCACTTCCGGATCAATAAAATTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGCCCCATTACCCG‑GCGCGCGGG‑AAACCGTCCCCCCCTGCGAACCCTCCCCCTTACggcccggggg                                                                                                                                              >  2:614710/1‑231 (MQ=255)
                                                                                                    acacTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGAct                                                                                                                                        <  1:133632/239‑1 (MQ=255)
                                                                                                        tctGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCCCTTGCGGTCACGGGGCTGACTCTGt                                                                                                                                    >  1:278677/1‑239 (MQ=255)
                                                                                                          tGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGCTGCCACCGGGATCGGCCTTTTCCGGCGCCTATTTACCCGTGCG‑GCGGGAAAACCGGGC‑CCACGGCGAACCTTCCCCCTTTCCGTCACcgggtgaacacctatg                                                                                                                                  >  2:435902/1‑227 (MQ=255)
                                                                                                          tGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGttt                                                                                                                                  >  1:170619/1‑239 (MQ=255)
                                                                                                            gATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGGGATCGGCCTTTTCCGGCGGCCCATTAACCA‑GCGCGCGGG‑AAACCGTGCGCCCCGGCCAACTCTCCCACTTGCGGTCACGGGGcggcccccgcttc                                                                                                                                 >  2:150308/1‑226 (MQ=255)
                                                                                                                  ccAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGGGCGCCACGGCGAATCTTCCCCCTTCCGGTCACGGGGCGGACTCTGTTTTTCATccc                                                                                                                          >  1:681834/1‑237 (MQ=255)
                                                                                                                      gcgcGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGCCCCTGTTTGTCATCccgggc                                                                                                                      >  1:241147/1‑239 (MQ=255)
                                                                                                                      gcgcGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCccggcc                                                                                                                      >  1:855577/1‑237 (MQ=255)
                                                                                                                      gcgcGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCCCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACGCTGTTTGTCATAccgggc                                                                                                                      >  1:246921/1‑239 (MQ=255)
                                                                                                                      gcgcGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGGCTCATCACCCG‑GCGCGCGGG‑AAACCGGGCGCCACGGCGAATCTTCCCCCTTGCGGTCACGGGCGCGCCTCTGTTTgcactaccgggg                                                                                                                      >  2:720771/1‑228 (MQ=255)
                                                                                                                             ccttCGCTTCCGGACGGCTAGGGAAAATTGCCTTGCATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGGTGGCGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCTCCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCGTCCCGCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGa                                                                                                               <  2:1065765/236‑1 (MQ=255)
                                                                                                                                                        gTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTAATCGCGGGCGGGGCGCCGGCCACCGCCCGCACgagggttt                                                                                       >  1:305856/1‑231 (MQ=255)
                                                                                                                                                        gTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCCGCTTGTCATCGCGGGCCGGGCGACGGACACCGCAAGACCGACAtttt                                                                                       >  1:985224/1‑236 (MQ=255)
                                                                                                                                                        gTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTCCGCCACGGCGGATCTTCCCCCTTGCGGTCACGGGGCTGACTCTGTTTTGCATCGCGGGGGGGGCGCGGGCCCACCCCAGACCGAGATTTTccg                                                                                    >  2:651423/1‑238 (MQ=255)
                                                                                                                                                        gTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGGGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGCCATCGCGGGCGGGGTGACGGCCCCCGCCACGACGGAAtttt                                                                                       >  1:820040/1‑236 (MQ=255)
                                                                                                                                                           ccTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTTCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGCCAGCACGACATTTTCCT‑gg                                                                                 >  1:206028/1‑238 (MQ=255)
                                                                                                                                                                cATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGGC‑GCAACACTTCCGGATCAATAAAACTTTTACCACGGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGAcgccc                                                                           <  1:435902/239‑1 (MQ=255)
                                                                                                                                                                           gCCGTAGCCGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTcc                                                                <  1:720771/239‑1 (MQ=255)
                                                                                                                                                                              gTAGACGCGGGTTTCTTCCGCACGGCGACGAC‑GCAACCCTTCCGGATCAATAAAACTTTTACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCgg                                                             <  1:108813/239‑1 (MQ=255)
                                                                                                                                                                               tAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGGCTCTGTTTGTCATCCGGGGCGGGGGGACGGCCACCGCCACCCAGGACTTTTCCTGGGAGGCCAGCCCCC‑TTCCCGGt                                                            >  2:392546/1‑239 (MQ=255)
                                                                                                                                                                                            cttcCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                                               <  1:651423/239‑1 (MQ=255)
                                                                                                                                                                                                cGCACGCTGCCGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCCTCGCGGGCGTGGCGACGGCCACCGTCCGGACGAGATTTTCCT‑GGGCGCCCCCCCCCTTTCCCCGGCCCCGGGTCGTggttga                                           >  2:1062516/1‑235 (MQ=255)
                                                                                                                                                                                                 gcccGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCCCCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACg                                          <  1:129133/236‑1 (MQ=255)
                                                                                                                                                                                                        ggcgcc‑gCAACCCCTCCGGATCAAAAAAACTTTTCCCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTTCGGCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCa                                   <  2:855577/234‑1 (MQ=255)
                                                                                                                                                                                                           cgccccaccccTTCCGGATCAATAAACCTTTTACCCCTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                 <  1:684615/229‑1 (MQ=255)
                                                                                                                                                                                                           cgacccaccctttCCGGTTCAATAAAACTTTTCCCCCTGATGCCACCGTTATCGGCCTTTGCCGGCGTCTCATCCCCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAACCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                 <  2:241147/228‑1 (MQ=255)
                                                                                                                                                                                                           ccacagacatcctCCGGAATAATAAACCTTTTCCCCCGGATGCCCGCCTGTTCGGCCTTTTCCGGGGGTCCCTCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCCAAACTGCCCTCTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑TGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGa                                 <  2:753763/227‑1 (MQ=255)
                                                                                                                                                                                                                            ttCAATAAAACTTTTCCCCCTGATGCCCCCGTGATCGGCCTTTCCCGGCGGCTCATCCCCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGa                <  2:883390/238‑1 (MQ=255)
                                                                                                                                                                                                                            aTCAATAAAACGTTTCCCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCGCCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCgg                 <  1:904902/238‑1 (MQ=255)
                                                                                                                                                                                                                            aTCAAATAAACTTTTACCACGGATGCCCCCGTGAGCGGCCTGTTCCGGCGTCTCATCCCCCG‑GCGCGGGGG‑AACCCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGGTTGGCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCgg                 <  1:1062516/238‑1 (MQ=255)
                                                                                                                                                                                                                                aTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGAcgcccg                                                                          <  1:880940/177‑1 (MQ=255)
                                                                                                                                                                                                                                aTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGAcgcccg                                                                          >  2:880940/1‑177 (MQ=255)
                                                                                                                                                                                                                                aTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacata            >  1:883390/1‑239 (MQ=255)
                                                                                                                                                                                                                                aTAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGCTGTGGGTACGTTTATCAGAATCATCATCCCTGCGGacaa             >  2:598726/1‑237 (MQ=255)
                                                                                                                                                                                                                                 tAAAACTTTCACCCCTGGTGCCACCGTGTTCGGCCTTGTCCGGCGTCCCATCGCCCG‑GCGCGCGGG‑AAACCGTTGGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTTACTCTGTTTGCCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacatac           <  2:246921/239‑1 (MQ=255)
                                                                                                                                                                                                                                  aaaaCTTTTCCCAACGACGCCCCCTTGTCTGGCCTTTTCCGGCGTCCCCCCCCCCG‑GCACGCGGG‑AATCCGTCCCCCACGGCGAATCTTCCCACTTGCGGCCACGGGGCTGAATCTGTTTTGCAGCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGacataca          <  2:985224/239‑1 (MQ=255)
                                                                                                                                                                                                                                       ttttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTgg                                               <  1:563067/197‑1 (MQ=255)
                                                                                                                                                                                                                                         ttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGT‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGATTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGTTTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGGTTATCATAATCATCATCACTCCGGACATACATCACttt   >  1:722382/1‑239 (MQ=255)
                                                                                                                                                                                                                                         ttACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACCCGGGTTGTGGGTACGTTTATCCGAATCATCATCCCTCCGGACATACATCCCtt    >  1:176661/1‑238 (MQ=255)
                                                                                                                                                                                                                                           aCCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCCCTTGCGGTCCCGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGCCGGCCCCCGTCAGCACGCGATTTTCTT‑GTCCGCCCCCCCCCTTTCCCGGTACCCGGGTTTTGGTTCCGTTTTTCAGAACCCTCATCCCAGCGGCCCACCATCCCTTTg  >  2:623193/1‑238 (MQ=255)
                                                                                                                                                                                                                                           aCCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGCTCGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGGTGTGGGGACTTTTTTCAGAATCACCATCCCTGCGGCCAAAAATCACTTTg  >  2:618114/1‑238 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                             
TTCCGTGCCTGACGCTTTTGTCAATTCCGCTTCATCCACCACATGGTACGCTTTGCGGTTTGCTGCCATCCAGCGCAAAGCTTCTTTAAAACGCGGCGTTACACTCTGGATAAACCAGGCGCGAACAATCGCTTCCGGACGGCTCTGGAAAAGTGCCTGACATGCATTTTCGCCGTAGACGCGGGTTTCTTCCGCACGCTGACGAC‑GCAACACTTCCGGATCAATAAAACTTTTACCACTGATGCCACCGTGATCGGCCTTTTCCGGCGTCTCATCACCCG‑GCGCGCGGG‑AAACCGTGCGCCACGGCGAATCTTCCCACTTGCGGTCACGGGGCTGACTCTGTTTGTCATCGCGGGCGGGGCGACGGCCACCGTCAGCACGAGATTTTCCT‑GGACGCCCGCCCCCTTTCCCGGTACGCGGGTTGTGGGTACGTTTATCAGAATCATCATCACTGCGGACATACATCACTTTG  >  NC_000913/2718020‑2718491

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