Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,108 A→T Q591L (CAG→CTG)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,1080AT98.5% 209.1 / ‑6.0 66Q591L (CAG→CTG) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base T (30/35);  minor base G (0/1);  total (30/36)
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.45e-01

TGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGA  >  NC_000913/1862‑2354
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
tGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                                                                                                     >  1:309423/1‑251 (MQ=255)
tGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTAATCCCATCTGTTg                                                                                                                                                                                                                                                     >  1:327438/1‑251 (MQ=255)
 gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                                                                                                     <  1:191830/250‑1 (MQ=255)
 gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACGTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                                                                                                                                                                      <  1:370354/249‑1 (MQ=255)
        aaTGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAATCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCTCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTACCATCTGTTGCGTCATgc                                                                                                                                                                                                                                             <  2:485518/251‑1 (MQ=255)
        aaTGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTGCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATgc                                                                                                                                                                                                                                             <  1:306942/251‑1 (MQ=255)
                   cctaaaTCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAt                                                                                                                                                                                                                                  <  1:488778/247‑1 (MQ=255)
                         ggtGGAAAACTGGCAGGAAGACCTGGCCCAGGCCAAAGAGCCGTTTAATGTCGGGCGCTTAATTCGCGTCGTGAAAGAATATCATCTGCTGGACCCGGTCATTGTTGACTGCACTTCCAGACAGGCAGTGGCGGAGCAATATGCCGACTTCGTGCGCGTAGGTTTCCACGTTGGCGCTCTGAGCAAAAAGGCCAACACCTCGTCGATGGATTACTACCAGCTGTTGCGTCATGCGGCGGAATAATCgcggt                                                                                                                                                                                                                            <  2:206870/249‑2 (MQ=255)
                                  cTGGCAGGAAGAACTGGCCCAAGCCAAAGAGCCGTTTAATATCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCCCGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGAGGGCGAAAAATCGGGGCGTAAATTcc                                                                                                                                                                                                                   >  2:398032/1‑251 (MQ=255)
                                     gCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCctct                                                                                                                                                                                                                >  1:385052/1‑251 (MQ=255)
                                       aGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAt                                                                                                                                                                                                              >  2:107739/1‑251 (MQ=255)
                                                tGGCGCACGCCACAGAGCCGTTTCATCTCGGGCGCTTCATTCGCCTCGTGAAAGAATATCATCTGCTGCACCCGGTCATTGTTGACTGCCCTTCCAGCCAGGCAGTGGCGGATCAATCTGCCGACTTCTTGCCCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCCCGTCGCTGGATTACTACCATCTGTTGCTGAATTGGGCGGAAAAATCACGGCGTACATTACTCTATGCCACCaa                                                                                                                                                                                                     >  2:607404/1‑250 (MQ=255)
                                                 ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACg                                                                                                                                                                                                   <  2:128922/251‑1 (MQ=255)
                                                 ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCGATGACACCAAc                                                                                                                                                                                                    <  2:79220/250‑1 (MQ=255)
                                                    gcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTg                                                                                                                                                                                                <  1:274815/251‑1 (MQ=255)
                                                     cAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTgg                                                                                                                                                                                               <  1:76037/251‑1 (MQ=255)
                                                            aaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGa                                                                                                                                                                                        <  1:454678/251‑1 (MQ=255)
                                                                    tttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTTGATTACCGGt                                                                                                                                                                                >  2:374643/1‑251 (MQ=255)
                                                                                     aaTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaa                                                                                                                                                                >  1:471191/1‑250 (MQ=255)
                                                                                           ccTCGTGGAAGAATATCATCTGCTGAACCCGGTCATTGTTGCCTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCCACGTTGGGGCTGGATTACCGGTTTTTGGGAACCTGCAAAATCTGc                                                                                                                                                         >  2:360149/1‑251 (MQ=255)
                                                                                              cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGGCGAATTCTTGCGCGAAGGTTTCCAAGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTAACCTCTGTTGCGTCCTGCGGCGGGAAAATCGCGGCGTAAATTCCTCTCTGACACCAACGTTGGGGCTGGGTTACCGGTTCTTGAAAACCTGCAAAATCTggtcc                                                                                                                                                      >  1:350956/1‑247 (MQ=255)
                                                                                                      aaTATCATCTGCTGAACCCGGTCATTGTTGCCTGCACTTCCAGCCCAGCAGTGGCGGATCTATATGTCGACTTCTAGCGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCCGAAAAATCGCGGTGTAAATTCATTCATGACACCAAAGATGGGGCTGGATTAACGGGTATTGAGAACCTGCAAAATATGCTCAATGctggg                                                                                                                                              >  2:105358/1‑247 (MQ=255)
                                                                                                       atatCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtg                                                                                                                                             <  1:444457/251‑1 (MQ=255)
                                                                                                        tatCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGGAAATTCCTCTATGCCACCAACGTTGGGGCTGGATTAACGGCTATTGAGGACCTGCAAAATATGCttcatgccggtgg                                                                                                                                            >  2:267696/1‑239 (MQ=255)
                                                                                                         atcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgat                                                                                                                                           <  1:57667/251‑1 (MQ=255)
                                                                                                           catcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgatga                                                                                                                                         <  1:77990/251‑1 (MQ=255)
                                                                                                            atcTGCTGTACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCCGTGGCGGATCAATATGACGACTTCTTGCGCGAAGGTTTCCCCGTTGTCACCCCGAACAAAAAGGGCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGGCAAGCAGCGGAAAAATTGCGGCGTAAAGTCATCTATGAAACACAAGATGGGGCTGGATTACCGGTTATTGAGAAACTGCAGCATCTGCTCAAtacgggggatgg                                                                                                                                         >  2:545917/1‑239 (MQ=255)
                                                                                                                gctgAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAGGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATtg                                                                                                                                     <  1:314137/250‑1 (MQ=255)
                                                                                                                   gAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGACCCTGCAAAATCTGCTCAATGCTGGTGATGAATtgatg                                                                                                                                  <  1:100849/250‑1 (MQ=255)
                                                                                                                     aCCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAg                                                                                                                               <  1:473525/251‑1 (MQ=255)
                                                                                                                      cccGGTCATTGTTGACTGCACTTCCAGCGAGGCATGTGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAATGCCAACACCTCGTCGATGGATTACTACCGGCGGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGGTGGGGCTGGATTACCGGGGAGGGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGt                                                                                                                              <  1:228698/251‑1 (MQ=255)
                                                                                                                       acggccATTGTTGACTGAACTTCCAGCCAGGCAGTGGCGGACCCATATGGCGACTTCTTGCGCGAAGGTTTCCACCTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCATCATGCGGCGGAAAACTCGCGGGGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGACAAATCTGCTCAATGCTGGTGATGACTTGCAGAAGtt                                                                                                                             >  2:528503/6‑251 (MQ=255)
                                                                                                                          gTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACGTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCTGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAGCCTGCAAAATCGGCTCAATGCTGGTTATGAATTGATGGAGTtctc                                                                                                                          <  2:36198/251‑1 (MQ=255)
                                                                                                                            cATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTcc                                                                                                                         <  1:572217/250‑1 (MQ=255)
                                                                                                                             aTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGTATTGATGAAGTTCTCCgg                                                                                                                       >  2:521732/1‑251 (MQ=255)
                                                                                                                               ggtttaCTGCACTTCCATCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGTATTACCGGTGATTGAGAACCTGCAAAATCTGCTCAATGCTGGTTATTAATGGATGAAGTtcggcggt                                                                                                                      <  2:30299/245‑7 (MQ=255)
                                                                                                                                 ttgACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAtt                                                                                                                   >  2:160704/1‑251 (MQ=255)
                                                                                                                                  tgACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGGGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAAAGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGAAttc                                                                                                                  >  1:535943/1‑251 (MQ=255)
                                                                                                                                    aCTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCCCGGTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCCACGTTGGGGCTGGATTACCGGTTCTTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                                >  2:129073/1‑251 (MQ=255)
                                                                                                                                    aCTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                                >  1:197974/1‑251 (MQ=255)
                                                                                                                                      tGCACTTCCAGCCAGGCAGTGGCGGATTAATATGCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGGTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGGTTACCGGTTATTGGGAACCTGCAAAATCTGCTCAATGATGGTGATGAATTGATGAAGTTCTCCGGCAttctttc                                                                                                              >  2:390980/1‑251 (MQ=255)
                                                                                                                                       gCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                                                                                                             <  1:284524/251‑2 (MQ=255)
                                                                                                                                          cTTCCAGCCAGGCAGTGGCGGATCAATATGCCCACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTACTCTATGACACCAAAGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAACATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGt                                                                                                          >  1:118569/1‑251 (MQ=255)
                                                                                                                                                 ccagGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCTTCGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCTTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGTGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                                                                   <  1:278152/251‑1 (MQ=255)
                                                                                                                                                      cAGTGGCGGATCAATATGCCGACTTCTTGCGCGGAGGGTTCCACGGTGTCCCGCCGAACAAAAAAGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGa                                                                                                                                                                                                            >  2:61587/1‑142 (MQ=255)
                                                                                                                                                      cAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGa                                                                                                                                                                                                            <  1:61587/142‑1 (MQ=255)
                                                                                                                                                          ggCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGACGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCt                                                                                          >  1:89774/1‑251 (MQ=255)
                                                                                                                                                          ggCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCt                                                                                          <  1:374643/251‑1 (MQ=255)
                                                                                                                                                                      atGCCGACTTCTTGCGCGAATGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGTGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTTCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGtt                                                                                <  1:251784/250‑1 (MQ=255)
                                                                                                                                                                                     gcgAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGt                                                                <  1:521732/251‑1 (MQ=255)
                                                                                                                                                                                         aGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtct                                                            <  2:201273/251‑1 (MQ=255)
                                                                                                                                                                                             ttCCACGTTGTCACGCCGAACAAAAAGGCCAATACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCATGTTAGACGAAGGCATGAGTTTCTCCgt                                                        >  2:112968/1‑250 (MQ=255)
                                                                                                                                                                                                   gTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGAc                                                  >  2:132251/1‑251 (MQ=255)
                                                                                                                                                                                                     tGTCACGCCGAACAAAACGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGAct                                                 >  2:83810/1‑249 (MQ=255)
                                                                                                                                                                                                        cACGCCGAACAAAAAGGCCAACCCCTCGTCGTTGGATTACTACCATCTGTTGCGTCATGCGGCGGACAAATCGCGGCGTAAATTCCTATATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCACAATCTGCTCAATGCTGGTCCTGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGActactc                                             >  2:395204/1‑246 (MQ=255)
                                                                                                                                                                                                                 cAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCggg                                    <  1:79413/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGaa                                  <  2:306512/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGaaa                                 <  1:553686/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCCGCGGAAACATCGCGGCGTAAATTCCTCCATGACACCAAAGTTGCAGCTGGATTACCGGCTATTGAGAAACTGAAAAATCTGCTCAATGCTGGTGCTGACTTGCTGAAGTTCTCCGGCATTCTTTCAGGTTAGCTTTCTTATATCTTCGGCAAGTTCGACGAAGCCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATggg                             >  1:321319/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cgATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAAcc                 <  1:66435/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gATGGATTACTACCATCTGTTGCGTCAAGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACAAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCgg               >  1:461446/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGa              <  1:390980/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tactacCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCcgc         <  1:107739/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 acCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgat    <  2:471191/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgatg   >  1:393869/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgatga  <  1:40464/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
TGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGA  >  NC_000913/1862‑2354

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