Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 23,104 A→G V238V (GTA→GTG)  ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091323,1040AG86.0% 120.5 / 8.4 59V238V (GTA→GTG) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (3/5);  new base G (27/22);  total (31/28)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

CAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGTAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCG  >  NC_000913/22867‑23348
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
cAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccact                                                                                                                                                                                                                                             <  1:113525/250‑2 (MQ=255)
   aaCTGAAGCCCACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCAAATCTCGCTGGTGATCTGGACCACTACGc                                                                                                                                                                                                                                         >  1:94221/1‑251 (MQ=255)
            cAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGActc                                                                                                                                                                                                                                <  1:10627/251‑1 (MQ=255)
               catcatCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGCCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgc                                                                                                                                                                                                                             >  1:367100/1‑251 (MQ=255)
                  catcCGTGCGCAGGGAAAAATCATCGGCAACAGGCTCCTGCACAAGTGCGCGAAGCC‑GTTACACAGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTCTTACGCAAAAACTTCTCCGGCCATCGACGTTGCTTTTCAGGCGGTCGATTAGGATGCACTGAAAGCAAAATTTGTCGTAATCAACGTTAACGGCCCAATCTCGCTGGTTATCTGTACCACTACGCCGTGGACTctgcatg                                                                                                                                                                                                                          <  2:369210/251‑4 (MQ=255)
                         gcgcTGGGCAAAAGCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACcg                                                                                                                                                                                                                   <  1:319312/251‑1 (MQ=255)
                          cgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACcgc                                                                                                                                                                                                                  >  1:292481/1‑251 (MQ=255)
                             ttggCAAAATCATCGGCATCGGTCACCTGCACCAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTAGTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGACGATCGATCAGGCTGAACTGAAAGCAAAATTTGCCGTAAGCAACGCTACCGGCCCAATCACGCTGGTTATATGGAACACTACGCCGTGGACTCGCCCTACGAACCgcgc                                                                                                                                                                                                                >  2:288086/3‑250 (MQ=255)
                                gCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAAtc                                                                                                                                                                                                            >  1:84339/1‑251 (MQ=255)
                                     atcatcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTAt                                                                                                                                                                                                       <  1:380803/251‑1 (MQ=255)
                                            gCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGACTGCCAACCGAGCAATCTCTATTGCAca                                                                                                                                                                                                 >  2:78572/1‑249 (MQ=255)
                                              aaCGGTCACCTGCACAAAGGCGCGAAGCCGGTA‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAcctg                                                                                                                                                                                               >  1:295584/1‑248 (MQ=255)
                                                                                cACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCCCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCTCTGGTGCAGATCGACGGTCa                                                                                                                                                             <  1:176412/251‑1 (MQ=255)
                                                                                cACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTc                                                                                                                                                              <  1:381459/250‑1 (MQ=255)
                                                                                 aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAg                                                                                                                                                            <  1:242473/251‑1 (MQ=255)
                                                                                      tGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGt                                                                                                                                                       >  1:123931/1‑251 (MQ=255)
                                                                                       gcggtGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTg                                                                                                                                                      <  1:28281/247‑1 (MQ=255)
                                                                                        cGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGACCGTGa                                                                                                                                                     <  2:80197/251‑1 (MQ=255)
                                                                                          ttGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCGGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGAtt                                                                                                                                                   >  1:318516/1‑251 (MQ=255)
                                                                                          ttGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGAtt                                                                                                                                                   >  1:283002/1‑251 (MQ=255)
                                                                                              ctgacGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAAAACGTCTCCGTCCATCGACGGTGCTTTTCAGGCGGTCGATCCGGATGCACTGACAGCAACATTTGCCGTAAGCAACGTTAACGGCCCTATCTCGCTGGTGATCTGGACAACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTTTTGCACCATATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGcagctcttct                                                                                                                                                 >  2:270170/5‑240 (MQ=255)
                                                                                                gCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTAAGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCACGCAATCTCTATTTCACCTGATTTAGACTATGCGCTGGTGCAGAACGACGGGCAGGCCGTGGTTCTGGCg                                                                                                                                             >  1:193842/1‑251 (MQ=255)
                                                                                                        gcgcTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGACTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCACTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGa                                                                                                                                        <  1:78572/248‑1 (MQ=255)
                                                                                                          ccTGGAGGAAGCGGAATTTGAGTTTTACGACAAAACTTCTCCTGCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTTGCGAAAGAGCTgg                                                                                                                                   <  2:300165/250‑1 (MQ=255)
                                                                                                               cggaagcggcagTTGAGTATTACGACAAAACTTCTCCGTCCATAGACGTTGCTTTTCAGGCGGTCGATCAGGATGCAATGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGAACCTGATTTCGACTATGCGCTGGTGAAGATCGACGTCCAGGCCGTGATTCTGGCGAAAGATCTGGTTGa                                                                                                                               >  2:362690/1‑250 (MQ=255)
                                                                                                                   agcggaagTTGAGTATTACGACAAAACTTCTCCGTGCATCTACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCCCGCCGTTGACTCTGCCTGCCAACCGCGCAATCTCTATTGCAGCAGATTTCGACTATGCGCTGGTGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCg                                                                                                                          <  2:81508/251‑1 (MQ=255)
                                                                                                                      ggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTa                                                                                                                        >  2:193650/1‑250 (MQ=255)
                                                                                                                         agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                                                                                                    <  1:144319/251‑1 (MQ=255)
                                                                                                                          gTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCa                                                                                                                   <  1:113558/251‑1 (MQ=255)
                                                                                                                            tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGAGGTCGATCAGGCTGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGAGATCTGGACAACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATATCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGCCGGACAGCCCGTGATTCTGGCGAGAGATCTGGTTGGAAGCGGAGTGCAGc                                                                                                                 >  2:850/1‑251 (MQ=255)
                                                                                                                                          aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc‑c                                                                                                  >  1:230837/1‑251 (MQ=255)
                                                                                                                                          aaaCTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc‑c                                                                                                  <  1:70087/251‑1 (MQ=255)
                                                                                                                                                   cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACAcc                                                                                         >  1:137541/1‑251 (MQ=255)
                                                                                                                                                   cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACAcc                                                                                         <  1:73616/251‑1 (MQ=255)
                                                                                                                                                              gTTGATTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGAAACGTTAACGGCCCACTCTCGCTGGGGATCTGGACCACTACGACGTGGCCTCTGCCTGCGAACCTCGCAATCTCTATTGCACCTTATTTCGACTATGCGCTGGCGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTTAAAGAGTAATGCAGCTTCTCGGCGTAGC‑CGATTaccacaatttcgacac                                                                              >  2:25393/1‑237 (MQ=255)
                                                                                                                                                                tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTg                                                                                                      <  1:293802/226‑1 (MQ=255)
                                                                                                                                                                tGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGagatcggaagagcgacgggtcggg                                                                              >  2:293802/1‑227 (MQ=255)
                                                                                                                                                                    tttcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTaa                                                                         <  1:211986/247‑1 (MQ=255)
                                                                                                                                                                      tcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTa                                                                                             >  1:207447/2‑228 (MQ=255)
                                                                                                                                                                      tcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTa                                                                                             <  2:207447/227‑1 (MQ=255)
                                                                                                                                                                        aGGCGGTCGATCAGGGTGCACTGAAAGCAAAATTTGCCGTACGCAACGTTAACCGCCCAAGCTCGCTGGTGCTCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTGTGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTaaaatgt                                                                    <  1:346754/251‑4 (MQ=255)
                                                                                                                                                                            ggTCGCTCAGGATGCACTGAAAGCAAAATTTGCCGTACGCAACGTTAACGGCCCAATCTCGCTGTTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTACACCTGATTTCGACTCTGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCGTGCGACATATCTGTTTGAAGAAGGAATACAGAGGATCGGCGTAAC‑CGATTCCACCATtctcaacacgtaaaaaggagcgg                                                                >  2:225540/2‑232 (MQ=255)
                                                                                                                                                                             gTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑Cggg                                                              >  1:6508/1‑250 (MQ=255)
                                                                                                                                                                               cGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑c                                                           <  1:126546/251‑1 (MQ=255)
                                                                                                                                                                                 aTCAGGATGCACTGAAAGCAAAATTTGCCGGAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAATATCTTGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑Ctt                                                         <  1:170099/251‑1 (MQ=255)
                                                                                                                                                                                    aGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTg                                                        >  1:359574/1‑249 (MQ=255)
                                                                                                                                                                                     ggATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGAAATGCAGCGGATCGGCGTGAC‑CGATTACACCATTTTCAGCACGGTAACAGGTG‑CGGAG‑CTTGAgt                                                     >  2:62424/1‑250 (MQ=255)
                                                                                                                                                                                          cACTGAAAGCAAGATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGGTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAACACGGTAAAAGGTG‑CGGAG‑GTTGAGgtgctg                                                >  1:259725/1‑251 (MQ=255)
                                                                                                                                                                                             tGAAAGCAAAATTTCCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTgcgc                                             >  1:311499/1‑251 (MQ=255)
                                                                                                                                                                                              gAAAGCAAAATTTGCCGGAAGCAACGGTAACGGCTCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACGCTGCCTGCGAACCGCGCAATCTCTATTGCACTTGATTTCGACGATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGGAATGCAGCGTATCGGCGTGAC‑CTATTACACCAGTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGGTGCTgcg                                              <  2:378096/249‑1 (MQ=255)
                                                                                                                                                                                                   caaaattttccgatatcaactgtAACGGCCCAATCTCGCTGGTGATCTGGACCTCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGAGTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGGTTTAca                                       <  2:157882/229‑2 (MQ=255)
                                                                                                                                                                                                           tccGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGCTTTACCCATCCGt                                <  2:168669/249‑1 (MQ=255)
                                                                                                                                                                                                              gTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGCTTTACCCATCCGTTTAt                            <  1:103580/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑AGATTAAACCATTCTCCACACGGAAAAAGGTGACGG‑G‑CTTGAGTTGCTTCGCTTTACCCATCCGTTTATGGGCtt                      >  2:273366/1‑251 (MQ=255)
                                                                                                                                                                                                                     aCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTAAAAGGTG‑CGGAG‑CTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTc                     <  1:368873/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAACGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCAGCACGGTACAAGGGG‑CGGAG‑ATTGAGTTGCTGCGCTTTACCCATCCGTTTATGGACTTCg                    >  1:42938/1‑251 (MQ=255)
                                                                                                                                                                                                                                  aaTCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGAGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCCTTCTCAGAACGGTAAAAGGAG‑CGGGG‑TTTGAGATGCCGCGCAGTACCCATCCGTTTATGGGCCTCGACGTACCGGcc         >  2:145304/1‑249 (MQ=255)
                                                                                                                                                                                                                                     ctcGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGATAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGGAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATcc     >  1:84614/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGAGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGTCGCAACGCG‑TTACAACATTCTCAGCACGGAAAAAGG‑G‑ACGAGATTTGTGTTGGTGGGCTTTACCCATCCGTTTATGAGCtttggcggtctggaactacacg  >  2:303397/1‑231 (MQ=255)
                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                       
CAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTT‑CACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAC‑CGATTACACCATTCTCGGCACGGTAAAAGGTG‑CGGAG‑CTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCG  >  NC_000913/22867‑23348

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