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,1040AG95.6% 201.8 / ‑0.5 69V238V (GTA→GTG) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base G (29/36);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.22e-01

CATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGC  >  NC_000913/22858‑23350
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
cATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGAAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCt                                                                                                                                                                                                                                                     >  1:281383/1‑251 (MQ=255)
   ggaattCAAAACTGAAGCCAACATCAACCGCGCGCTGGCCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGACGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGAGCTGGa                                                                                                                                                                                                                                                  <  2:269656/247‑1 (MQ=255)
    gACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGTCAAAAGCATCGGCAACGGTCACCTGCCCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATGTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGac                                                                                                                                                                                                                                                 <  2:500961/251‑1 (MQ=255)
       ttCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                                                                                                                                                                                              >  2:217186/1‑251 (MQ=255)
        tCAAAACGGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACGTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGaccac                                                                                                                                                                                                                                              <  2:393554/250‑1 (MQ=255)
               tGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGt                                                                                                                                                                                                                                      >  1:8002/1‑251 (MQ=255)
                gAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTg                                                                                                                                                                                                                                     <  1:65444/251‑1 (MQ=255)
                  aGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                                                                                                    >  2:251561/1‑250 (MQ=255)
                               cgcgcgCTGGGCAAACTCACCGGTAACGGTCACCTGTGCAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAACCTGCTCCGTCCATCGACGTTG‑TTTCTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCTCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcctgcgaa                                                                                                                                                                                                                      <  1:20772/251‑4 (MQ=255)
                                   cgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTAAGCCGTGGACTCTGCCTGCGAACcgc                                                                                                                                                                                                                  >  1:496499/1‑251 (MQ=255)
                                       gggCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCC‑AATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAAt                                                                                                                                                                                                             <  2:237385/251‑1 (MQ=255)
                                        ggCAAAATCATCGGCATCTGGCCCCTGCACAAAGGCGCGATGCCGGTACACTGGTGCGGTGACTGCCGTTCTGCGCTGGCGGAAGCTGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGGCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACTTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAATCGCGCAAt                                                                                                                                                                                                             <  2:400878/251‑1 (MQ=255)
                                                                cTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACt                                                                                                                                                                                      >  1:482778/1‑250 (MQ=255)
                                                                  tcacaAAGGCGCGAAGCCGGTACACTGGTGCGTTGTCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTATCCGTCCATCGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTATCGGCCCAAGCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCGGCCTGCGAACCGCGCAATCTCGTTTGCACTGGATTGCGACTAt                                                                                                                                                                                    <  2:503381/249‑1 (MQ=255)
                                                                      aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                               >  2:115229/1‑251 (MQ=255)
                                                                      aaaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCt                                                                                                                                                                               >  1:379488/1‑251 (MQ=255)
                                                                       aaGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCCCCTGATTTCGCCTATGCGCTg                                                                                                                                                                              >  2:125876/1‑251 (MQ=255)
                                                                                   cggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACAGAAAGCAAAATTTGCCGTCAGAAAAGTTAACGGCCCAATCTCGCTGGTGATCTGTACCACTACGCCGTGGCATCTGCCTGCAAACCGCGCAATCTCTATTCCACCTGATTTCGACAATGCTCTGGAGCAGTTCGAc                                                                                                                                                                  >  2:274235/6‑251 (MQ=255)
                                                                                    ggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                                                                                                 <  1:265486/247‑1 (MQ=255)
                                                                                    ggtacACTGGTGCGTTGACTGCCGTTCGGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACg                                                                                                                                                                 <  2:149249/247‑1 (MQ=255)
                                                                                         aCTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGATGGTGCAGATCGACGGTc                                                                                                                                                              >  1:368787/1‑249 (MQ=255)
                                                                                                     aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                                                                                                <  2:86781/251‑1 (MQ=255)
                                                                                                     aCTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGAGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTg                                                                                                                                                >  2:306621/1‑251 (MQ=255)
                                                                                                             tCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAg                                                                                                                                         <  1:251561/250‑1 (MQ=255)
                                                                                                                gcgcGGGCTGAAGCTGAAGTTTAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTTAAAGCAAAATTTGCCTTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTTTATTGCATCTGATTTCGACTATGCGCTGGTGCAGATCGTCGTTCAGGTCGTGATTCTGGCGAAAGATCt                                                                                                                                     <  1:360000/251‑1 (MQ=255)
                                                                                                                          aagcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAACATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGc                                                                                                                           >  1:230879/1‑251 (MQ=255)
                                                                                                                               gaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGATTCTGGCAGAAGATCTGTTTGaagacgcaat                                                                                                                      >  2:450379/1‑243 (MQ=255)
                                                                                                                                 agTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGc                                                                                                                    >  1:117650/1‑251 (MQ=255)
                                                                                                                                      attactACGACAAATCTTTTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGTATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTTGTGATCTGGACCACTACGCCGTGGACTCTGCCTTCGAACCGCGCAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                                                                                               <  1:54215/246‑1 (MQ=255)
                                                                                                                                      aggaTTACGACAAAACTTCTCCGTCCAACGTCGTTGCTGTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCCCTGGTGAGCTGGACCACTACGCCGTGGACTCTGCCGGCGAACCGCGCAATCTCTATTGCACCTGATTTCGAATATGCGCTGGTGCAGATCGACGGTCAGGTCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTACTGCAGCGt                                                                                                               <  2:351610/248‑1 (MQ=255)
                                                                                                                                      aGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGt                                                                                                               <  2:425430/251‑1 (MQ=255)
                                                                                                                                                  aaaCTTCTCCGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTggcc                                                                                                   >  1:228682/1‑248 (MQ=255)
                                                                                                                                                          ccggccATCGACGTTTCTTTT‑CAGGCGGACGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGACCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCGGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGGATCTGCGTGACCGATTac                                                                                             <  2:225698/245‑1 (MQ=255)
                                                                                                                                                           cGTCCATCGACGTTGCTTTT‑CAGGCGGTCGGTCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctg                                                                                                                                                                                                                              >  1:277680/1‑119 (MQ=255)
                                                                                                                                                           cGTCCATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACAcc                                                                                          <  1:78738/251‑1 (MQ=255)
                                                                                                                                                           cGTCCAGCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctg                                                                                                                                                                                                                              <  2:277680/119‑1 (MQ=255)
                                                                                                                                                            gtccctctaCGTTGCTTTT‑CAGGCGGTCGACCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTCTTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCTAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCa                                                                                         <  2:281383/243‑1 (MQ=255)
                                                                                                                                                              ccattggcGTTGCTTTT‑CAGGCGGTCGATTAGGATGCACAGAAAGCAAAATTTTCCGTAAGCAACGTTAACGGCCCAATCTATCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCTCCGGAGTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCAt                                                                                        <  2:278228/243‑1 (MQ=255)
                                                                                                                                                               cATCGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTTACCGATTACAcccttt                                                                                      >  2:6420/1‑247 (MQ=255)
                                                                                                                                                                agcGACGGTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                                                                                     <  1:267733/249‑1 (MQ=255)
                                                                                                                                                                  cGACGTTGCTTTT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtct                                                                                     <  1:356569/249‑1 (MQ=255)
                                                                                                                                                                  cGACGTTGCGTTT‑CAGTCGGTCGATCAGGATGCACTGAAAGCAAAAGTTTCCGTAAGAAACGTTAACGGCCCAATCTCGCGGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAACTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATtctca                                                                                   <  2:460447/251‑2 (MQ=255)
                                                                                                                                                                    aCGTTGCTTGT‑CAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCGGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagc                                                                                 <  1:74963/251‑4 (MQ=255)
                                                                                                                                                                         gCTTTT‑CAGGCGGTCGATCAGGATGCCCTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACgg                                                                             <  2:32924/250‑1 (MQ=255)
                                                                                                                                                                              ttcAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTaaa                                                                         >  1:445739/3‑250 (MQ=255)
                                                                                                                                                                                cAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAg                                                                       <  1:243926/250‑1 (MQ=255)
                                                                                                                                                                                   gcggTCGATCAGGAGGCACTAAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACGACGTCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTGCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGTCGAAATATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGc                                                                   <  2:45915/248‑1 (MQ=255)
                                                                                                                                                                                         gATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCt                                                             >  1:338687/1‑251 (MQ=255)
                                                                                                                                                                                            cAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTg                                                           <  1:362273/250‑1 (MQ=255)
                                                                                                                                                                                              ggATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCCTTCTCAGCACGGGAAAAGGTGCGGAGCGTGAgt                                                        >  1:265810/1‑250 (MQ=255)
                                                                                                                                                                                                  gCACTGAAAGCAAAATTTTCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCGGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGTCGGGCAGGCCGTGATTCTGTCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTAGCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTTCGGAGCTTGAgttgct                                                    <  2:399939/251‑1 (MQ=255)
                                                                                                                                                                                                      tGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCAGTGGACTATGCCTGCGAACCGCGCACTCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTTGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATAGGCGTGCCCGACTACACCATTCTCAGCACGGCAAAAggggaggggcttggactgtcgcgc                                                >  2:215072/1‑229 (MQ=255)
                                                                                                                                                                                                        aaaGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGCTTACACCATTCTCAGCACGGTAAAAGGTg                                                                    >  2:423572/1‑229 (MQ=255)
                                                                                                                                                                                                        aaaGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTg                                                                    <  1:423572/229‑1 (MQ=255)
                                                                                                                                                                                                        aaaGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACGACGCCGTCGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAATGTGCGGAGCTTGAGTTGCTGCGCtt                                              <  2:123954/251‑1 (MQ=255)
                                                                                                                                                                                                           gCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAGTCTCTATTGCACCAGATTTCGACTATGCGTTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑Ac                                          >  1:47707/1‑251 (MQ=255)
                                                                                                                                                                                                            cAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTCCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGAGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGAACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑Acc                                         >  1:225499/1‑251 (MQ=255)
                                                                                                                                                                                                            cAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGTCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑Acc                                         <  1:255105/251‑1 (MQ=255)
                                                                                                                                                                                                             aaaaTTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGGTTACCCAATTTTCAGCACGGTAAAAGGGGCGGAGCTTGAGTTGCTGCGCTTT‑Accc                                        >  2:101602/1‑251 (MQ=255)
                                                                                                                                                                                                               aaTTTGCCGTAAGCAACGTTAACGGCCCAATCGCGCTGGTTATCTGTACCACTACGCTGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGTATGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACTGGAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCAt                                      <  1:53061/251‑1 (MQ=255)
                                                                                                                                                                                                                      cGTAAGCAACGTTAACGGCCCAATCTCGCTGGCGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTAAACCATTCTCAGCACGGTAAAAGGTGCGGGGCTTGAGTTGCTGCGCTTT‑ACCCATCCGttt                                >  2:25972/1‑250 (MQ=255)
                                                                                                                                                                                                                               cGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGGGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCg                      >  1:423073/1‑251 (MQ=255)
                                                                                                                                                                                                                                  tAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACAATTCTCAGCACGGTAAAAGGTGCGGAGCtt                                                            >  1:346162/1‑211 (MQ=255)
                                                                                                                                                                                                                                     cGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTc                <  1:278557/251‑1 (MQ=255)
                                                                                                                                                                                                                                     cGGCCCAATCTCGCTGGGGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTc                <  2:220551/251‑1 (MQ=255)
                                                                                                                                                                                                                                       gCCCAATCTCGCTGGTGAGCTGGACCCCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCg              <  2:397928/251‑1 (MQ=255)
                                                                                                                                                                                                                                                cGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGTCGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc     <  1:11877/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCTGGTGATCTGGACCATGACGCCGTGTACGCTGCCTGCGAACCGCGCAATCTCTAGTGCTCCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTTTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCg    <  1:328995/251‑1 (MQ=255)
                                                                                                                                                                                                                                                   tGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGAATACACCATTCTCAGCACGGTACAAGGTGCGGAGCTTGAGTTGATGCGCTTTCACCC‑TCCGTTTATGGGCTTCGGCGATCCGGCAATCCTCGGc  >  2:395486/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
CATGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTC‑CAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTT‑ACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGC  >  NC_000913/22858‑23350

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: ‑