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

TGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGA  >  NC_000913/1862‑2354
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
tGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                                                                                                     >  1:309423/1‑251 (MQ=255)
tGCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTAATCCCATCTGTTg                                                                                                                                                                                                                                                     >  1:327438/1‑251 (MQ=255)
 gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTg                                                                                                                                                                                                                                                     <  1:191830/250‑1 (MQ=255)
 gCTCACCAATGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACGTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGtt                                                                                                                                                                                                                                                      <  1:370354/249‑1 (MQ=255)
        aaTGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAATCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCTCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTACCATCTGTTGCGTCATgc                                                                                                                                                                                                                                             <  2:485518/251‑1 (MQ=255)
        aaTGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTGCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATgc                                                                                                                                                                                                                                             <  1:306942/251‑1 (MQ=255)
                   cctaaaTCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAAt                                                                                                                                                                                                                                  <  1:488778/247‑1 (MQ=255)
                         ggtGGAAAACTGGCAGGAAGACCTGGCCCAGGCCAAAGAGCCGTTTAATGTCGGGCGCTTAATTCGCGTCGTGAAAGAATATCATCTGCTGGACCCGGTCATTGTTGACTGCACTTCCAGACAGGCAGTG‑GCGGAGCAATATGCCGACTTCGTGCGCGTAGGTTTCCACGTTGGCGCTCTGAGCAAAAAGGCCAACACCTCGTCGATGGATTACTACCAGCTGTTGCGTCATGCGGCGGAATAATCgcggt                                                                                                                                                                                                                            <  2:206870/249‑2 (MQ=255)
                                  cTGGCAGGAAGAACTGGCCCAAGCCAAAGAGCCGTTTAATATCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCCCGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGAGGGCGAAAAATCGGGGCGTAAATTcc                                                                                                                                                                                                                   >  2:398032/1‑251 (MQ=255)
                                     gCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCctct                                                                                                                                                                                                                >  1:385052/1‑251 (MQ=255)
                                       aGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCCTCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAt                                                                                                                                                                                                              >  2:107739/1‑251 (MQ=255)
                                                tGGCGCACGCCACAGAGCCGTTTCATCTCGGGCGCTTCATTCGCCTCGTGAAAGAATATCATCTGCTGCACCCGGTCATTGTTGACTGCCCTTCCAGCCAGGCAGTG‑GCGGATCAATCTGCCGACTTCTTGCCCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCCCGTCGCTGGATTACTACCATCTGTTGCTGAATTGGGCGGAAAAATCACGGCGTACATTACTCTATGCCAC‑Caa                                                                                                                                                                                                     >  2:607404/1‑250 (MQ=255)
                                                 ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACg                                                                                                                                                                                                   <  2:128922/251‑1 (MQ=255)
                                                 ggCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCGATGACAC‑CAAc                                                                                                                                                                                                    <  2:79220/250‑1 (MQ=255)
                                                    gcAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTg                                                                                                                                                                                                <  1:274815/251‑1 (MQ=255)
                                                     cAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTgg                                                                                                                                                                                               <  1:76037/251‑1 (MQ=255)
                                                            aaGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGa                                                                                                                                                                                        <  1:454678/251‑1 (MQ=255)
                                                                    tttAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTTGATTACCGGt                                                                                                                                                                                >  2:374643/1‑251 (MQ=255)
                                                                                     aaTTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCaa                                                                                                                                                                >  1:471191/1‑250 (MQ=255)
                                                                                           ccTCGTGGAAGAATATCATCTGCTGAACCCGGTCATTGTTGCCTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CCACGTTGGGGCTGGATTACCGGTTTTTGGGAACCTGCAAAATCTGc                                                                                                                                                         >  2:360149/1‑251 (MQ=255)
                                                                                              cGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGGCGAATTCTTGCGCGAAGGTTTCCAAGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGCTTACTAACCTCTGTTGCGTCCTGCGGCGGGAAAATCGCGGCGTAAATTCCTCTCTGACAC‑CAACGTTGGGGCTGGGTTACCGGTTCTTGAAAACCTGCAAAATCTggtcc                                                                                                                                                      >  1:350956/1‑247 (MQ=255)
                                                                                                      aaTATCATCTGCTGAACCCGGTCATTGTTGCCTGCACTTCCAGCCCAGCAGTG‑GCGGATCTATATGTCGACTTCTAGCGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCCGAAAAATCGCGGTGTAAATTCATTCATGACAC‑CAAAGATGGGGCTGGATTAACGGGTATTGAGAACCTGCAAAATATGCTCAATGctggg                                                                                                                                              >  2:105358/1‑247 (MQ=255)
                                                                                                       atatCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtg                                                                                                                                             <  1:444457/251‑1 (MQ=255)
                                                                                                        tatCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGGAAATTCCTCTATGCCAC‑CAACGTTGGGGCTGGATTAACGGCTATTGAGGACCTGCAAAATATGCttcatgccggtgg                                                                                                                                            >  2:267696/1‑239 (MQ=255)
                                                                                                         atcatcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgat                                                                                                                                           <  1:57667/251‑1 (MQ=255)
                                                                                                           catcTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGtgatga                                                                                                                                         <  1:77990/251‑1 (MQ=255)
                                                                                                            atcTGCTGTACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCCGTG‑GCGGATCAATATGACGACTTCTTGCGCGAAGGTTTCCCCGTTGTCACCCCGAACAAAAAGGGCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGGCAAGCAGCGGAAAAATTGCGGCGTAAAGTCATCTATGAAACACAA‑GATGGGGCTGGATTACCGGTTATTGAGAAACTGCAGCATCTGCTCAAtacgggggatgg                                                                                                                                         >  2:545917/1‑239 (MQ=255)
                                                                                                                gctgAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTAGGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATtg                                                                                                                                     <  1:314137/250‑1 (MQ=255)
                                                                                                                   gAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGACCCTGCAAAATCTGCTCAATGCTGGTGATGAATtgatg                                                                                                                                  <  1:100849/250‑1 (MQ=255)
                                                                                                                     aCCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAg                                                                                                                               <  1:473525/251‑1 (MQ=255)
                                                                                                                      cccGGTCATTGTTGACTGCACTTCCAGCGAGGCA‑TGTGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAATGCCAACACCTCGTCGATGGATTACTACCGGCGGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGGTGGGGCTGGATTACCGGGGAGGGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGt                                                                                                                              <  1:228698/251‑1 (MQ=255)
                                                                                                                       acggccATTGTTGACTGAACTTCCAGCCAGGCAGTG‑GCGGACCCATATGGCGACTTCTTGCGCGAAGGTTTCCACCTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCATCATGCGGCGGAAAACTCGCGGGGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGACAAATCTGCTCAATGCTGGTGATGACTTGCAGAAGtt                                                                                                                             >  2:528503/6‑251 (MQ=255)
                                                                                                                          gTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACGTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCTGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAGCCTGCAAAATCGGCTCAATGCTGGTTATGAATTGATGGAGTtctc                                                                                                                          <  2:36198/251‑1 (MQ=255)
                                                                                                                            cATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTcc                                                                                                                         <  1:572217/250‑1 (MQ=255)
                                                                                                                             aTTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGTATTGATGAAGTTCTCCgg                                                                                                                       >  2:521732/1‑251 (MQ=255)
                                                                                                                               ggtttaCTGCACTTCCATCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCGCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGTATTACCGGTGATTGAGAACCTGCAAAATCTGCTCAATGCTGGTTATTAATGGATGAAGTtcggcggt                                                                                                                      <  2:30299/245‑7 (MQ=255)
                                                                                                                                 ttgACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAtt                                                                                                                   >  2:160704/1‑251 (MQ=255)
                                                                                                                                  tgACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGGGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAAAGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGAAttc                                                                                                                  >  1:535943/1‑251 (MQ=255)
                                                                                                                                    aCTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCCCGGTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CCACGTTGGGGCTGGATTACCGGTTCTTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                                >  2:129073/1‑251 (MQ=255)
                                                                                                                                    aCTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctt                                                                                                                >  1:197974/1‑251 (MQ=255)
                                                                                                                                      tGCACTTCCAGCCAGGCAGTG‑GCGGATTAATATGCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGGTTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGGTTACCGGTTATTGGGAACCTGCAAAATCTGCTCAATGATGGTGATGAATTGATGAAGTTCTCCGGCAttctttc                                                                                                              >  2:390980/1‑251 (MQ=255)
                                                                                                                                       gCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCAttctttca                                                                                                             <  1:284524/251‑2 (MQ=255)
                                                                                                                                          cTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCCACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAAAACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTACTCTATGACAC‑CAAAGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAACATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGt                                                                                                          >  1:118569/1‑251 (MQ=255)
                                                                                                                                                 ccagGCAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCTTCGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCTTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGTGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCttt                                                                                                   <  1:278152/251‑1 (MQ=255)
                                                                                                                                                      cAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGGAGGGTTCCACGGTGTCCCGCCGAACAAAAAAGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGa                                                                                                                                                                                                            >  2:61587/1‑142 (MQ=255)
                                                                                                                                                      cAGTG‑GCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGa                                                                                                                                                                                                            <  1:61587/142‑1 (MQ=255)
                                                                                                                                                          g‑gCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGACGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCt                                                                                          >  1:89774/1‑251 (MQ=255)
                                                                                                                                                          g‑gCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCt                                                                                          <  1:374643/251‑1 (MQ=255)
                                                                                                                                                                      atGCCGACTTCTTGCGCGAATGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGTGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTTCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGtt                                                                                <  1:251784/250‑1 (MQ=255)
                                                                                                                                                                                     gcgAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGt                                                                <  1:521732/251‑1 (MQ=255)
                                                                                                                                                                                         aGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTtct                                                            <  2:201273/251‑1 (MQ=255)
                                                                                                                                                                                             ttCCACGTTGTCACGCCGAACAAAAAGGCCAATACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCATGTTAGACGAAGGCATGAGTTTCTCCgt                                                        >  2:112968/1‑250 (MQ=255)
                                                                                                                                                                                                   gTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGAc                                                  >  2:132251/1‑251 (MQ=255)
                                                                                                                                                                                                     tGTCACGCCGAACAAAACGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGAct                                                 >  2:83810/1‑249 (MQ=255)
                                                                                                                                                                                                        cACGCCGAACAAAAAGGCCAACCCCTCGTCGTTGGATTACTACCATCTGTTGCGTCATGCGGCGGACAAATCGCGGCGTAAATTCCTATATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCACAATCTGCTCAATGCTGGTCCTGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGActactc                                             >  2:395204/1‑246 (MQ=255)
                                                                                                                                                                                                                 cAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCggg                                    <  1:79413/251‑1 (MQ=255)
                                                                                                                                                                                                                   aaaaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGaa                                  <  2:306512/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaaGGCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGaaa                                 <  1:553686/251‑1 (MQ=255)
                                                                                                                                                                                                                        gCCAACACCTCGTCGATGGATTACTACCATCTGTTGCGTCATGCCGCGGAAACATCGCGGCGTAAATTCCTCCATGACAC‑CAAAGTTGCAGCTGGATTACCGGCTATTGAGAAACTGAAAAATCTGCTCAATGCTGGTGCTGACTTGCTGAAGTTCTCCGGCATTCTTTCAGGTTAGCTTTCTTATATCTTCGGCAAGTTCGACGAAGCCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATggg                             >  1:321319/1‑251 (MQ=255)
                                                                                                                                                                                                                                     cgATGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAAcc                 <  1:66435/250‑1 (MQ=255)
                                                                                                                                                                                                                                      gATGGATTACTACCATCTGTTGCGTCAAGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑AAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGGGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCgg               >  1:461446/1‑251 (MQ=255)
                                                                                                                                                                                                                                       aTGGATTACTACCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGa              <  1:390980/251‑1 (MQ=255)
                                                                                                                                                                                                                                             tactacCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCcgc         <  1:107739/250‑1 (MQ=255)
                                                                                                                                                                                                                                                 acCATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgat    <  2:471191/251‑1 (MQ=255)
                                                                                                                                                                                                                                                  ccATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgatg   >  1:393869/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   cATCTGTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAgatga  <  1:40464/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       
TGCTCACCAATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGA  >  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: ‑