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

AATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCT  >  NC_000913/1870‑2361
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
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)
                                                                                                                                                                                                                                             tctgTTGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATc       >  1:191307/4‑251 (MQ=255)
                                                                                                                                                                                                                                                  tGCGTCATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCTGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCAGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACTACTCTGGCGCGGGAAATGGGTTATACCGAACCGGATCCGCGAGATGATCTTTCt  >  1:521266/1‑251 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                      
AATGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTG‑GCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTACCATCAGTTGCGTTATGCGGCGGAAAAATCGCGGCGTAAATTCCTCTATGACAC‑CAACGTTGGGGCTGGATTACCGGTTATTGAGAACCTGCAAAATCTGCTCAATGCAGGTGATGAATTGATGAAGTTCTCCGGCATTCTTTCTGGTTCGCTTTCTTATATCTTCGGCAAGTTAGACGAAGGCATGAGTTTCTCCGAGGCGACCACGCTGGCGCGGGAAATGGGTTATACCGAACCGGACCCGCGAGATGATCTTTCT  >  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: ‑