Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,760 T→C A808A (GCT→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,7600TC96.4% 159.3 / ‑3.9 55A808A (GCT→GCC) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/1);  new base C (26/27);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.60e-01

ACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTAT  >  NC_000913/2513‑2998
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
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)
                                                          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)
                                                                                   ggcTGCCGCGTGAAGATGGCCGAATGGGATGGTAAGGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAAGGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTTCAGTGCCAATATGAGCGTCg                                                                                                                                                          <  2:163236/249‑1 (MQ=255)
                                                                                               aaGATTGCCGCAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCAGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGACGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGGCGGGGTTGAtgtgg                                                                                                                                              >  2:434451/1‑250 (MQ=255)
                                                                                                     gCCGCAGTGGATGGGAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGCGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCgggg                                                                                                                                        >  2:302608/1‑251 (MQ=255)
                                                                                                       cGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAACATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCAGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGAAGTTACAGCTGACGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATACGAGAGTCGGGTTTGATGTGCTCGGggcg                                                                                                                                      >  2:121437/1‑251 (MQ=255)
                                                                                                                 ggTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCAGTTGGTGCTGCGCGGGTATGGTACGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTAACCTCTCATGGAAGTTAGGCGTCTGCCATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGAGTCGGGTTTGATGGGATCGGGGCGGCGGTGacac                                                                                                                            >  1:297402/1‑251 (MQ=255)
                                                                                                                 cgactcccttccgatctCAAAGTG‑AAAAATGGCGAAAACGCCCTTGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGacac                                                                                                                            <  2:225532/235‑1 (MQ=255)
                                                                                                                    aatgacccgattgtaaaaGTA‑AAAAATGGCGAAAATGCCCTGGCCTTCTATAGCCACTATGATCAGCCGCTGCCGTTGGTAATGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTGCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCt                                                                                                                          <  2:431870/235‑1 (MQ=255)
                                                                                                                                  cAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGacac                                                                                                                            >  1:225532/1‑234 (MQ=255)
                                                                                                                                   aaaGTG‑AAAAATGGCGAAAACGCCTTGGCCTTTTATATCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGc                                                                                                          <  2:168246/251‑1 (MQ=255)
                                                                                                                                   aaaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATgg                                                                                                                                                                                                    <  2:35042/161‑1 (MQ=255)
                                                                                                                                   aaaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATgg                                                                                                                                                                                                    >  1:35042/1‑161 (MQ=255)
                                                                                                                                    aaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                                                                                              >  2:354137/1‑166 (MQ=255)
                                                                                                                                    aaGTG‑AAAAATGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAg                                                                                                                                                                                              <  1:354137/166‑1 (MQ=255)
                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACa                                                                                                                                                                                                       <  1:136578/152‑1 (MQ=255)
                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACa                                                                                                                                                                                                       >  2:136578/1‑152 (MQ=255)
                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTg                                                                                                                                                                         <  2:255686/182‑1 (MQ=255)
                                                                                                                                          aaaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTg                                                                                                                                                                         >  1:255686/1‑182 (MQ=255)
                                                                                                                                           aaaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGgagg                                                                                                   >  1:24703/1‑250 (MQ=255)
                                                                                                                                            aaaTGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCgg                                                                                                      >  1:272209/1‑247 (MQ=255)
                                                                                                                                               tgtcGAAAACGCCCTGTCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCTCGGATATGGTGCGGTCAATGACGTTACAGCTTCAGGTGTCTTTGCCGATCTGCTACTTACCCTCTCATGTAAGTTAGTACTCTCACATTGTTAAAGTTTAt                                                                                                                                                                                         <  2:220924/158‑1 (MQ=255)
                                                                                                                                               tGGCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTAt                                                                                                                                                                                         >  1:220924/1‑161 (MQ=255)
                                                                                                                                                ggCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAg                                                                                              <  1:245993/251‑1 (MQ=255)
                                                                                                                                                 gCGAAAACGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGt                                                                                             <  1:275111/251‑1 (MQ=255)
                                                                                                                                                     aaaCGCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCTGAGATGTAGTCACg                                                                                         <  1:123219/251‑1 (MQ=255)
                                                                                                                                                         gCCCTGGCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTg                                                                                     >  1:226037/1‑251 (MQ=255)
                                                                                                                                                             gggCCTTTTATAGCCACTATTATCAGCCGCTGCCGTGGGTGCTGCGCGGATAGGGTGCGGGCAATGACGTTACAGCTTCCGGTGTCTTTGCCGATCTGTTACGTCCCCTCTCATGGAAGTTAGGAGTCGGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGAGGTGCTCGGGGCGGCGGTGACACCTGTTGAGGGTGCATTGCTCGGAGATTTAGTCACGGTTGAggc                                                                                 <  2:287142/250‑1 (MQ=255)
                                                                                                                                                              ggCCTTTTATAGCCACTATTATCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAggcg                                                                                <  1:1014/251‑1 (MQ=255)
                                                                                                                                                                                 tatCAGCCGTGGCCGTTGGTGCTGCGCGGAGATGGTGCGGGCAATGCCGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTCCCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGTCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTTCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGCCATTCAGcctc                                                              <  2:70531/250‑1 (MQ=255)
                                                                                                                                                                                 tatCAGCCGCTGCCGTTGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGtctc                                                              >  1:3633/1‑250 (MQ=255)
                                                                                                                                                                                    cAGCCGCTGCCGTTGGTGCTGCGCGGATATGGGGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATACGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTcaaca                                                          <  1:452474/251‑1 (MQ=255)
                                                                                                                                                                                            gCCGTTGGGGCTGCGCGGAGAGGGGGCGGGCAATGACGTTACAGCTGCCTGTGTCTGTGCCGATCTGCTACGGACCCCCTCATGGAAGTTAGGATTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCGCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCg                                                    <  2:291734/249‑1 (MQ=255)
                                                                                                                                                                                                 tGGTGCTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCt                                              <  1:93634/250‑1 (MQ=255)
                                                                                                                                                                                                          gcgGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTGCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATaa                                    <  1:394727/251‑1 (MQ=255)
                                                                                                                                                                                                                atGGTTCGGGCCATGACCTTACAGCTGCTGGTGTCTTTTCCGATCTGCTAGGTACCCTCGCATGGAAGTTAGGAGTCTGACAAGGTTAAAGTGTATGCCCCGGCTTCCATTGCAAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCAGAGATGTATTCACGGGTGACGCGGCAGAGACATTCAGTCGGAACATCTGCGGACGGTTTGCCGATAAGCTGcc                              <  2:253128/251‑1 (MQ=255)
                                                                                                                                                                                                                    tGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGGGGCGGCAGAGACATTCAGTCTCAAAAACCTCGGACGCTTTGCCGATAAGCTGCCGTc                           >  1:223455/1‑250 (MQ=255)
                                                                                                                                                                                                                           aatggcgttagagttgtcgttgtCTTTGCCGATCGGCTACGGACCCTCTCATGGAAGTGAGGAGGCTGACATGGTTAAAGTGTAGGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTAGGGGCGTCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATACAGTCCCAACAACCTCGGACGCTTTGCCTATAAGCTGCCGTCAGAACCAc                   <  2:416639/231‑1 (MQ=255)
                                                                                                                                                                                                                            aTGACGTTACAGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACg                  >  1:136238/1‑251 (MQ=255)
                                                                                                                                                                                                                                      aGCTGCCGGTGTCTTTGCCGATCTGCTACGTACCCTCTCATGGAAGATAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGGTGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGGAGCTTTGCCGGTACTCTGCCGTCAGAACCACggggcaatat         >  2:453675/1‑243 (MQ=255)
                                                                                                                                                                                                                                            cGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGCGATGTAGTCACGGTTGAGGCGGCAGAGACCTTCAGTCTCAACAACCCCGGACGCTTTTCCGGAAAGCAGCCGTCAATAACACGGGAAAATAACGTTTAt  >  2:255399/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                              
ACGATCTCTTTGCCGCGCGCGTGGCGAAGGCCCGTGATGAAGGAAAAGTTTTGCGCTATGTTGGCAATATTGATGAAGATGGCGTCTGCCGCGTGAAGATTGCCGAAGTGGATGGTAATGATCCGCTGTTCAAAGTG‑AAAAATGGCGAAAACGCCCTGGCCTTCTATAGCCACTATTATCAGCCGCTGCCGTTGGTACTGCGCGGATATGGTGCGGGCAATGACGTTACAGCTGCCGGTGTCTTTGCTGATCTGCTACGTACCCTCTCATGGAAGTTAGGAGTCTGACATGGTTAAAGTTTATGCCCCGGCTTCCAGTGCCAATATGAGCGTCGGGTTTGATGTGCTCGGGGCGGCGGTGACACCTGTTGATGGTGCATTGCTCGGAGATGTAGTCACGGTTGAGGCGGCAGAGACATTCAGTCTCAACAACCTCGGACGCTTTGCCGATAAGCTGCCGTCAGAACCACGGGAAAATATCGTTTAT  >  NC_000913/2513‑2998

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