Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,597 T→C V754A (GTC→GCC)  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,5970TC95.8% 123.8 / ‑1.9 48V754A (GTC→GCC) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (23/23);  total (24/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.18e-01

TTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATA  >  NC_000913/2358‑2837
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
ttCTGGTATGGATGTAGCGCGTAAGCGATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGt                                                                                                                                                                                                                                          <  1:119108/249‑1 (MQ=255)
             gtagcgcGTAAGCTATTGATTCTCGCTCGTGACACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGACGATGGGGCCTGCCGCGGGGAGGTTGCCggagg                                                                                                                                                                                                                           >  2:50490/4‑247 (MQ=255)
                 cgcgTAAGCTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGAt                                                                                                                                                                                                                       <  2:358612/251‑1 (MQ=255)
                           aTTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGCAGATTGCCGAAGTGGATGGTAATGATc                                                                                                                                                                                                             >  1:291734/1‑251 (MQ=255)
                             ttaTTCTCGGTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCTCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATcc                                                                                                                                                                                                            <  2:492039/248‑1 (MQ=255)
                               ggtcttgcgagtggaaCGGGACGCGTGCTGGGGCTGGCGGATAGTGACATTGAGCCTGTGCTGGCCGCAGAGTGTAGCGCTGGGGGGTATGTGGCCGCTTTGATGGCGAATCTGTCACAGCTCGACGATCTCTTGGCCGCGCGCGGGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTCGTCAATATTGAGGAAAATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATAGTaagtagcctct                                                                                                                                                                                                         <  2:259179/237‑10 (MQ=255)
                               aTTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTGTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCt                                                                                                                                                                                                         <  1:50490/251‑1 (MQ=255)
                                    cgctcgTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTTTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTc                                                                                                                                                                                                     <  1:181212/250‑1 (MQ=255)
                                                gggACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATg                                                                                                                                                                                       <  1:423138/251‑1 (MQ=255)
                                                  gACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGGG‑AAAAATggg                                                                                                                                                                                     >  1:119142/1‑250 (MQ=255)
                                                       gAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGaaaa                                                                                                                                                                                <  1:42596/251‑1 (MQ=255)
                                                               gCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTgg                                                                                                                                                                        <  1:100529/251‑1 (MQ=255)
                                                                   gCGGAAATTGAAATGGAACCTGTGCGGACCGCAGAGTTTTACGCTGAGGGTGATGGTGCCGCTTTTATTGCGAAGCGGTCACAGCTCGACGAACTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGGTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCGGTTCAAAGTG‑AAAAAGGGCGAAAACGCCCTGGCCtt                                                                                                                                                                    <  2:39345/251‑1 (MQ=255)
                                                                    cGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCttt                                                                                                                                                                   >  1:154563/1‑250 (MQ=255)
                                                                       atatTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGGTGCCGCTTTTAGGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACTCCCTGGccgtttat                                                                                                                                                                <  1:403650/251‑7 (MQ=255)
                                                                           tgaaattgaaCCTGTGCTGCCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTTGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGc                                                                                                                                                             <  1:249885/250‑1 (MQ=255)
                                                                                  gaaCCTGTGCTAGCCGCAGAGTTTAACGAGGAGGTGGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCCCGTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGGAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACtatta                                                                                                                                                     <  2:28912/251‑1 (MQ=255)
                                                                                         tgctgcCGGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGACTCTGTCCCAGCACGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGGTTGCCGAAGTGGATGGTAATGGACCGCTGTTCACAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCg                                                                                                                                              >  2:452474/1‑251 (MQ=255)
                                                                                         tgctgcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGGCGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAGAACGACCTGGCCTTTTATAGCCACTATTCTCAGCCg                                                                                                                                              >  1:131538/1‑251 (MQ=255)
                                                                                            tgcCCGCAGAGTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTg                                                                                                                                           >  1:457083/1‑251 (MQ=255)
                                                                                               ccGCAGAGTTTAAAGCTGAGGGTGATGTTGCCGCTTTTATGGCGCATCTGTCACAGCTCGACGGTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGGACAGTTTTGCGATATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGACGAAGTGGATGGTAATGATCCGCTGTTCCAAGGG‑ACAAATGGCGCAAACTCCCTGGCCTTTTTGAGCCACTATTAACAGcagactgcg                                                                                                                                        >  2:215774/1‑243 (MQ=255)
                                                                                                   agagTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTgg                                                                                                                                    <  1:173796/251‑1 (MQ=255)
                                                                                                   agagTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTgg                                                                                                                                    >  1:475592/1‑251 (MQ=255)
                                                                                                   agagTTTAACGCTGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCGGGCCTTTTATAACCACTATTATCAGCCGCGGCCGGTgg                                                                                                                                    >  2:100320/1‑251 (MQ=255)
                                                                                                               tgAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AACAATGGCGAAAACGCCATGGCCTTTTCTAGCCACTATTATCAGCCGCTGCCGTTGGTGCTgcgggggt                                                                                                                        >  2:245993/2‑246 (MQ=255)
                                                                                                                         gTTGCCGCTTTTATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCggg                                                                                                              <  1:59206/251‑1 (MQ=255)
                                                                                                                                 ttttATGGCGAATCTGTCACGGCTCGCCGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGTAGGTGGTGCCTGCCGCGTGAAGATTGCCGACGTGGGTGGTAATGATCCGCTGTTCAAAGGG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGGGCTGCGCgggtaagatcggaaaggcgtcg                                                                                                      >  2:286157/1‑231 (MQ=255)
                                                                                                                                 ttttATGGCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGATAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGata                                                                                                                       <  1:286157/234‑1 (MQ=255)
                                                                                                                                       ggCGAATCTGTCACAGCTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTAc                                                                                                  >  1:110531/1‑249 (MQ=255)
                                                                                                                                                     agcTCGACGATCTCGTTACCGCGCGCGTGGCGAAGGCCCGTGATGCAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAGGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGTTATGGTGCGGGCAATgcagttccatctgcggtggtcttt                                                                                  >  2:93634/3‑228 (MQ=255)
                                                                                                                                                           acgaTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCACTATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGTTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGAt                                                                            >  1:278142/1‑251 (MQ=255)
                                                                                                                                                           acgaTCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGAt                                                                            <  2:457083/251‑1 (MQ=255)
                                                                                                                                                                        cgcgcgcgTGGCGAAGGCCCGGGATGAAGGAAAAGTTTTGCGCTATGTGGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTAccc                                                               <  2:234362/251‑1 (MQ=255)
                                                                                                                                                                             gcgTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCa                                                          <  2:110531/251‑1 (MQ=255)
                                                                                                                                                                                gggCGAAGGCCGGTGATGAAGGAAAAGTTTTGCGTTATGTGGGCAATATTGATGAAGATGGGGCCTGCCGCTTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAAGGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTAACAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTTCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATgg                                                       <  2:26808/250‑1 (MQ=255)
                                                                                                                                                                                 ggCGAAGGCCCGTGAGGAAGGAAAAGTGTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATgg                                                       <  1:320511/250‑1 (MQ=255)
                                                                                                                                                                                      agggccgggatgaAGGAAAAGTTTTGAGCTATGTTGGCAATCTTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCTGCTGTTCACAGTG‑AAAAATGGCGAAAACGCCCTGGACTTTTATAGCCACTATTATCAGCCGCAGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGAAGTTACCGCTGCCGGTGTCTTTGCCGAGCTGATACGGACCCTCTCATGGCAGtt                                                  >  2:156832/9‑250 (MQ=255)
                                                                                                                                                                                              gatgaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGCAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGGGTCTg                                         >  2:138113/1‑251 (MQ=255)
                                                                                                                                                                                                 gaAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGCAGATGGTGCCTGCCGCGTGAAGATTGCAGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑CAAACTGGCGAAAAAGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGAGGGCAATGACGGTCCAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGGATCTGAc                                       >  2:456506/1‑250 (MQ=255)
                                                                                                                                                                                                   agtgaactttggtcggtATGTTGGAATTATTGTTGAATATGGTGCCGGCCGCGGGAATATGGCTGAAGTGGATGGTAATGATCCGCTGTT‑AAAATGTAAAAAGTACGAAACGGCGCTGTCCTGTTATGGCCACTATTATCAGCAGCTGCCGTTGGTTCTGCTCTTGTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATATGCTACGTACCCTCTCATGGACGTTAGGAGTCTGACATg                                    <  2:12012/235‑1 (MQ=255)
                                                                                                                                                                                                         aGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACGTGGTTaa                               >  2:85710/1‑250 (MQ=255)
                                                                                                                                                                                                          gTTTTGCGCTATGTTGGCAATATTGATGAATATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGAc                                                                                                       >  1:26128/1‑177 (MQ=255)
                                                                                                                                                                                                          gTTTTGCGCTATGTTGGCAATATTGATGAATATGGCGTCTGCCGCGTGAAGATGGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAAGGGCGAAAACGCCCTGGCCGTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGAc                                                                                                       <  2:26128/177‑1 (MQ=255)
                                                                                                                                                                                                                     tGTTGGCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTACCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGACATGGTTAAAGTTTATGCCCCg                  >  1:367696/1‑251 (MQ=255)
                                                                                                                                                                                                                         ggCAATATTGATGAAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCtt              <  1:408043/251‑1 (MQ=255)
                                                                                                                                                                                                                                    tgaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAat   >  1:486221/1‑251 (MQ=255)
                                                                                                                                                                                                                                     gaAGATGGTGCCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGGAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAat   >  1:158633/1‑250 (MQ=255)
                                                                                                                                                                                                                                      gaGATGGTGCCTGCCGCGTGACGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGCTATGGTGCGGGCAATGACGTTACAGCTGCCGGGGTCTTTGCCGATCTGCTACGTACCCTCTCATGGACGTTAGGCGTCTGACATGGTTAACGTTTATGCCCCGGCTTCCAGTGCCAata  >  2:123219/2‑250 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                 
TTCTGGTATGGATGTGGCGCGTAAACTATTGATTCTCGCTCGTGAAACGGGACGTGAACTGGAGCTGGCGGATATTGAAATTGAACCTGTGCTGCCCGCAGAGTTTAACGCCGAGGGTGATGTTGCCGCTTTTATGGCGAATCTGTCACAACTCGACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATA  >  NC_000913/2358‑2837

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