Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,986,009 A→C 8.1% intergenic (+37/+153) ilvE → / ← PP_3512 branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,986,0090AC 8.1% 499.9 / 25.1 174intergenic (+37/+153)ilvE/PP_3512branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (92/67);  new base C (9/6);  total (101/73)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

Only 100 of 175 total aligned reads displayed.

CCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACA  >  NC_002947/3985861‑3986158
                                                                                                                                                    |                                                                                                                                                     
ccACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAAt                                                                                                                                                      >  2:3286036/1‑150 (MQ=255)
 cACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATg                                                                                                                                                     >  1:4346392/1‑150 (MQ=255)
  aCGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATgg                                                                                                                                                    <  1:3258563/150‑1 (MQ=255)
    gACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCa                                                                                                                                                  >  2:2633215/1‑150 (MQ=255)
     aCCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCAt                                                                                                                                                 <  1:3894236/150‑1 (MQ=255)
      ccTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATg                                                                                                                                                >  2:3318958/1‑150 (MQ=255)
        tGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGcc                                                                                                                                              >  2:2959254/1‑150 (MQ=255)
         ggAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGccc                                                                                                                                             <  1:527743/150‑1 (MQ=255)
          gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGcccc                                                                                                                                            <  1:857958/150‑1 (MQ=255)
             aaGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGtt                                                                                                                                         <  2:4347369/150‑1 (MQ=255)
                gTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGttttt                                                                                                                                      >  2:2965523/1‑150 (MQ=255)
                    gCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAt                                                                                                                                  <  2:934032/150‑1 (MQ=255)
                    gCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAt                                                                                                                                  <  2:3202873/150‑1 (MQ=255)
                     cccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGttttcataca                                                                                                                                 <  2:2530123/150‑7 (MQ=38)
                      ccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtgt                                                                                                                                <  2:218270/150‑1 (MQ=255)
                         gTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAATCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTGTTGTGtt                                                                                                                             >  2:838587/1‑150 (MQ=255)
                            aCCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTac                                                                                                                          <  2:396487/150‑1 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTacac                                                                                                                        >  2:4060283/1‑150 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTacac                                                                                                                        >  1:2644341/1‑150 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTacac                                                                                                                        >  1:2648906/1‑150 (MQ=255)
                              ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGAGGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCCCTGGCAGGCCCCGTTTTTTCATGTGTTTacac                                                                                                                        >  2:3694952/1‑150 (MQ=255)
                               cAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACg                                                                                                                       >  1:1332431/1‑150 (MQ=255)
                               cAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACg                                                                                                                       >  1:1332417/1‑150 (MQ=255)
                                      tctACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCt                                                                                                                >  1:1671141/1‑150 (MQ=255)
                                        tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                              >  1:858116/1‑150 (MQ=255)
                                        tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                              >  1:628187/1‑150 (MQ=255)
                                        tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                              >  1:628097/1‑150 (MQ=255)
                                        tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                              >  2:3783020/1‑150 (MQ=255)
                                        tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTATTTTCATGTGTTTACACGGCCCCCTCg                                                                                                              >  1:3126508/1‑150 (MQ=255)
                                         aCAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGt                                                                                                             >  2:4160903/1‑150 (MQ=255)
                                         aCAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGt                                                                                                             <  2:1350221/150‑1 (MQ=255)
                                          cAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTa                                                                                                            <  1:3048919/150‑1 (MQ=255)
                                           aaCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAg                                                                                                           >  1:427452/1‑150 (MQ=255)
                                                  tGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCggg                                                                                                    >  1:1457218/1‑150 (MQ=255)
                                                     ccGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGttt                                                                                                 <  1:1680733/150‑1 (MQ=255)
                                                        gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATGAGCTGTATGAAATCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAcc                                                                                              >  2:2689231/1‑150 (MQ=255)
                                                         tATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAccc                                                                                             <  2:4171304/150‑1 (MQ=255)
                                                          aTCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcg                                                                                            <  2:704045/150‑1 (MQ=255)
                                                           tCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                           <  2:3047298/150‑1 (MQ=255)
                                                           tCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                           >  2:3685525/1‑150 (MQ=255)
                                                           tCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                           <  1:649447/150‑1 (MQ=255)
                                                             cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                                         >  1:3695687/1‑150 (MQ=255)
                                                             cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGa                                                                                         >  1:3063700/1‑150 (MQ=255)
                                                             cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTAAACCcgctc                                                                                         <  1:4493176/150‑3 (MQ=255)
                                                             cAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGATTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTAAACCcgctc                                                                                         <  1:1501607/150‑3 (MQ=255)
                                                              agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGaa                                                                                        <  1:1093281/150‑1 (MQ=255)
                                                                agCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCGATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATa                                                                                      <  2:695168/150‑1 (MQ=255)
                                                                   ggcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                   <  2:2412474/150‑1 (MQ=255)
                                                                   ggcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                   <  2:1179229/150‑1 (MQ=255)
                                                                    gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGcc                                                                                  <  2:3438578/150‑1 (MQ=255)
                                                                        tGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCacac                                                                              >  1:1082815/1‑150 (MQ=255)
                                                                         gTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACAcc                                                                             >  2:2913648/1‑150 (MQ=255)
                                                                         gTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACAcc                                                                             >  1:2975107/1‑150 (MQ=255)
                                                                          tAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACAccc                                                                            <  1:1742617/150‑1 (MQ=255)
                                                                          tAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACAccc                                                                            >  1:3560052/1‑150 (MQ=255)
                                                                              ggCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgcca                                                                        <  2:2242166/150‑1 (MQ=255)
                                                                               gCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccag                                                                       >  1:3119825/1‑150 (MQ=255)
                                                                               gCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCgccag                                                                       >  1:2186595/1‑150 (MQ=255)
                                                                                   cGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGGGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaa                                                                   >  2:1897361/1‑150 (MQ=255)
                                                                                            ggATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTg                                                          <  2:4102113/150‑1 (MQ=255)
                                                                                             gATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTgg                                                         <  1:1945986/150‑1 (MQ=255)
                                                                                               tCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                       <  1:3149108/150‑1 (MQ=255)
                                                                                                 gTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGgaga                                                     >  2:365364/1‑150 (MQ=255)
                                                                                                  tCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGAt                                                    >  1:2698139/1‑150 (MQ=255)
                                                                                                  tCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGAt                                                    >  1:2516688/1‑150 (MQ=255)
                                                                                                  tCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGAt                                                    >  2:2588823/1‑150 (MQ=255)
                                                                                                        tCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGt                                              >  2:1902754/1‑150 (MQ=255)
                                                                                                           ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGt                                           <  1:4128815/150‑1 (MQ=255)
                                                                                                             tGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCa                                         >  1:2126950/1‑150 (MQ=255)
                                                                                                                 aCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGc                                     <  1:3285251/150‑1 (MQ=255)
                                                                                                                   gCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAg                                   >  2:3350972/1‑150 (MQ=255)
                                                                                                                     aTCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGcc                                 <  1:1308235/150‑1 (MQ=255)
                                                                                                                             tATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGc                         >  2:3741943/1‑150 (MQ=255)
                                                                                                                              aTGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCg                        >  1:3123159/1‑150 (MQ=255)
                                                                                                                               tGAAAACGGGGCCTGCCAGCAAGGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                       <  2:1561039/150‑1 (MQ=255)
                                                                                                                               tGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                       >  2:1693817/1‑150 (MQ=255)
                                                                                                                               tGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCgg                       >  1:18770/1‑150 (MQ=255)
                                                                                                                                  aaaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCa                    >  2:1338175/1‑150 (MQ=255)
                                                                                                                                   aaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCaa                   >  2:4389910/1‑150 (MQ=255)
                                                                                                                                     cGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGc                 >  2:2986159/1‑150 (MQ=255)
                                                                                                                                      ggggCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGcc                >  1:2852252/1‑150 (MQ=255)
                                                                                                                                      ggggCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGcc                >  1:2026873/1‑150 (MQ=255)
                                                                                                                                      ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGcc                >  1:974249/1‑150 (MQ=255)
                                                                                                                                       gggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc               <  2:2486035/150‑1 (MQ=255)
                                                                                                                                         gCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcgc             <  1:363614/150‑1 (MQ=255)
                                                                                                                                          cATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcgcg            >  1:3824689/1‑150 (MQ=255)
                                                                                                                                          cATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcgcg            >  1:2428021/1‑150 (MQ=255)
                                                                                                                                          cATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCcgcg            >  1:1028868/1‑150 (MQ=255)
                                                                                                                                            tGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGcc          >  1:1393808/1‑150 (MQ=255)
                                                                                                                                             gCCAGCAATGGCATGCCCCGTATTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCTTGCCcgctaat         <  2:228241/150‑5 (MQ=255)
                                                                                                                                              ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg        >  2:4002547/1‑150 (MQ=255)
                                                                                                                                              ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg        <  1:3318958/150‑1 (MQ=255)
                                                                                                                                              ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg        <  1:3066017/150‑1 (MQ=255)
                                                                                                                                              ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg        <  1:796943/150‑1 (MQ=255)
                                                                                                                                              ccAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGGAGGAGCGGGTTGACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCgg        <  2:3136313/150‑1 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGAt                                                    <  2:3288745/105‑1 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGAt                                                    >  1:3288665/1‑105 (MQ=255)
                                                                                                                                               cAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCggg       >  1:1786483/1‑150 (MQ=255)
                                                                                                                                                  cAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGa    <  1:4389910/150‑1 (MQ=255)
                                                                                                                                                    aTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACa  <  1:3294552/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGACA  >  NC_002947/3985861‑3986158

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