Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 17,436 C→T intergenic (+426/‑53) 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,4360CT90.0% 100.2 / 2.3 40intergenic (+426/‑53)sokC/nhaAsmall regulatory RNA antitoxin SokC/Na(+):H(+) antiporter NhaA
Reads supporting (aligned to +/- strand):  ref base C (1/3);  new base T (19/17);  total (20/20)
Fisher's exact test for biased strand distribution p-value = 6.05e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.41e-01

CAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAACGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAA  >  NC_000913/17193‑17678
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
cAAGCTGGATTATTTTTGAAATATTGCCCTACCAATTACCGCCGCCTGCCAAGAAATTAATATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCggg                                                                                                                                                                                                                                             <  1:223063/251‑1 (MQ=255)
aaaGCTGGATTATTTTGTAAAAATTTTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTGCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCATGGATTATCTGGGCTTAGGGGCACTCGTGAGCGCTTACAGTCGTCAAAAAACGCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCgg                                                                                                                                                                                                                                              <  2:303611/250‑1 (MQ=255)
        attattTTTGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAAt                                                                                                                                                                                                                                      <  1:341175/251‑1 (MQ=255)
              tttGAAATATTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGTTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGAAAAATAGTaa                                                                                                                                                                                                                                 >  1:293822/1‑250 (MQ=255)
                tGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAACGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTaaaaa                                                                                                                                                                                                                              <  1:190485/250‑1 (MQ=255)
                      atTGTCCTAACAAATAACGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGACACTCGTGAGCGCTTACAGCCGTCAAAATACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGATATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAGACTAGTAAACACGATCt                                                                                                                                                                                                                        >  1:244306/1‑251 (MQ=255)
                               acaaatcatgCCGACTGACAAGAAATTACTTATTACCCTTCCTAATTAACCCCTCACAAATACTCACCTTAAGCTATGATTTTCTGGGCTTAGGCTCACTAGTGACCGCTTACAGCCGTCAAAAACATCATCTCACCGCTGACGGCGCAAATTCTTCAATGGATCGTAAAAAAAGAATTATTCCTACACTATAATCTGATTTTAACAATGATTTATGCGGGGTAAACTAGTAAAAAAGCTCTATCCACCTg                                                                                                                                                                                                               >  2:236474/10‑251 (MQ=255)
                               aaaaataacGCCGACTGACAAGAAATTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCAACTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTg                                                                                                                                                                                                               <  2:427193/243‑1 (MQ=255)
                                      acGCCGACTGACAAGAAAGTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAgaga                                                                                                                                                                                                        <  1:483827/250‑1 (MQ=255)
                                                  aagaaaTTAATTATTACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGa                                                                                                                                                                                            <  2:37814/248‑1 (MQ=255)
                                                    aaaaTTAATTATTACTTTACCTAATTAATCCCTCAGAAAACCTCACCTTAAGCTATTATTATCTGGGCTTAGGGTCACTCGTGATCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCATATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCTGGGTAAAATAATAAAAACGATCTATTCCCCTGAAAGAGAAATAAAAAGTGaaa                                                                                                                                                                                          <  2:62093/250‑1 (MQ=255)
                                                           attattACTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGATATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAACCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTCACGATGATTTGTGCGGGGTAAATTAGTAAAAACGATCTTTTCACCTGAAAGAGAATaaaaacgacaaacatatgc                                                                                                                                                                                  >  2:444916/1‑237 (MQ=255)
                                                                 aCTTTTCCTAATTAATCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATACAAAGTGAAACATCTGCATCGAt                                                                                                                                                                             >  1:270606/1‑251 (MQ=255)
                                                                          aattaatCCCTCAAAAATCCTCACCTTAAGCTATGATTAACTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCCCCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATGTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCa                                                                                                                                                                    <  2:99928/251‑1 (MQ=255)
                                                                            ttaatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGt                                                                                                                                                                  >  2:262957/1‑251 (MQ=255)
                                                                              aatCCCTCAGAAATCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTATTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAACTAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACCTCTGCATCGATTCTTTAGCAGTgg                                                                                                                                                                >  2:123799/1‑250 (MQ=255)
                                                                                     cagaaaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTc                                                                                                                                                          <  1:406917/246‑1 (MQ=255)
                                                                                         aaTCCTTACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCTTGAGCGCTTACAGCCGTAAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCTCCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAg                                                                                                                                                     <  2:401456/251‑1 (MQ=255)
                                                                                         aaTCCTCACCTTAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGCGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCggg                                                                                                                                                       >  1:217932/1‑249 (MQ=255)
                                                                                                  cctAAGCTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCGATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCGCGGGAGGCATTATTc                                                                                                                                            <  2:302391/249‑1 (MQ=255)
                                                                                                        cTATGATTATCTGGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATc                                                                                                                                       >  2:367861/1‑250 (MQ=255)
                                                                                                                           aGGGTCACTCGTGAGCGCTTACAGCCGTCTAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCg                                                                                                                     >  1:93143/1‑250 (MQ=255)
                                                                                                                           aGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAAACTCATCTCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCgg                                                                                                                    >  1:103665/1‑250 (MQ=255)
                                                                                                                                                    cggcaaaaaaCTCATCGCACCGCTGATGGCGCAAATTCTTCAATGGATCGTAAAAAACGAGTTATTCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAATGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAAc                                                                                           <  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)
                                                                                                                                                                                                                        tCCTACACTATAATCTGATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCTCCACCAGTGGATGGTATCCAGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGGTAACTCGAAATcc                          >  2:264562/1‑250 (MQ=255)
                                                                                                                                                                                                                                        gATTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGACATAAACAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGCTTATGGCCACTAGCGGAGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCCGCTCCGCGTTGGTTCACTCGAAATCCACACAAACATGCTGtt          >  2:256397/1‑251 (MQ=255)
                                                                                                                                                                                                                                         aTTTTAACGATGATTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGACATCAACAAAAACATGCTGtt          >  2:12789/1‑250 (MQ=255)
                                                                                                                                                                                                                                                 gatgatTTGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGACTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAATAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGGGACGCCGGTTCAGCTCCGGGGTGGTTCACTCGAAATCAACAAAAAACTGCTTTTATGGATaa  >  2:344266/1‑250 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                  
CAAGCTGGATTATTTTTGAAATATTGGCCTAACAAGCATCGCCGACTGACAACAAATTAATTATTACTTTTCCTAATTAATCCCTCAGGAATCCTCACCTTAAGCTATGATTATCTAGGCTTAGGGTCACTCGTGAGCGCTTACAGCCGTCAAAAACGCATCTCACCGCTGATGGCGCAAATTCTTCAATAGCTCGTAAAAAACGAATTATTCCTACACTATAATCTGATTTTAACGATGATTCGTGCGGGGTAAAATAGTAAAAACGATCTATTCACCTGAAAGAGAAATAAAAAGTGAAACATCTGCATCGATTCTTTAGCAGTGATGCCTCGGGAGGCATTATTCTTATCATTGCCGCTATCCTGGCGATGATTATGGCCAACAGCGGCGCAACCAGTGGATGGTATCACGACTTTCTGGAGACGCCGGTTCAGCTCCGGGTTGGTTCACTCGAAATCAACAAAAACATGCTGTTATGGATAA  >  NC_000913/17193‑17678

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