Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,116 T→C Y594H (TAT→CAT)  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,1160TC98.4% 187.3 / NA 64Y594H (TAT→CAT) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (29/34);  minor base A (1/0);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.76e-01

AATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCT  >  NC_000913/1870‑2361
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
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)
                                                                                                                                                                                                                                             tctgTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATc       >  1:191307/4‑251 (MQ=255)
                                                                                                                                                                                                                                                  tGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCt  >  1:521266/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
AATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACACCAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCT  >  NC_000913/1870‑2361

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: