Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,860 T→A A508A (GCT→GCA)  thrA → fused aspartate kinase/homoserine dehydrogenase 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,8600TA97.1% 214.7 / NA 72A508A (GCT→GCA) thrAfused aspartate kinase/homoserine dehydrogenase 1
Reads supporting (aligned to +/- strand):  ref base T (1/1);  major base A (27/41);  minor base C (2/0);  total (30/42)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTAC  >  NC_000913/1615‑2103
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
aTTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCt                                                                                                                                                                                                                                                 <  2:648993/251‑1 (MQ=255)
 gtggcgacatggcgtatggatctggtGAACGCTCAAACGCTGTCGTGGTAAATAACGATGATGCGACCACTGGTGTGCGCGTTTCTCATCAGATGCTGTTCAATCCCGATCAGGTTATCGATGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGCAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTc                                                                                                                                                                                                                                                <  2:324984/226‑1 (MQ=255)
     cGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGCTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCa                                                                                                                                                                                                                                            <  2:301528/251‑1 (MQ=255)
        cATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑Tg                                                                                                                                                                                                                                        >  1:58181/1‑251 (MQ=255)
             cataGGGATCTTCGGAAAGCTCAATCTCTGTCGGGGTAAATAACGAAGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCAt                                                                                                                                                                                                                                   <  2:125033/248‑1 (MQ=255)
                 gggATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGCTCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGGAGCGTCCACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACCCGAAGGCAATACCCACCAA‑TGTGCATGGcc                                                                                                                                                                                                                               >  2:463057/1‑251 (MQ=255)
                  ggATCTTCTGAAAGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGATATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCt                                                                                                                                                                                                                              >  1:112668/1‑251 (MQ=255)
                        tctGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATc                                                                                                                                                                                                                         <  1:253589/250‑1 (MQ=255)
                           gAACGCTCAATCGCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGTCTAAATCTg                                                                                                                                                                                                                       <  1:107304/249‑1 (MQ=255)
                            aaCGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGa                                                                                                                                                                                                                     >  1:666637/1‑250 (MQ=255)
                             aCGCTCAATCTATGTCGTGGGAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATACTGGTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGAGACGGTGGCGTTGACGGTACGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGATGAAGAATAAACATATCGCCTTAAGTGTCTGCGGTGTTGACAACTCGACGGACCTGCTCACCAA‑Tgggcagggcctacaagtggta                                                                                                                                                                                                                    >  2:420680/1‑230 (MQ=255)
                                  cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                                                                                <  1:227325/249‑1 (MQ=255)
                                  cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTgg                                                                                                                                                                                                              <  1:256849/251‑1 (MQ=255)
                                  cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTg                                                                                                                                                                                                               <  1:51793/250‑1 (MQ=255)
                                  cAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCGGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTTTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACt                                                                                                                                                                                                                <  2:374563/249‑1 (MQ=255)
                                          gTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAAAATAAACATATCGACTTACGTGTCTGCGGGGTTTACAACTCGAAGGCACTGCTCACACA‑TGTGCATGGCCTAAATCTGGAAAACTGCCAGga                                                                                                                                                                                                         >  2:576554/1‑248 (MQ=255)
                                               ggTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCTACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTgg                                                                                                                                                                                                 <  1:626377/251‑1 (MQ=255)
                                                    aTAACGATGAAGCGACCACTGGCATGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGGGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTATCTGCGGTGTTGCCAACCCGAAGTACCTGCTCACCAAATGGAC‑TGGACCAAATCTGGAAAACTGGCAGgaaaaaacggcccc                                                                                                                                                                                             >  2:444113/1‑239 (MQ=255)
                                                      aaCGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGGACTGGCGCAAGc                                                                                                                                                                                          <  1:501096/251‑1 (MQ=255)
                                                        cGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCa                                                                                                                                                                                        <  1:290569/251‑1 (MQ=255)
                                                             atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAt                                                                                                                                                                                                                          <  1:3764/212‑1 (MQ=255)
                                                             atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAt                                                                                                                                                                                                                          >  2:3764/1‑212 (MQ=255)
                                                             atgCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAagag                                                                                                                                                                                   <  1:268555/251‑1 (MQ=255)
                                                                  aCCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATTGCCTAAATCTGGAAAACTGGCCGGAAGAACTGGCGCAAGCCAAAGAGCCGtt                                                                                                                                                                              >  2:20377/1‑251 (MQ=255)
                                                                     aCTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTaa                                                                                                                                                                           <  1:20377/251‑1 (MQ=255)
                                                                        gtcgttctcgTTACTCTGCAGATGCTGTTCAAGTCCGATCAGGTTATCGACGTGTTTGTGATTGGCGTCGGTGTCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACTTGTCTTCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCGAAATCTTGAAAACTGGCAGGAAGAACTTGCGTAAGCCAAATAGCCGTTTAAt                                                                                                                                                                          <  2:34004/241‑1 (MQ=255)
                                                                           gTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCgg                                                                                                                                                                     <  1:249673/251‑1 (MQ=255)
                                                                                 gTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTAGCGTCGGTGGCGTTGGCGGTGCGATGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACAAATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAAACTGGAAAACTGGCAGGAAGAACTGGACCAAGCCAAAGAACCGTTTAATATCGGGCGCtt                                                                                                                                                               >  2:477947/1‑251 (MQ=255)
                                                                                     cTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATCTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTAGGGCGCTTAAtt                                                                                                                                                           >  2:163745/1‑251 (MQ=255)
                                                                                       catctGATGCTGGTCAATTCCGATCAGGTTATCGAAGTGTTTGGGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAATAATAAACATGTCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCGTACCAA‑TGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCTATAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                          <  2:655808/250‑1 (MQ=255)
                                                                                       catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGGCGGTGGCGTTGGCGGTGCGTTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTc                                                                                                                                                          <  1:388932/250‑1 (MQ=255)
                                                                                       catcaGATGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                                                                                         <  1:247855/251‑1 (MQ=255)
                                                                                       catcaGAGGCTGTTCAATACCGATCAGGTTATCGAAGTGGTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCg                                                                                                                                                         <  1:617089/251‑1 (MQ=255)
                                                                                           aGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTc                                                                                                                                                     <  1:609608/251‑1 (MQ=255)
                                                                                            ggtgtggTTCAATACCGAGCAGGTTTTCGAAGTGGTTGTGATTGGCTGCGGTGGCGTGGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACGGGCAGGAAGAACTGGCGCAAGCCAAAGAGCTGTTTAATCTCTGGCGCTTAATTCGCCTCg                                                                                                                                                    <  2:400836/245‑1 (MQ=255)
                                                                                                         aCCGATCAGGTTCTCGTAGTGTTTGTGATTCACGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAGAGCTGGCTGGGGAATAACCATATCGGCTTACGTGTCTGCGGTGTTGCCAACTCGCAGGCACTGCTCCCCAC‑TGAGCATGGCCTAAATCTGGaa                                                                                                                                                                                                                    >  2:73224/1‑174 (MQ=255)
                                                                                                         aCCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGaa                                                                                                                                                                                                                    <  1:73224/174‑1 (MQ=255)
                                                                                                                aGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGAGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATctgctg                                                                                                                                >  1:41477/1‑251 (MQ=255)
                                                                                                                   ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGCATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACAAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTAGAGCAAGCAAAAGCGCAGTTTTATATCGGGCGCTTTATTCGCCTCGTGAAAGAATATCATCTGATGAAc                                                                                                                             >  2:29525/1‑251 (MQ=255)
                                                                                                                   ttATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAc                                                                                                                             >  1:439176/1‑251 (MQ=255)
                                                                                                                    tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                                                                                          >  2:125305/1‑205 (MQ=255)
                                                                                                                    tATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAAt                                                                                                                                                                          <  1:125305/205‑1 (MQ=255)
                                                                                                                     aTCGAAGTGTTTGGGATTGGCGGCGGGGGAGTAGGCGGTGCGCTTCTGGTGCAACTGAATCGTCAACAAAGCGTGCTGATGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCATTGCTCACCAA‑TGTGCATGGCCTAAATTTGGAAAACTGGCAGGAAGAACTGGCGCAAGACAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAAcc                                                                                                                            <  2:375278/250‑1 (MQ=255)
                                                                                                                       cGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCg                                                                                                                          <  1:156721/250‑1 (MQ=255)
                                                                                                                            tgtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCCAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGCAAGAATATCATCTGCTGAACCCGGTCAtt                                                                                                                    >  2:563405/1‑251 (MQ=255)
                                                                                                                             gtTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCAttg                                                                                                                   <  1:545447/251‑1 (MQ=255)
                                                                                                                                   gATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACt                                                                                                             >  2:176529/1‑251 (MQ=255)
                                                                                                                                   gATTGGAGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCCCCCCGTGAAAAAATATCATCTTCTGACCCCGGTCATTGTTGACt                                                                                                             >  2:280962/1‑251 (MQ=255)
                                                                                                                                    aTTGGTGTCTGTGGCGTTGGGGGTTCGTTGATGGAGCAACTGAGGCGTCAATAACGCTGGCTGAAGAAGAAACATATCGTCTGACGTGTCTGCTGTGTTGCCAACTCGAAGGCACAGCACACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAAATAGCCGTTTAATCACGGGCGCGTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGGCATTGTTGACt                                                                                                             <  2:543366/250‑1 (MQ=255)
                                                                                                                                                     tGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGACTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGCACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGccg                                                                                           >  2:261741/1‑249 (MQ=255)
                                                                                                                                                                   tgGAGCAACTGAAGCGTCAACAAAGCTGGCTGCAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGt                                                                                                                                                                                                                                       >  2:352448/1‑97 (MQ=255)
                                                                                                                                                                   tgGAGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGt                                                                                                                                                                                                                                       <  1:352448/97‑1 (MQ=255)
                                                                                                                                                                      aGCAACTGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc                                                                          <  1:611862/251‑1 (MQ=255)
                                                                                                                                                                      aGCAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGcc                                                                          <  1:344392/251‑1 (MQ=255)
                                                                                                                                                                        cAACTGAAGCGTCAACAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATgacgg                                                                        >  1:89258/1‑247 (MQ=255)
                                                                                                                                                                            tGAAGCGTCAACAAAGCTGTCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGACCAATATGCCgccctc                                                                    >  2:673986/1‑246 (MQ=255)
                                                                                                                                                                               aGCGTCAACAAATATGGCTGAAGAATAATCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACGGGCGCAAGCCAACGAGACGTTTAATCTCGGGCGCTTAATTAGCCTCGTGAACGAATATCATCTGCTGAACCCGGTCCTTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGGTCAATATGCCGACtttctt                                                                 >  2:378470/1‑247 (MQ=255)
                                                                                                                                                                                          aGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGttt                                                      <  1:373944/251‑1 (MQ=255)
                                                                                                                                                                                                    agaaTAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATTAGGCGAATTCTTGCGCGAAGGTTTCCACGTTGTc                                            >  1:476267/1‑251 (MQ=255)
                                                                                                                                                                                                         aaaCATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGcc                                       <  1:463057/251‑1 (MQ=255)
                                                                                                                                                                                                         aaaCATAACGACTTACGTGTCGGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGcc                                       <  2:163702/251‑1 (MQ=255)
                                                                                                                                                                                                                 cgattgaCGTGTCGTCGGTGGTGCCAACTAGAAGGAACTTCGAACCAA‑TGTGCCTGGCCTAAATCTGGAAAACTTGCAGGAAGAACTGGCGTAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTTAGTCGCCTCGTGAAAGAAGATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCTGAGCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGGCTCGCCGAACaaa                                <  2:289873/244‑1 (MQ=255)
                                                                                                                                                                                                                  gACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGGTGTCACGCCGAACcaaa                               >  1:378867/1‑250 (MQ=255)
                                                                                                                                                                                                                   aCTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCGCTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACaaaa                               >  1:329781/1‑249 (MQ=255)
                                                                                                                                                                                                                     ttACGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGCCTTCTTGCGCGAAGGGTTCCACGTTGTCACGCCGAACAAAAAGGc                           >  1:113085/1‑251 (MQ=255)
                                                                                                                                                                                                                       aCGTGTCTGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGcc                          <  1:110532/250‑1 (MQ=255)
                                                                                                                                                                                                                              tGCGGTGTTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCAACCCGCACAAAAAAGCCAACAccac                  >  2:599376/1‑249 (MQ=255)
                                                                                                                                                                                                                                    gtTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCGGGAAAACTGGCAGGAAGAACTGGCGAAAGCCAAAGTGCCGCTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCAGCTGCCGAACCCGGTCATTGTTGACTGCACTTCCAGCCATGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGTCCAACACCTCGTCGAt            <  2:41477/251‑1 (MQ=255)
                                                                                                                                                                                                                                    gtTGCCAACTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGAt            <  1:176529/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcca                                                                                                    >  1:68169/1‑155 (MQ=255)
                                                                                                                                                                                                                                            cTCGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTcca                                                                                                    <  2:68169/155‑1 (MQ=255)
                                                                                                                                                                                                                                              cGAAGGCACTGCTCACCAA‑TGTGCATGGCCTAAATCTGGAAAACTGGCAGGAAGAACTGGCTCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCTTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATtactac  <  1:132325/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
ATTGTCGCCATTGCTCAGGGATCTTCTGAACGCTCAATCTCTGTCGTGGTAAATAACGATGATGCGACCACTGGCGTGCGCGTTACTCATCAGATGCTGTTCAATACCGATCAGGTTATCGAAGTGTTTGTGATTGGCGTCGGTGGCGTTGGCGGTGCGCTGCTGGAGCAACTGAAGCGTCAGCAAAGCTGGCTGAAGAATAAACATATCGACTTACGTGTCTGCGGTGTTGCCAACTCGAAGGCTCTGCTCACCAA‑TGTACATGGCCTTAATCTGGAAAACTGGCAGGAAGAACTGGCGCAAGCCAAAGAGCCGTTTAATCTCGGGCGCTTAATTCGCCTCGTGAAAGAATATCATCTGCTGAACCCGGTCATTGTTGACTGCACTTCCAGCCAGGCAGTGGCGGATCAATATGCCGACTTCCTGCGCGAAGGTTTCCACGTTGTCACGCCGAACAAAAAGGCCAACACCTCGTCGATGGATTACTAC  >  NC_000913/1615‑2103

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: ‑