Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,148,875 T→A 7.4% intergenic (+35/‑88) PP_3647 → / → PP_3648 oxidoreductase/carboxymuconolactone decarboxylase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,148,8750TA 7.4% 946.0 / 47.4 305intergenic (+35/‑88)PP_3647/PP_3648oxidoreductase/carboxymuconolactone decarboxylase family protein
Reads supporting (aligned to +/- strand):  ref base T (125/155);  new base A (11/12);  total (137/168)
Fisher's exact test for biased strand distribution p-value = 8.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.23e-01

Only 100 of 310 total aligned reads displayed.

CGGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTA  >  NC_002947/4148726‑4149024
                                                                                                                                                     |                                                                                                                                                     
cggTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCt                                                                                                                                                       <  2:2214030/150‑1 (MQ=255)
     aCCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGccc                                                                                                                                                  <  1:8392398/150‑1 (MQ=255)
       cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgc                                                                                                                                                <  2:249846/150‑1 (MQ=255)
       cACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgc                                                                                                                                                <  2:3046002/150‑1 (MQ=255)
        aCGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgcg                                                                                                                                               <  2:7595432/150‑1 (MQ=255)
        aCGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCcgcg                                                                                                                                               <  1:2911846/150‑1 (MQ=255)
         cgcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGa                                                                                                                                              <  1:792007/150‑1 (MQ=255)
         cgcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGa                                                                                                                                              >  2:1655654/1‑150 (MQ=255)
         cgcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGa                                                                                                                                              <  1:1792154/150‑1 (MQ=255)
         cgcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGa                                                                                                                                              <  2:2387591/150‑1 (MQ=255)
          gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             >  2:7566569/1‑150 (MQ=255)
          gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             >  2:7521252/1‑150 (MQ=255)
          gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             <  2:6358206/150‑1 (MQ=255)
          gcgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGaa                                                                                                                                             <  1:3332817/150‑1 (MQ=255)
           cgATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAg                                                                                                                                            <  1:1219270/150‑1 (MQ=255)
            gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                           <  1:2515514/150‑1 (MQ=255)
            gATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGa                                                                                                                                           <  2:5114933/150‑1 (MQ=255)
             aTCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCAAGCCCGCGAAGAt                                                                                                                                          <  1:5890390/150‑1 (MQ=255)
              tCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATg                                                                                                                                         <  1:4752975/150‑1 (MQ=255)
               cATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACACGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGc                                                                                                                                        <  1:4663149/150‑1 (MQ=255)
                 ttGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGccc                                                                                                                                      >  2:7846572/1‑150 (MQ=255)
                          cgGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACACGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTGGCCCGCGAAGATGCCCCTGTTTTCa                                                                                                                             <  1:5175421/150‑1 (MQ=255)
                           ggTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGTTACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                            <  2:5424378/150‑1 (MQ=255)
                           ggTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                            <  2:5728401/150‑1 (MQ=255)
                           ggTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                            <  2:4420780/150‑1 (MQ=255)
                            gttaagaaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGTTGTTGTACCACATGGGCAAGCTCGGCGCCCTGCCGCAACTTGGCTTATACATCTTGAGTCTTTGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                           <  1:3856809/147‑1 (MQ=255)
                                  aaCCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTcgc                                                                                                                     <  2:1420792/150‑1 (MQ=255)
                                          gtTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCg                                                                                                             >  1:8316315/1‑150 (MQ=255)
                                            tGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTTATACATCTTGAGTCGATGTTGGAATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtc                                                                                                           <  1:260371/150‑1 (MQ=255)
                                             gCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                          <  2:4955866/150‑1 (MQ=255)
                                             gCCTGAGCAGCCGATGACCCCGAGGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTGGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtct                                                                                                          <  2:1262134/150‑1 (MQ=255)
                                              ccTGCGCAGCCGATGACCCCGAGGGGGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtctc                                                                                                         <  1:5442093/150‑1 (MQ=255)
                                              ccTGAGCCGCCGATGACCCCGATGGTGTACCACACGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGtctc                                                                                                         <  1:4966863/150‑1 (MQ=255)
                                                tGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCAc                                                                                                       >  2:5820184/1‑150 (MQ=255)
                                                          aTGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGc                                                                                             <  2:5884979/150‑1 (MQ=255)
                                                             aCCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGc                                                                                          >  2:6590257/1‑150 (MQ=255)
                                                             aCCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGc                                                                                          >  2:515819/1‑150 (MQ=255)
                                                              ccccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGcc                                                                                         >  1:8024510/1‑150 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAgg                                                                                                                           >  1:4332413/1‑115 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                        >  2:3962452/1‑150 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                        >  2:7723877/1‑150 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                        >  2:1903828/1‑150 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                        >  2:2059870/1‑150 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                        >  1:5199619/1‑150 (MQ=255)
                                                               cccGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTGGCCCGCGAAGATGCCCCTGTGTTCAGGGGGCGCGCTAATCGTCTCACCCTCGGTTGCGGccc                                                                                        >  2:8530852/1‑150 (MQ=255)
                                                                    tGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCtttt                                                                                   <  2:810287/150‑1 (MQ=255)
                                                                     ggTGTACCACACGGGCACGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGCGTCGAGGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCttttt                                                                                  <  1:4694465/150‑1 (MQ=255)
                                                                      gtgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCgg                                                                                                 <  2:7203612/134‑1 (MQ=255)
                                                                       tgtACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCtttt                                                                                   <  1:5546833/147‑1 (MQ=255)
                                                                            cacaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCg                                                                           <  1:2277918/150‑1 (MQ=255)
                                                                             acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCgg                                                                          >  1:2201988/1‑150 (MQ=255)
                                                                             acaAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCgg                                                                          >  1:7274586/1‑150 (MQ=255)
                                                                              caAGGGCAAGCTTGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAg                                                                                                                            <  2:4650282/99‑1 (MQ=255)
                                                                              caAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGa                                                                         <  2:6908341/150‑1 (MQ=255)
                                                                                   gCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACt                                                                    <  1:3621979/150‑1 (MQ=255)
                                                                                    cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACTTCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTg                                                                   <  2:356790/150‑1 (MQ=255)
                                                                                    cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTg                                                                   <  1:7329620/150‑1 (MQ=255)
                                                                                    cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTg                                                                   <  1:3692985/150‑1 (MQ=255)
                                                                                    cAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTGTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTg                                                                   >  2:6107478/1‑150 (MQ=255)
                                                                                        cTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTGGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTCCTGccca                                                               >  2:3451252/1‑150 (MQ=255)
                                                                                         tCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACAt                                                              <  1:1955364/150‑1 (MQ=255)
                                                                                           ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGa                                                            <  1:970003/150‑1 (MQ=255)
                                                                                           ggCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCACGGTTGCGGCCCCATTTTTATGCCGGAGTACTGCACATGa                                                            >  1:4015139/1‑150 (MQ=255)
                                                                                            gcgcCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAg                                                           <  2:7453052/150‑1 (MQ=255)
                                                                                                 cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATg                                                      >  2:3084583/1‑150 (MQ=255)
                                                                                                 cTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATg                                                      >  2:3242301/1‑150 (MQ=255)
                                                                                                  tGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGa                                                     >  2:6176142/1‑150 (MQ=255)
                                                                                                   gccgccACTTGGCTGATACATCTAGAGTCGAAGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaa                                                    <  2:3808106/150‑1 (MQ=255)
                                                                                                     cgccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGaaaa                                                  <  1:3113929/150‑1 (MQ=255)
                                                                                                       ccACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGATAAGt                                                >  2:2640052/1‑150 (MQ=255)
                                                                                                        cACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTa                                               <  2:5622935/150‑1 (MQ=255)
                                                                                                         aCTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTGTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTAc                                              >  1:7095887/1‑150 (MQ=255)
                                                                                                           ttGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGa                                            >  1:8609986/1‑150 (MQ=255)
                                                                                                           ttGGCTGATACATCTTGAGTCGATGTAGGCATCATCGCGGGCTTGCCCGCGAAGATGCCCCTCTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGa                                            >  2:1782103/1‑150 (MQ=255)
                                                                                                              gCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGaaaa                                         <  2:2933818/150‑1 (MQ=255)
                                                                                                              gCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATACCGGAGTACTGCACATGAGCAATGAAAAGTACGaaaa                                         <  1:7702145/150‑1 (MQ=255)
                                                                                                                tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAACAgg                                       >  2:8407762/1‑150 (MQ=255)
                                                                                                                tGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAgg                                       >  2:5295596/1‑150 (MQ=255)
                                                                                                                   tACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTc                                                                                                                              >  1:3472367/1‑60 (MQ=255)
                                                                                                                       tCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaa                                >  1:6755823/1‑150 (MQ=255)
                                                                                                                        cTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGaaa                               <  1:1528148/150‑1 (MQ=255)
                                                                                                                                   gTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCCCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCa                    >  2:1136177/1‑150 (MQ=255)
                                                                                                                                    ttGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAg                   <  1:7555894/150‑1 (MQ=255)
                                                                                                                                         aTCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCACTGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACt              >  2:7301883/1‑150 (MQ=255)
                                                                                                                                           cTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTgg            >  1:1678085/1‑150 (MQ=255)
                                                                                                                                           cTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTgg            <  2:1233508/150‑1 (MQ=255)
                                                                                                                                           cTTCGCGGGCAAGCCCGCGAAGATGCAACTGTTTTCAGGGGTCGCGCTAATAGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGTGCCTGAAAATCCGCACCCAGGTACTgg            >  2:479247/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg           >  1:525424/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg           >  2:6230104/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTggg           >  1:8075426/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGc          <  1:6262231/150‑1 (MQ=255)
                                                                                                                                             tCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGc          >  2:1447518/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGCAAGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGc          >  1:1315573/1‑150 (MQ=255)
                                                                                                                                               gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa        <  1:1737017/150‑1 (MQ=255)
                                                                                                                                               gcgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGa        <  1:6503209/150‑1 (MQ=255)
                                                                                                                                                cgGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAg       <  2:1755908/150‑1 (MQ=255)
                                                                                                                                                  ggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGa     >  1:1864874/1‑150 (MQ=255)
                                                                                                                                                  ggCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGa     <  2:8328094/150‑1 (MQ=255)
                                                                                                                                                    cTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTGGCGGCCCCTTGTTTATGCCGGAGTACTGCACCTGAGCCCTGCACCGGACGGAGCGGGCCTGACGCGCCGCCCCCCGGGCCTGGGCGAGGTCt   >  2:6527543/1‑150 (MQ=255)
                                                                                                                                                     ttGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTa  <  1:5820184/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGGTGACCACGCGATCATTGTCGGCGCGGTGAAGAACCTGGTGTTGCCTGAGCAGCCGATGACCCCGATGGTGTACCACAAGGGCAAGCTCGGCGCCCTGCCGCCACTTGGCTGATACATCTTGAGTCGATGTTGGCATCTTCGCGGGCTTGCCCGCGAAGATGCCCCTGTTTTCAGGGGTCGCGCTAATCGTCTCACCCTCGGTTGCGGCCCCTTTTTTATGCCGGAGTACTGCACATGAGCAATGAAAAGTACGAAAAGGGCCTGAAAATCCGCACCCAGGTACTGGGCGAGGACTA  >  NC_002947/4148726‑4149024

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