Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 23,297 C→T L303L (CTG→TTG)  ileS → isoleucine‑‑tRNA ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_00091323,2970CT91.4% 84.0 / ‑1.8 39L303L (CTG→TTG) ileSisoleucine‑‑tRNA ligase
Reads supporting (aligned to +/- strand):  ref base C (2/1);  new base T (15/17);  total (19/20)
Fisher's exact test for biased strand distribution p-value = 6.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAACGTCTTCAAAGCGAACGAC  >  NC_000913/23052‑23533
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
cACTGAAAGCAAGATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGGTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAACACGGTAAAAGGTGCGGAGGTTGAGgtgctg                                                                                                                                                                                                                                           >  1:259725/1‑251 (MQ=255)
   tGAAAGCAAAATTTCCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTgcgc                                                                                                                                                                                                                                        >  1:311499/1‑251 (MQ=255)
    gAAAGCAAAATTTGCCGGAAGCAACGGTAACGGCTCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACGCTGCCTGCGAACCGCGCAATCTCTATTGCACTTGATTTCGACGATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGGAATGCAGCGTATCGGCGTGACCTATTACACCAGTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGGTGCTgcg                                                                                                                                                                                                                                         <  2:378096/249‑1 (MQ=255)
         caaaattttccgatatcaactgtAACGGCCCAATCTCGCTGGTGATCTGGACCTCTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGAGTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGGTTTAca                                                                                                                                                                                                                                 <  2:157882/229‑2 (MQ=255)
                 tccGTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGt                                                                                                                                                                                                                          <  2:168669/249‑1 (MQ=255)
                    gTAAGCAACGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTAt                                                                                                                                                                                                                      <  1:103580/251‑1 (MQ=255)
                          aaCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACAGATTAAACCATTCTCCACACGGAAAAAGGTGACGGGCTTGAGTTGCTTCGCTTTACCCATCCGTTTATGGGCtt                                                                                                                                                                                                                >  2:273366/1‑251 (MQ=255)
                           aCGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTc                                                                                                                                                                                                               <  1:368873/251‑1 (MQ=255)
                            cGTTAACGGCCCAATCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAACGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTACAAGGGGCGGAGATTGAGTTGCTGCGCTTTACCCATCCGTTTATGGACTTCg                                                                                                                                                                                                              >  1:42938/1‑251 (MQ=255)
                                        aaTCTCGCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGAGTAATGCAGCGTATCGGCGTGACCGATTACACCCTTCTCAGAACGGTAAAAGGAGCGGGGTTTGAGATGCCGCGCAGTACCCATCCGTTTATGGGCCTCGACGTACCGGcc                                                                                                                                                                                                   >  2:145304/1‑249 (MQ=255)
                                           ctcGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGATAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGGAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATcc                                                                                                                                                                                               >  1:84614/1‑251 (MQ=255)
                                              gCTGGTGATCTGGACCACTACGCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGAGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGTCGCAACGCGTTACAACATTCTCAGCACGGAAAAAGGGACGAGATTTGTGTTGGTGGGCTTTACCCATCCGTTTATGAGCtttggcggtctggaactacacg                                                                                                                                                                                            >  2:303397/1‑231 (MQ=255)
                                                                   gCCGTGGACTCTGCCTGCGAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATg                                                                                                                                                                       <  1:93284/251‑1 (MQ=255)
                                                                                   gcgaaCCGCGCAATCTCTATTGCATCTGATTTCGACTATGCGCTTGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTGCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCg                                                                                                                                                        <  1:148750/247‑1 (MQ=255)
                                                                                   gccAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGt                                                                                                                                                       >  1:189423/1‑251 (MQ=255)
                                                                                          gcgcAATCGCTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACAcc                                                                                                                                                <  1:160371/251‑1 (MQ=255)
                                                                                             cactctcTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACcgcg                                                                                                                                            <  1:235491/248‑1 (MQ=255)
                                                                                                tctcTATTGCACCTGATTTCGACTATGCGCTGGTGCAGATCGACGGTCCGGCCGTGTTTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGAGTGACAGATAACACCAATCTCCGCAAGGGAAAAGGGGTGGAGCGTGAGTTGCTCGCCTTCACCCATCCGTTTAAGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGATACCAGAGATGCCTGTACCgatcccattcccacccgcacc                                                                                                                                          >  2:310830/1‑231 (MQ=255)
                                                                                                        gCACCTGACTTCGAGTATGAGTTGGTGCAGATCGACGGTCAGGTCGTGATTCTGGCGAAAGAGCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACTCCATTCCCAGCCCGGTAAAAGGTGCGGAGCTTGAGTGGCTGCGCTTTGCCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTGCCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCAACg                                                                                                                                  <  2:283002/250‑1 (MQ=255)
                                                                                                          acctgATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCCCACCGCGCCTGGCCACGGcc                                                                                                                               >  1:128006/5‑251 (MQ=255)
                                                                                                                       tATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGGTCACACCGCGCCTGGCCACGGCCCGGACGACTAtgt                                                                                                                  >  1:313630/1‑251 (MQ=255)
                                                                                                                               ggTGAAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGCTCTGGTTGAAAGCGTAATGCAGCATATCGGCGTGACCGACTACCACATTCTCAGCACAGCAAAAGGGGCGGCGCTTGAGATTCTGAGCTTTACCCCTCCGTTTATAGGCTTCGACGTTCCGGCAGTCCGCGGCGATACCGTTACCCTGTATTCCTGTGGCGGTGCTGTTCCACCCTCGCCTGGCCAccgacaggacgctcatggtctaggac                                                                                                          >  2:116168/1‑226 (MQ=255)
                                                                                                                                            ggTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTCCCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCt                                                                                             <  1:362690/251‑1 (MQ=255)
                                                                                                                                            ggTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGccg                                                                                             >  1:306341/1‑250 (MQ=255)
                                                                                                                                             gTCAGGCCTTGGTTCTGGCGAAAGATCTGGTTGAAAGCGGAATGACGCGTATCGGCGTGACCTATTACACCAGACTCAGCACGGTAAAATGTGCGGAGCTTGAGGTGCGTCGCGTTGCCCATCCGTGTATGGGCTTCGACGGTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTAGGTGATCGGTCAGAAATACTGCCTg                                                                                            <  2:283388/251‑1 (MQ=255)
                                                                                                                                                    cGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCATGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGACTGGAAACCg                                                                                     >  1:79210/1‑251 (MQ=255)
                                                                                                                                                                         ggTTTAAAGCGTAATGCAGCGTTTCGGAGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTTCCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCg                                                                    <  2:6508/247‑1 (MQ=255)
                                                                                                                                                                                   gTAATGCTGCGTATCGGCGGGACCTATTACACCTATCTCAGCTCGGTAAATGGGGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGAGCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTATCCCGGTTGGCCCGGACGGCACTTATCt                                                      <  2:230837/251‑1 (MQ=255)
                                                                                                                                                                                          aGCGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGc                                               >  1:255810/1‑251 (MQ=255)
                                                                                                                                                                                            cGTATCGGCGTGACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCAc                                             <  1:273366/251‑1 (MQ=255)
                                                                                                                                                                                                   gCGTGACCGATTACAGCATGCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGTCAATCCTCGGCGATCACGTTACCCTGGAGGCCGGTACCGGGGCCGTTCGCACCGCGCCTGGCCACGGCCCGGACTACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCg                                      <  1:352618/251‑1 (MQ=255)
                                                                                                                                                                                                       gACCGATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGc                                  >  1:229073/1‑251 (MQ=255)
                                                                                                                                                                                                            gATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGt                                                                                                                                                                                  >  1:291884/1‑104 (MQ=255)
                                                                                                                                                                                                            gATTACACCATTCTCAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGt                                                                                                                                                                                  <  2:291884/104‑1 (MQ=255)
                                                                                                                                                                                                                       tctcAGCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAAc                    <  1:116168/250‑1 (MQ=255)
                                                                                                                                                                                                                        ctcagCACGGTAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATAGGCTTCGACGTTCCGGCAATCCTCGGCGATCCCGTTACCCTGGATGCCGGCACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGCCGACTATGTGATCGGACAGAAATAAGGCATAGAAACCGCTAAACCGGTTGGACCGGCCGGAACTTATATGACGGGCACTTATACGCCGCTGGATGGCgggcacg                   >  2:204492/5‑244 (MQ=255)
                                                                                                                                                                                                                                  tAAAAGGTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAACGTCTTCAAAGCg        >  1:92049/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGGTCCGGCAATCCTCGGCGCTCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCCCGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAACGTCTTCAAAGCGAACg    <  2:84614/251‑1 (MQ=255)
                                                                                                                                                                                                                                        gTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAACGTCTTCAAAGCGAACGAc  >  1:84239/1‑251 (MQ=255)
                                                                                                                                                                                                                                        gTGCGGAGCTTGAGTTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAACGTCTTCAAAGCGAACGAc  <  1:193650/251‑1 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
CACTGAAAGCAAAATTTGCCGTAAGCAACGTTAACGGCCCAATCTCGCTGGTAATCTGGACCACCACGCCGTGGACTCTGCCTGCCAACCGCGCAATCTCTATTGCACCAGATTTCGACTATGCGCTGGTGCAGATCGACGGTCAGGCCGTGATTCTGGCGAAAGATCTGGTTGAAAGCGTAATGCAGCGTATCGGCGTGACCGATTACACCATTCTCGGCACGGTAAAAGGTGCGGAGCTTGAGCTGCTGCGCTTTACCCATCCGTTTATGGGCTTCGACGTTCCGGCAATCCTCGGCGATCACGTTACCCTGGATGCCGGTACCGGTGCCGTTCACACCGCGCCTGGCCACGGCCCGGACGACTATGTGATCGGTCAGAAATACGGCCTGGAAACCGCTAACCCGGTTGGCCCGGACGGCACTTATCTGCCGGGCACTTATCCGACGCTGGATGGCGTGAACGTCTTCAAAGCGAACGAC  >  NC_000913/23052‑23533

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