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,1040AG98.5% 218.2 / NA 67V238V (GTA→GTGileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (40/26);  total (41/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.98e-01

TGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCG  >  NC_000913/22860‑23348
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
tGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTCCGACAAAACTTCTCCGTCCATCGACGCTGCTTTTCAGGCGGTCGATCAGGATGCCCTGAAAGCACAATTTGCCGTAATCACCGTTAACGGCCAAATCTCGCTGGTAAttagg                                                                                                                                                                                                                                                >  2:154179/1‑247 (MQ=255)
   aCTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGCAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTTGacc                                                                                                                                                                                                                                             >  1:199055/1‑251 (MQ=255)
               aaGCCAACAACATCCGCGCGCTGGGCAAAATCAGCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                                                                                                 <  2:484539/251‑1 (MQ=255)
                aGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTgg                                                                                                                                                                                                                                 >  1:5974/1‑250 (MQ=255)
                 gCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCAGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAAAGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGAc                                                                                                                                                                                                                               >  1:12686/1‑251 (MQ=255)
                    atcaccatcCGCGCGCTGGGCAAAATCCTCGGAAACGGTAACCTGCAAAAAGGCGAGAAGCCGGTAAATGGGTGCGTTGACTGCCGTGCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGActct                                                                                                                                                                                                                            <  2:422941/246‑1 (MQ=255)
                    aaCATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGCCAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGActct                                                                                                                                                                                                                            >  1:298042/1‑251 (MQ=255)
                       atcatcCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTcagcc                                                                                                                                                                                                                         >  1:242823/1‑251 (MQ=255)
                               cgcgcTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcctgcgaccc                                                                                                                                                                                                                 >  1:103800/1‑246 (MQ=255)
                                    tGGGCAAAATCATCGGCAACGGTCTCCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCgcgc                                                                                                                                                                                                             >  1:123181/1‑250 (MQ=255)
                                           aaTCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTa                                                                                                                                                                                                     <  1:239649/251‑1 (MQ=255)
                                            atcatcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTAt                                                                                                                                                                                                    <  1:227793/251‑1 (MQ=255)
                                               atcGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGc                                                                                                                                                                                                 >  1:143641/1‑251 (MQ=255)
                                                 cGGCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAc                                                                                                                                                                                               >  1:296228/1‑251 (MQ=255)
                                                   gCAACGGTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCAcct                                                                                                                                                                                             <  1:394655/251‑2 (MQ=255)
                                                        ggTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGAttt                                                                                                                                                                                        >  1:117997/1‑251 (MQ=255)
                                                        ggTCACCTGCACAAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGAttt                                                                                                                                                                                        >  2:98103/1‑251 (MQ=255)
                                                                   caAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcc                                                                                                                                                                                                                         >  1:12921/1‑207 (MQ=255)
                                                                   caAAGGCGCGAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTctgcc                                                                                                                                                                                                                         <  2:12921/207‑1 (MQ=255)
                                                                        gcgcGAAGACGGTACACTGGTGCGTTGACTGCCGTCCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAACACAAAATTTGCCGTAAGCAACGTTAAAGGCCCAGTCTCGCTGGTGATCTGTCCCCCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCAATTTCCCCTGCTTTCGACTATGCGCTGGTg                                                                                                                                                                        >  2:32458/1‑251 (MQ=255)
                                                                           cgAAGCCGGTACACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGCCTATGCGCTGGTGc                                                                                                                                                                       >  1:89218/1‑249 (MQ=255)
                                                                               gccggtacACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGc                                                                                                                                                                       >  1:233736/8‑245 (MQ=255)
                                                                                  ggtacACTGGTGCTTTGACTGCCGTTCTGCGCTTGCGGAAGCGGGCGTTGCGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGCTCAGGCTTCACTGAAAGCAAAATTTACCGTAAGCAACGTTAACGGCCCCCTCTCGCTGTTGATCTATACCCCTACTCCGTTGACTCTGCCTCCCAACCGCTCACTCTCTTTTTCCCCTTTTTTCGTCTATGCGCTGGTGCAggtccacc                                                                                                                                                              >  2:35812/5‑244 (MQ=255)
                                                                                                     gtcggTTCTGCGCTTGCGGAAGCGGAAGTTGAGTATTACGACAAAACCTCTCCGTCCATCGACGTGGAGTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTCGCCGTAAGCAACGTGAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAAGCTCTACTGCACCTGATTTCGTCTAGGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTgga                                                                                                                                           <  2:166884/247‑2 (MQ=255)
                                                                                                                   ggCGTAAGCGGAAGTTGTGTATTACGACAAAACTTCTCCGTCCAGCGACGTTGCTTTTCAGGCGGTCGATCATGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCTTGTACTCTGCCTGCGAACCGCGACATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCATATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGtt                                                                                                                              <  2:407693/250‑1 (MQ=255)
                                                                                                                      ggaagcgtaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCCCGCTGGTGATCTGGACCACTACGCCGTGGAGTCTGACTGCGAAACGGCCAATCTCTATTGCAACTGATTTCGACTACGCGCTGGTGCAGGTCGACGGTCAGGCCGTGATTTTGGCGAAAgctagagttgac                                                                                                                           >  2:153141/1‑239 (MQ=255)
                                                                                                                          gcggaagTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGt                                                                                                                      >  1:269122/1‑251 (MQ=255)
                                                                                                                                gTTGTGTATTACGACAAAACTTCCCCGTCCAACGACGTTGCTTTTCAGGCGGTAGATCAGGATGCACTGAAAGCAAAATTCGCCGTAAGCAACGTTAACGGACCAATCTCGCAGGGGATCTGGACCACTACGCCGTGGACTCTGCCTGCGACCCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTTGTCGAAAGCGTAATGCa                                                                                                                <  2:106332/251‑1 (MQ=255)
                                                                                                                                 ttGAGTATTACGCCAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGCCGGTCCGGCCGTGATTCTTGCGAAAAATCTGTTTGAAAGCGTAAtccgg                                                                                                               >  2:355536/1‑247 (MQ=255)
                                                                                                                                  tGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGCCGGTCAGGCCGTGTCTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGc                                                                                                              >  2:206350/1‑251 (MQ=255)
                                                                                                                                          ccGACAAAACTTCGCCGTCCAGCGTCGTTGCTTTGCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTTACGGCCCAATCGCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGc                                                                                                      <  2:221645/250‑1 (MQ=255)
                                                                                                                                          aCGACAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCACCGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCgg                                                                                                       <  2:273930/250‑1 (MQ=255)
                                                                                                                                              cAAAACTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTgg                                                                                                  >  2:151927/1‑250 (MQ=255)
                                                                                                                                               aaaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc                                                                                                 >  1:189799/1‑251 (MQ=255)
                                                                                                                                               aaaaCTGCACCGTCCAGCGACGTTGCTTTTCAGGCGGTCTATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAc                                                                                                 <  2:23536/251‑1 (MQ=255)
                                                                                                                                                aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGGTCTGGTTGAAAGCGTAATGCAGCGTGTCGGCGTGAcc                                                                                                >  2:127087/1‑251 (MQ=255)
                                                                                                                                                aaaCTTCTCCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGAcc                                                                                                <  1:6549/251‑1 (MQ=255)
                                                                                                                                                   cTTCTCCGTACATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCACAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACAGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGCTCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCggt                                                                                             >  2:436398/1‑249 (MQ=255)
                                                                                                                                                       tcCGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGCTTac                                                                                          >  1:260096/1‑250 (MQ=255)
                                                                                                                                                         cGTCCATCGACGTTGCTTTTCAGGCGGTCGATCAGGCTGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGAAGCGTATCGGCGTGACCGATTACAcc                                                                                       >  2:150008/1‑251 (MQ=255)
                                                                                                                                                                 gACGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcag                                                                               <  2:346915/251‑3 (MQ=255)
                                                                                                                                                                   cGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcagcc                                                                             >  1:5520/1‑247 (MQ=255)
                                                                                                                                                                   cGTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCa                                                                             <  1:366276/251‑1 (MQ=255)
                                                                                                                                                                    gTTGCTTTTCAGGCGGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctcaaccc                                                                            >  2:266037/1‑246 (MQ=255)
                                                                                                                                                                               ggcggTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTg                                                                                                                                                   <  1:475583/165‑1 (MQ=255)
                                                                                                                                                                                         cAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTg                                                        <  2:123181/250‑1 (MQ=255)
                                                                                                                                                                                         cAGGAAGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTg                                                        <  2:269122/250‑1 (MQ=255)
                                                                                                                                                                                            gATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctca                                                                                >  1:76022/1‑222 (MQ=255)
                                                                                                                                                                                            gATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATtctca                                                                                <  2:76022/223‑2 (MQ=255)
                                                                                                                                                                                            gATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTAGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGCAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGGGCTTGAggt                                                    >  1:160067/1‑249 (MQ=255)
                                                                                                                                                                                               gCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGGAAAAGGTGCGGAGCTTGAgttgct                                                 >  1:64731/1‑251 (MQ=255)
                                                                                                                                                                                                cACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGGAAAAGGTGCGGAGCTTGAGttgctg                                                >  1:411951/1‑251 (MQ=255)
                                                                                                                                                                                                      aaGCAAAATTTGCAGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGAGCTGGACCACTACGCCGTGGTCTCTGTCTGCGAACCGCGCAATCTCTAGTGCACCTGATTTCGATTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGAACTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCttt                                          <  2:298042/251‑1 (MQ=255)
                                                                                                                                                                                                        gCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCACGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTAc                                        >  2:461066/1‑251 (MQ=255)
                                                                                                                                                                                                           aaaTTTGCAGGAATCAACGGTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAAGCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGAGCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCGCAGCACGGTAAAAGGTGCGGAGCTTGATTTGCTGCGCTTTACCCa                                     <  2:209806/251‑1 (MQ=255)
                                                                                                                                                                                                                tGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCg                                 <  1:197707/250‑1 (MQ=255)
                                                                                                                                                                                                                tGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGATGCGGAGCTTGAATTGCTGCGCTTTACCCATCCGt                                >  1:234697/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACCCCATTCTCAGCACGGTGAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTAt                            >  2:116874/1‑251 (MQ=255)
                                                                                                                                                                                                                    gTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGAATCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGCTTTCGACTATGCGCTGGTGCAGATCGACGGTCAAGCCGTGATTCTGGCGAAAGATCTGGTTGAACGCGTAATGCAGCGTATCGGCGTGACCGATTCCACCATTCTCAGCGCGGTCAAAGGTGCGGAGATTTAGTTGCTGCGCTTTACCCCTCCGttttt                            >  2:288516/1‑249 (MQ=255)
                                                                                                                                                                                                                     tAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATg                           >  1:148243/1‑251 (MQ=255)
                                                                                                                                                                                                                       aGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATACACCATTCTCAGCACGGTAAAAGGGGCGGGGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGc                        >  1:227496/1‑251 (MQ=255)
                                                                                                                                                                                                                          aaCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGGGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCATATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTAGCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGTAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCtt                      <  2:401874/251‑1 (MQ=255)
                                                                                                                                                                                                                            cGTTAACGTCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCg                    <  1:291916/251‑1 (MQ=255)
                                                                                                                                                                                                                             gTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGGAATGCAGCGTATCGGCGTGGCCGATTACACCCATCTCAGCACGGGAAAAGGGGCGGAGCGTGGATTGCTGCGCTTTACCCATCCGTTTATGGGCTTCgc                   >  1:238813/1‑250 (MQ=255)
                                                                                                                                                                                                                                    gCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCg            >  2:8719/1‑251 (MQ=255)
                                                                                                                                                                                                                                             cGCTGGTGATCTGGACCCCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCg  <  2:227496/251‑1 (MQ=255)
                                                                                                                                                                                                                                             cGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTc   <  1:254483/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                    
TGGACTTCAAAACTGAAGCCAACATCATCCGCGCGCTGGGCAAAATCATCGGCAACGGTCACCTGCACAAAGGCGCGAAGCCAGTTCACTGGTGCGTTGACTGCCGTTCTGCGCTGGCGGAAGCGGAAGTTGAGTATTACGACAAAACTTCTCCGTCCATCGACGTTGCTTTCCAGGCAGTCGATCAGGATGCACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCG  >  NC_000913/22860‑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: