Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,121,033 T→C I83V (ATT→GTT)  glcA ← glycolate/lactate:H(+) symporter GlcA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,121,0330TC98.0% 165.1 / NA 49I83V (ATT→GTT) glcAglycolate/lactate:H(+) symporter GlcA
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base C (18/30);  total (19/30)
Fisher's exact test for biased strand distribution p-value = 3.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

GACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACATTTGGGTCCAGGT  >  NC_000913/3120799‑3121273
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
gACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCAACGGCGCACCAAAGCCAGCCGCTCATTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGCACTGCCCGCTGGCAGCGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGccct                                                                                                                                                                                                                                  >  2:547820/1‑249 (MQ=255)
cacccacccccccCGCTTATAACGGTTTGAAGCCCAGGCCCCCCAGCAGCCCACCGGTAATCGCCACCGGCGCACCAAATCCAGCCCCTCCTTCAAGCAACGCACCAAAGGAGAAACCAATCAGTAACCCCTGCAAACGCTGATCGCCGGTGAGGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCAt                                                                                                                                                                                                                                  <  2:440554/244‑1 (MQ=255)
   acaGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCCACGCACCCACGGAGAAACCAATCCGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACCGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAAt                                                                                                                                                                                                                               >  1:582510/1‑251 (MQ=255)
     cgCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACGGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATcc                                                                                                                                                                                                                             <  1:403434/250‑1 (MQ=255)
         ccgccgCGTATACCGGTTGGAACCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCataa                                                                                                                                                                                                                         <  1:43869/251‑1 (MQ=255)
             cgcgTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCataaataa                                                                                                                                                                                                                     <  1:199445/251‑1 (MQ=255)
                  ataACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGCACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGcc                                                                                                                                                                                                                 >  2:474647/1‑250 (MQ=255)
                    aaCGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACCGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCTGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATa                                                                                                                                                                                                              <  1:547820/251‑1 (MQ=255)
                     aCGGTTTGACGCCCAGGCCCACCAGCAGCGCACCGGTAACCGCCCCCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACTCTGATCGTCGGTGATGTAGATAACCGAGCGGCGGATAATAACGAACTGCCCGCGGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAg                                                                                                                                                                                                             <  2:311825/251‑1 (MQ=255)
                       ggttcgacgccaaGGCCCACAAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACAAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCAGCGGATAATCTCGAACTGCCCGCTGGCAACGGGTAATTTATACAGGAACACCGCCGCGACAATAACCCACGATATTGGCCATAATCCATAAATaaagccaaagcc                                                                                                                                                                                                           >  2:23706/13‑251 (MQ=255)
                        gTTTGAAGCCCAGGCCCCCCAGCAGCGCACCGGTAACCGCCACCGGCGCACCAAATCCAGCCGCCCCTTCCAGCAACGCACCAAAGGAGAAACCAACCAGTAACCCCTGAAAACGCTGATCGCCGGGGATGGAGATAACCGAGTTGCGGATAATATCGAACCGCCCGCGGGCAACGGTTAATTTATACAGGAAAACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGccc                                                                                                                                                                                                          <  2:100476/251‑1 (MQ=255)
                            gAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACCCCTGCAAACGCTGATCGTCGGTGATGGAGATACCCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCgcag                                                                                                                                                                                                      <  1:584301/251‑1 (MQ=255)
                                cccAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCgcagcagc                                                                                                                                                                                                  >  1:514336/1‑251 (MQ=255)
                                  aaGGCCCCCCAGCAGCGCCCCGAAAAGCGCCACCGGCGCACAAACGCCAGCCCCTCCTCCCAGCAACGCACCAAAGGATAAACAAATCAGAAACCCCTGAAAACGCTGATCGGCGGTGAGGGAGATAACCGATCTGGGAATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAAGAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCaa                                                                                                                                                                                                <  2:296715/250‑1 (MQ=255)
                                     ccccACCAGCAGCTCACCGGTAACCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATg                                                                                                                                                                                             <  1:417679/250‑1 (MQ=255)
                                         ccaacaagcGCACCGATAACCCCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAAAAACCAATCAGTAACCCCGTCAAACGCTGACCGCCGGTAATGGAGAAACCCGAGCTGCGGATAATATCGAACTGCCCGCTGGCACCGGTTAATTTAACCAGGAACACCGCCGCGACAATACCCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCat                                                                                                                                                                                         <  2:2150/245‑1 (MQ=255)
                                          ccaGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGAGGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCata                                                                                                                                                                                        <  1:243741/251‑1 (MQ=255)
                                          ccaGCAGCCCACCGGTAACCTCCACCGGCCCACCAAACCCAGCCGCCCCTTCCAGCAACGCACCAAAGAAGAAACCAATCAGTAACCCCTGAAAACGCTGAGCGTCGGGGATGGAGATAACCGAGCTGCGGATAATATCGACCTGCCCGCTGTCACCGTTTAATTTATAAAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCata                                                                                                                                                                                        <  2:268316/251‑1 (MQ=255)
                                                gcgcACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCGCCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATc                                                                                                                                                                                  <  1:181248/251‑1 (MQ=255)
                                                          aaccGCCACGGGCGCACCAAAGCCAGCCGCTCCTACCAGCAACGCACCAAAGGAGAAACCAATCAGAACCACCTGAAAACGCGGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATAACGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTaa                                                                                                                                                                        <  2:150875/248‑1 (MQ=255)
                                                            ccGCCACCGGCACAACAAAGCCAGCCCCTCCTTCAAGCAACGCACAAAAGGAGAAACCAATCAGTAACACCTGCAAACTCGGATCGTCGGTGATGGAGATAACCGAGCTGCGAATAATATCGAACTGCCCGCTGGCACTGGTTAATTTAAACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAg                                                                                                                                                                      <  2:153880/250‑1 (MQ=255)
                                                                         cccAAACCCAACCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACCCCTGCAAACGCTGATGGTCGGTGATGGAGATACCCGATCTGCGGATAATATCAACCTGCCCGCTGGCAACGGTAAATTTATACAGAAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCaa                                                                                                                                                         <  2:513921/250‑1 (MQ=255)
                                                                              agccagccGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAgg                                                                                                                                                    <  1:552747/251‑1 (MQ=255)
                                                                                    ccGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGgatagat                                                                                                                                              <  1:239500/251‑1 (MQ=255)
                                                                                             ccAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAgg                                                                                                                                     >  2:20320/1‑251 (MQ=255)
                                                                                                 cAACGCACCAAAGGAAAAACCAATAATTACCCCCTGCAAACGCTGAGCGTCGTTGAGGGAGATAACCAAGCTGAGGATAATATCGAACTGCTCGCTGCCACCGGTTATTTTATACAGATACACCGCCGCGACAATAACCCACGCTATTGTCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCTGGATAGATAATATAAGGGTTa                                                                                                                                 <  2:582492/251‑1 (MQ=255)
                                                                                                            aGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCg                                                                                                                      >  2:172885/1‑251 (MQ=255)
                                                                                                                                   tGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCTGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATAc                                                                                               >  1:421878/1‑251 (MQ=255)
                                                                                                                                          gccgagcgtcgGGGATGGAGATAACCGAGCTGCGAATAAGACCGAACTGCCCGCGGGCAACGGTTACTTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAAGATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGt                                                                                        <  2:97408/245‑1 (MQ=255)
                                                                                                                                              aTCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCa                                                                                     >  1:163330/1‑250 (MQ=255)
                                                                                                                                                    ggTGATGGAGATAACCGAGCTGCGGATAATATCGAAATGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCAGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGGTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACga                                                                                 >  2:353725/1‑248 (MQ=255)
                                                                                                                                                      tGATGGAGATAACCGAGCTGAGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTAAGACGTAATACCGCGAGAGCAACgaaga                                                                              >  2:423277/1‑249 (MQ=255)
                                                                                                                                                      ggATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGTTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACgaagaag                                                                            <  1:474647/250‑1 (MQ=255)
                                                                                                                                                                   cGAGGTGCGGATAATATCGAACGTCCCGCTGGCAAGGGTTAATTTATACAGCAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAAGGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCGCCAGCGACATGTCCTTGCAGACGTAATGCCGCGAGTGCATGGAAGAATAATATTATCGGGa                                                               <  2:582510/251‑1 (MQ=255)
                                                                                                                                                                                    cGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCg                                              >  1:534489/1‑251 (MQ=255)
                                                                                                                                                                                          gcacGCTGGAAACGGTTAATTTATACAGGAACCCCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATGGCTCCAGCGACATGTCCTTTCAGTCGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGCCCAGAGCGGATAg                                         <  2:398349/247‑1 (MQ=255)
                                                                                                                                                                                              gCTGGCAACGGTTATTTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAAGCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCa                                    <  1:518811/251‑1 (MQ=255)
                                                                                                                                                                                                 ggCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGtc                                 >  1:271353/1‑251 (MQ=255)
                                                                                                                                                                                                  gCAACGGTTAATTTATAAAGGAACACCGCCGCGCCAATAATCCACGCTCTTGGCCATAATCCATAAATAAAGCCATAGCCCGCACAAGCAAATGCCATCTCAATCGGCATTTTAAAGGCGCATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTCTCAGACGAAATACAGCGAGTGCAACGAAGACGAATATTATCGGGATCAGGGCGCACAGAGCGGTTAGCCCCAGtcc                                >  2:402314/1‑251 (MQ=255)
                                                                                                                                                                                                       ggTTAATTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTccc                            >  1:505670/1‑250 (MQ=255)
                                                                                                                                                                                                          gaatgtataCAGGAACTCCGCCGCGACAATAACCCACGCTAGTGGCCATAATCCATAAATAAAGCCAGAGCCCGCAGCAGCAAATGCCATATGAATCGGCATTTTAAAGGCGAATAGTGCAATCAGGATAGATAATATAAGGGTTATTGCTTCAGCGACATGTCGTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATAGGATGGGGAGCAGGGCGACCAGAGCGTAGATCCCCAGTCCGCCCATCg                        <  2:237687/246‑1 (MQ=255)
                                                                                                                                                                                                           aaTTTATACAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCACCGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCgg                       >  1:400181/1‑251 (MQ=255)
                                                                                                                                                                                                                  aCAGGAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATACAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTCATACCGCGAGTGCAACgaag                                                                               >  2:342412/1‑188 (MQ=255)
                                                                                                                                                                                                                     ggAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACAtt             <  1:246033/251‑1 (MQ=255)
                                                                                                                                                                                                                      gAACACCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACAttt            >  1:384512/1‑251 (MQ=255)
                                                                                                                                                                                                                       accacCGCCGCGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACATTTg           <  2:426573/249‑1 (MQ=255)
                                                                                                                                                                                                                          acctccgcGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACATTTGGGt        <  2:495284/247‑1 (MQ=255)
                                                                                                                                                                                                                              ccgcGACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACATTTGGGTCCAg    <  2:534489/251‑1 (MQ=255)
                                                                                                                                                                                                                                gcgACAATAACCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCACGTCCTCCCATCGGCATATACATTTGGGTCCAGGt  >  1:106433/1‑251 (MQ=255)
                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                
GACACAGCCCCGCCGCGTATAACGGTTTGAAGCCCAGGCCCACCAGCAGCGCACCGGTAATCGCCACCGGCGCACCAAAGCCAGCCGCTCCTTCCAGCAACGCACCAAAGGAGAAACCAATCAGTAACACCTGCAAACGCTGATCGTCGGTGATGGAGATAACCGAGCTGCGGATAATATCGAACTGCCCGCTGGCAACGGTTAATTTATACAGGAACACCGCCGCGACAATAATCCACGCTATTGGCCATAATCCATAAATAAAGCCATAGCCCGCAGCAGCAAATGCCATATCAATCGGCATTTTAAAGGCGAATATTGCAATCAGGATAGATAATATAAGGGTTATTGCTCCAGCGACATGTCCTTTCAGACGTAATACCGCGAGTGCAACGAAGAAGAATATTATCGGGATCAGGGCGACCAGAGCGGATAGCCCCAGTCCTCCCATCGGCATATACATTTGGGTCCAGGT  >  NC_000913/3120799‑3121273

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