Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,873,336 C→T V93M (GTG→ATG)  dgoA ← 2‑dehydro‑3‑deoxy‑6‑phosphogalactonate aldolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,873,3360CT97.1% 228.9 / ‑1.3 70V93M (GTG→ATG) dgoA2‑dehydro‑3‑deoxy‑6‑phosphogalactonate aldolase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (38/29);  total (39/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.67e-01

GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATG  >  NC_000913/3873091‑3873580
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
gcccgcccCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCACACGGCGGAtcatttcac                                                                                                                                                                                                                                                 >  1:340125/1‑251 (MQ=255)
           cccagccTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCTAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGTCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTg                                                                                                                                                                                                                                      <  2:144302/246‑1 (MQ=255)
                   gCGTCTAGCCACTGCTCCAGGTTTTCTGGCTGCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGGAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAAc                                                                                                                                                                                                                              <  1:127991/251‑1 (MQ=255)
                       cTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATg                                                                                                                                                                                                                          >  2:150023/1‑251 (MQ=255)
                          tCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCCTGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATg                                                                                                                                                                                                                          >  1:231812/1‑248 (MQ=255)
                                     aGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCAt                                                                                                                                                                                                            >  1:132570/1‑251 (MQ=255)
                                           tCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGc                                                                                                                                                                                                      >  1:326060/1‑251 (MQ=255)
                                                cGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGc                                                                                                                                                                                                                <  2:142837/236‑1 (MQ=255)
                                                cGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGc                                                                                                                                                                                                                >  1:142837/1‑236 (MQ=255)
                                                  ccACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCa                                                                                                                                                                                               <  1:251697/250‑1 (MQ=255)
                                                   cACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCTCTATGGATATTGGGCGTAACGATGAGCt                                                                                                                                                                                                                      <  2:223199/227‑1 (MQ=255)
                                                   cACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCt                                                                                                                                                                                                                      >  1:223199/1‑227 (MQ=255)
                                                      gccgccAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCg                                                                                                                                                                                            >  1:377725/1‑250 (MQ=255)
                                                        cgccAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGAcc                                                                                                                                                                                         >  2:93554/1‑251 (MQ=255)
                                                         gccAACGGCAAAGACTGCGATGTCCGATGGCAATACAGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCt                                                                                                                                                                                        >  1:359031/1‑251 (MQ=255)
                                                            aaCGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGCTGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGtt                                                                                                                                                                                     >  2:22751/1‑251 (MQ=255)
                                                             aCGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTc                                                                                                                                                                                    >  2:405019/1‑251 (MQ=255)
                                                                  aaaGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGtt                                                                                                                                                                               >  1:199498/1‑251 (MQ=255)
                                                                    aGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTc                                                                                                                                                                             >  1:210593/1‑251 (MQ=255)
                                                                                       cAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGCCGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCaa                                                                                                                                                          <  1:107259/251‑1 (MQ=255)
                                                                                              gCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGc                                                                                                                                                   <  1:168086/251‑1 (MQ=255)
                                                                                                 tttAACGCCTTCATGTATTGCGCACCAAAAGCCGATGACGGAAATATTTTCCGCGCCTGCGCGCCCGCCTCGAGCGCGGTAACGGCTTCGGTCGCCGTCGCGCCGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCCGATCTTTTCACTATGGATATTCGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGTACTGTTCAGGTTTCACTACCGTCCCTGCGCCACTCAACGCCtt                                                                                                                                                >  2:43405/1‑251 (MQ=255)
                                                                                                     aCGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCg                                                                                                                                            >  1:345471/1‑251 (MQ=255)
                                                                                                      cGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGTGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGc                                                                                                                                           <  1:112637/251‑1 (MQ=255)
                                                                                                        cTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCg                                                                                                                                         >  1:61727/1‑251 (MQ=255)
                                                                                                             aTGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGc                                                                                                                                    >  1:305358/1‑251 (MQ=255)
                                                                                                               gTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCa                                                                                                                                   <  1:379307/250‑1 (MQ=255)
                                                                                                                tATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGcc                                                                                                                                                                                                               <  2:19335/173‑1 (MQ=255)
                                                                                                                tATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGcc                                                                                                                                                                                                               >  1:19335/1‑173 (MQ=255)
                                                                                                                 aTTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGTCGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCTCCAATCAACGCCTTGTCGCCGTACGCATCa                                                                                                                                <  1:308420/251‑1 (MQ=255)
                                                                                                                   ggCGGACGAAGAGGCGATGACGGAAATATTTTGAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTGGCCGTCGCGCAGCGGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAAc                                                                                                                              <  2:276917/250‑1 (MQ=255)
                                                                                                                     cGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGa                                                                                                                            <  1:15538/251‑1 (MQ=255)
                                                                                                                          cAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGTTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCCCCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACTATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCATGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCg                                                                                                                       <  2:112159/251‑1 (MQ=255)
                                                                                                                          cAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATAAGCTGACACCCCATCCTGGCGAGCGCTTCGACTGGTTCAGGTTTCAGTACCGAACCGGCGCCAATCAACCCCTTTTCGCCGTACGCTTCAACGATGGCg                                                                                                                       >  2:230854/1‑251 (MQ=255)
                                                                                                                             aaGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCACCGATGGCggg                                                                                                                     >  1:231853/1‑250 (MQ=255)
                                                                                                                             aaGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGTCGGGa                                                                                                                    <  2:89656/251‑1 (MQ=255)
                                                                                                                                ccGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGAAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGACTTGTCGCCGTAAGCATCAACGATGGCgggcagg                                                                                                                 >  2:276987/1‑247 (MQ=255)
                                                                                                                                ccGATAACGTAAATATTTTCAGCGTCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGTTCGCCGTCGCGCAGCAGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATCTTGGGCGTAACGATGAGCTGACAGCCCACCCTGGCGAGCGCATTGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCACTCAACGCCTTGTCGCCGTAAGCATCAACGATGGCGGGAAtt                                                                                                                 <  2:155316/251‑2 (MQ=255)
                                                                                                                                 cGATGACGGAACTATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGAAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATAGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTAAGCATCAACGATGGCggggatac                                                                                                                >  2:234220/1‑246 (MQ=255)
                                                                                                                                  gATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCt                                                                                                               <  1:22751/251‑1 (MQ=255)
                                                                                                                                   aTGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGCAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCtt                                                                                                              >  1:266603/1‑251 (MQ=255)
                                                                                                                                     gACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCttt                                                                                                             <  1:110202/250‑1 (MQ=255)
                                                                                                                                         gagatatTTTCAGCGCCTGCGAGCCCGCTTCGAGCTCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATAATTTCACTATGGATATTGGGAGTAACGATGAGCTGACAGCCCATCCAGGCGAGCGCATCAACCTGATCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGCCGCCGGACGCATCAACGATGGCGGGAATGGttttgccc                                                                                                        >  2:322327/4‑246 (MQ=255)
                                                                                                                                         gAAATATTTGTAGCTCCGGCGCGCCCGCTTCGAGCGCGGTAAAGGCCTCGGTCGCCGTCGCAGAGCCGGGGCAGACGGTCCTGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGGGCGCATCGACCTGGTCAGGTTTCAGGACCGTACCCGCGGCAATTATCTCCGTGTCGCCGTACGCATCAAGGATGTCGGGAATGCTTTGCTcc                                                                                                        <  2:176170/251‑1 (MQ=255)
                                                                                                                                                 ttCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTgg                                                                                                >  2:330394/1‑251 (MQ=255)
                                                                                                                                                       gcCTGCGCGCCCGCTTCGAGCGCGGTAAGGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGGATGCTTTGCTCCCATTGTGGGGAAtt                                                                                          >  2:296781/1‑251 (MQ=255)
                                                                                                                                                       gcCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGGTTTGCTCCCCTTGTGGGGAAtt                                                                                          >  2:234340/1‑251 (MQ=255)
                                                                                                                                                                 ccGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGTGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAtt                                                                                <  2:319928/251‑1 (MQ=255)
                                                                                                                                                                  cGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAttt                                                                               >  1:359263/1‑251 (MQ=255)
                                                                                                                                                                           cgcgGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGt                                                                      >  2:75811/1‑251 (MQ=255)
                                                                                                                                                                                tAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAAc                                                                 >  1:281464/1‑251 (MQ=255)
                                                                                                                                                                                 aaaGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGTAATTCAGCGGGATTTCAACCGCGGCGAAcc                                                                <  1:123879/251‑1 (MQ=255)
                                                                                                                                                                                 aaaGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAAcc                                                                >  1:347616/1‑251 (MQ=255)
                                                                                                                                                                                     gCTTCGGTCTCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGc                                                            <  2:346/251‑1 (MQ=255)
                                                                                                                                                                                          ggTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCaa                                                       <  2:61727/251‑1 (MQ=255)
                                                                                                                                                                                                          ccGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGCATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATg                                       >  2:89507/1‑251 (MQ=255)
                                                                                                                                                                                                          ccGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACAtc                                       >  2:115851/1‑250 (MQ=255)
                                                                                                                                                                                                           cGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATgc                                      >  1:235287/1‑251 (MQ=255)
                                                                                                                                                                                                            ggggCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATgcc                                     >  2:246153/1‑250 (MQ=255)
                                                                                                                                                                                                                 aGACGGGCATGCCGTAGCCCCCCGCACGGCGGATCTGTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGTGGGAATGCTTTGCTCCCATTGTGGTGAATTCAGCGTGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAga                                <  2:39142/251‑2 (MQ=255)
                                                                                                                                                                                                                  gACGGTCATGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCAACGCGCCAACATGCGCCAgcc                               >  1:259569/1‑250 (MQ=255)
                                                                                                                                                                                                                         aTGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgtc                        <  2:2984/251‑1 (MQ=255)
                                                                                                                                                                                                                         aTGCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgtc                        <  2:305358/251‑1 (MQ=255)
                                                                                                                                                                                                                           gCCGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCtcgtcg                       <  2:132570/250‑1 (MQ=255)
                                                                                                                                                                                                                            ccGTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCggg                     >  1:181239/1‑251 (MQ=255)
                                                                                                                                                                                                                              gTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCg                   <  1:258574/251‑1 (MQ=255)
                                                                                                                                                                                                                              gTAGCCCACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCg                   <  1:171848/251‑1 (MQ=255)
                                                                                                                                                                                                                                aGCCCACCGCACGGCGGATCATTTCCCTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGTAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGCACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGt                  >  1:4411/1‑250 (MQ=255)
                                                                                                                                                                                                                                   ccACCGCACGGCGGATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAAt               >  2:122186/1‑250 (MQ=255)
                                                                                                                                                                                                                                           cggcggATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGAGTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCaa      <  1:234340/251‑1 (MQ=255)
                                                                                                                                                                                                                                               ggATCATTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATg  <  2:317381/251‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    
GCCCGCCCCTGCACAACCTGCGTCTATCCACTGCGCCAGGTTTTCTGGCGTCACGCCGCCAACGGCAAAGACTGCGATGTCCGATGGCAATACCGCTTTTAACGCTTTGATGTATTGCGGACCAAAAGCCGATGACGGAAATATTTTCAGCGCCTGCGCGCCCGCTTCGAGCGCGGTAAAGGCTTCGGTCGCCGTCGCGCAGCCGGGGCAGACGGTCATGCCGTAGCCCACCGCACGGCGGATCACTTCACTATGGATATTGGGCGTAACGATGAGCTGACAGCCCATCCTGGCGAGCGCATCGACCTGTTCAGGTTTCAGTACCGTACCTGCGCCAATCAACGCCTTGTCGCCGTACGCATCAACGATGGCGGGAATGCTTTGCTCCCATTGTGGGGAATTCAGCGGGATTTCAACCGCGTCGAACCCGGCGTCAATCACCGCGCCAACATGCGCCAGCGCCTCGTCGGGCGTAATACCGCGCAAAATG  >  NC_000913/3873091‑3873580

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