Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,399 A→G intergenic (+389/‑90) sokC → / → nhaA small regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091317,3990AG88.1% 100.8 / 3.9 42intergenic (+389/‑90)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base A (2/3);  new base G (19/18);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

TTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTG  >  NC_000913/17166‑17639
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
ttATCTCTTTATCGCAATTATTGACGACAAGCTGGATTATTTTTGTAATCTTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGACGCAAATTCTTCAATGGCTCGTAAAAAACGAGTTATTCATACACtat                                                                                                                                                                                                                                   >  2:426814/1‑250 (MQ=255)
    ctctTTATCGCACATATTGACGACAAGCTGGATTATTTTTGAAATATTGTCATAACAAATAACGCCGCCTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCCCCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTAGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑CCTCATATCACCGCTGATGGCGCAAATTCTTCAATGGTTCGAAAACAAAGAGTTATTCCTACACTATCATCt                                                                                                                                                                                                                              >  2:378540/1‑251 (MQ=255)
          agcGCAAATATTGACGACAAGCTGGATTATTTGTGTAAGATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGTGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtgt                                                                                                                                                                                                                         <  2:204937/248‑3 (MQ=255)
          aTCGCAATTATTGACGACAAGATGGATTATTTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGAtttt                                                                                                                                                                                                                        >  2:323167/1‑251 (MQ=255)
                 ttattGACGACAAGCTGGATTATTGTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGACATCCTCACATTAAGCTATG‑ATTATCTGGACTTAGGGTCACTCGTGAGCGCTTACATCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACTAGTTATTCCTACACTATAATCTGATTTTAACgatg                                                                                                                                                                                                                 >  2:378945/1‑251 (MQ=255)
                           cAAGCTGGATTATTTTTGAAATATTGCCCTACCAATTACCGCCGCCTGCCAAGAAATTAAT‑ATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCggg                                                                                                                                                                                                      <  1:223063/251‑1 (MQ=255)
                           aaaGCTGGATTATTTTGTAAAAATTTTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTGCCTAATTAATCCCTCAGAAATCCTCACCTTAAGC‑ATGGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGTCGTC‑AA‑AA‑AAC‑GCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCgg                                                                                                                                                                                                       <  2:303611/250‑1 (MQ=255)
                                   attattTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAAt                                                                                                                                                                                               <  1:341175/251‑1 (MQ=255)
                                    ttattTTTGAAATCTTGTCCTAACAAATAACGCCGAATGAAAAGAACTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATA‑ATTATCTGGGATTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AACAA‑A‑C‑TCATCTCACCGCTGATGGCGCAACTTCTTCAATGGATCGTAAAAAACGTGTTATTCCTACACTATAACATGATTTTAACgatgctttattcggggaacaat                                                                                                                                                                                               >  1:365850/1‑232 (MQ=255)
                                         tttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGTTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGAAAAATAGTaa                                                                                                                                                                                          >  1:293822/1‑250 (MQ=255)
                                           tGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTaaaaa                                                                                                                                                                                       <  1:190485/250‑1 (MQ=255)
                                                 atTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGTC‑AA‑AATA‑C‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAGACTAGTAAACACGATCt                                                                                                                                                                                 >  1:244306/1‑251 (MQ=255)
                                                          acaaatcatgCCGACTGACAAGAAATTACTTATTACCCTTCCTAATTAACCCCTCACAAATACTCACCTTAAGCTATG‑ATTTTCTGGGCTTAGGCTCACTAGTGACCGCTTACAGCCGTC‑AA‑AA‑A‑CATCATCTCACCGCTGACGGCGCAAATTCTTCAATGGATCGTAAAAAAAGAATTATTCCTACACTATAATCTGATTTTAACAATGATTTATGCGGGGTAAACTAGTAAAAAAGCTCTATCCACCTg                                                                                                                                                                        >  2:236474/10‑251 (MQ=255)
                                                          aaaaataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCAACTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTg                                                                                                                                                                        <  2:427193/243‑1 (MQ=255)
                                                                 acGCCGACTGACAAGAAAGTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAgaga                                                                                                                                                                 <  1:483827/250‑1 (MQ=255)
                                                                             aagaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGa                                                                                                                                                     <  2:37814/248‑1 (MQ=255)
                                                                               aaaaTTAATTATTACTTTACCTAATTAATCCCTCAGAAAACCTCACCTTAAGCTATT‑ATTATCTGGGCTTAGGGTCACTCGTGATCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCATATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAATAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGaaa                                                                                                                                                   <  2:62093/250‑1 (MQ=255)
                                                                                      attattACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGATATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAACCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTCACGATGATTTGTGCGGGGTAAATTAGTAAAAACGATCTTTTCACCTGAAAGAGAA‑Taaaaacgacaaacatatgc                                                                                                                                           >  2:444916/1‑237 (MQ=255)
                                                                                            aCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATACAAAGTGAAACATCTGCATCGAt                                                                                                                                      >  1:270606/1‑251 (MQ=255)
                                                                                                     aattaatCCCTCAAAAATCCTCACCTTAAGCTATG‑ATTAACTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATGTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCa                                                                                                                             <  2:99928/251‑1 (MQ=255)
                                                                                                       ttaatCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGt                                                                                                                           >  2:262957/1‑251 (MQ=255)
                                                                                                         aatCCCTCAGAAATCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTATTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAACTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACCTCTGCATCGATTCTTTAGCAGTgg                                                                                                                         >  2:123799/1‑250 (MQ=255)
                                                                                                                cagaaaTCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                                                                   <  1:406917/246‑1 (MQ=255)
                                                                                                                    aaTCCTTACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCTTGAGCGCTTACAGCCGTA‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCTCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAg                                                                                                              <  2:401456/251‑1 (MQ=255)
                                                                                                                    aaTCCTCACCTTAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                >  1:217932/1‑249 (MQ=255)
                                                                                                                             cctAAGCTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCGCGGGAGGCATTATTc                                                                                                     <  2:302391/249‑1 (MQ=255)
                                                                                                                                   cTATG‑ATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                                                                                >  2:367861/1‑250 (MQ=255)
                                                                                                                                                      aGGGTCACTCGTGAGCGCTTACAGCCGTCTAA‑AA‑A‑C‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                                                                              >  1:93143/1‑250 (MQ=255)
                                                                                                                                                      aGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑AAC‑TCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgg                                                                             >  1:103665/1‑250 (MQ=255)
                                                                                                                                                                               cggcaaa‑aa‑a‑C‑TCATCGCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAATGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAAc                                                    <  1:209338/246‑1 (MQ=255)
                                                                                                                                                                                                  ctcACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGCAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGTTTATGGCCAATAGCGGCGAAACCAGGGGATGGTATc                                      >  1:437986/1‑250 (MQ=255)
                                                                                                                                                                                                    cACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCa                                     <  1:48024/249‑1 (MQ=255)
                                                                                                                                                                                                    cACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCGTTAGCAGTGATGCCTCGGGAGGCATTATGCATATCATTGCCGCGATCCGGGCGATGATTATGGCCAATAGCGGCGCAACCAGTCGATGGTATCAc                                    <  2:60793/250‑1 (MQ=255)
                                                                                                                                                                                                       cGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGAGCGGGGTAACATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCGGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACt                                >  1:497/1‑251 (MQ=255)
                                                                                                                                                                                                           gATGGCGCAAATTCTTCAATGGAGCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACttt                              <  2:93143/249‑1 (MQ=255)
                                                                                                                                                                                                                     aTTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCg                  <  1:323167/251‑1 (MQ=255)
                                                                                                                                                                                                                          tctgtggatCGTAAAAAACGAGTTATTCCTCCACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCCCCTTAAATAGAAATAAACAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGTAGGCATTATTCTTATCATTGCAGCTATCCTGGCGATGATTATGGCCATTATTGGTTCAATCATTGGATGGAATCACGACTTTCTGGAGACGCCGGTTcc             >  2:15718/9‑250 (MQ=255)
                                                                                                                                                                                                                            aatggatCGTAAAAAACGAGTTATTCCTACACTAGAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTAGCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGc           <  2:77699/245‑1 (MQ=255)
                                                                                                                                                                                                                            aaTAGCTCGTAAAAAACGAATGATTCCTACACTATAATCTGATTTTAACGATGAGTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCGGAAAGAGAAGTAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCTTCGGGAGTCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGc           <  1:301432/251‑1 (MQ=255)
                                                                                                                                                                                                                                     tAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAACTAGTAAAACCGATCTCTTCAACTGAAAGAGAAATAACAAGTGAAACCTCTGCATCGATTATTTAGCAGTGAAGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCACTAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGTAGACGCCGGTTCAGCTCCGggttg  >  1:36399/1‑251 (MQ=255)
                                                                                                                                                                                                                                             gAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCt          <  1:90029/235‑1 (MQ=255)
                                                                                                                                                                                                                                             gAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAAAAGCCGGTTCAGCt          >  2:90029/1‑235 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
TTATCTCTTCATCGCAATTATTGACGACAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATG‑ATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTC‑AA‑AA‑A‑C‑GCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTG  >  NC_000913/17166‑17639

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