Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,000,269 A→C 5.2% E330A (GAG→GCG)  PP_3526 → AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,000,2690AC 5.2% 695.4 / 20.0 202E330A (GAG→GCG) PP_3526AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (87/104);  new base C (8/3);  total (95/107)
Fisher's exact test for biased strand distribution p-value = 1.19e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 205 total aligned reads displayed.

TTACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAGGTCG  >  NC_002947/4000120‑4000417
                                                                                                                                                     |                                                                                                                                                    
ttACCTGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGa                                                                                                                                                      <  2:508137/150‑1 (MQ=255)
ttACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGa                                                                                                                                                      <  1:1231420/150‑1 (MQ=255)
 tACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAg                                                                                                                                                     >  1:3847694/1‑150 (MQ=255)
 tACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAg                                                                                                                                                     <  1:2336306/150‑1 (MQ=255)
   cATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCa                                                                                                                                                   >  2:1013951/1‑150 (MQ=255)
    aTGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCaa                                                                                                                                                  >  1:1410/1‑150 (MQ=255)
     tGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAAc                                                                                                                                                 <  1:2470276/150‑1 (MQ=255)
      ggAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACg                                                                                                                                                >  2:1086029/1‑150 (MQ=255)
      ggAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACg                                                                                                                                                >  1:2222921/1‑150 (MQ=255)
        aGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTTGCTCGCGCGACGGc                                                                                                                                              <  1:635813/150‑1 (MQ=255)
           tGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTgg                                                                                                                                           <  2:680891/150‑1 (MQ=255)
            gcgcATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTggg                                                                                                                                          >  1:4143257/1‑150 (MQ=255)
             cgcATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGc                                                                                                                                         <  2:1894570/150‑1 (MQ=255)
              gcATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCg                                                                                                                                        >  1:340451/1‑150 (MQ=255)
                aTCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAg                                                                                                                                      <  1:3788964/150‑1 (MQ=255)
                 tCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGt                                                                                                                                     <  2:1227695/150‑1 (MQ=255)
                  cACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTg                                                                                                                                    <  2:1329513/150‑1 (MQ=255)
                     cAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGAcc                                                                                                                                 <  2:843330/150‑1 (MQ=255)
                         cgccgccGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAgg                                                                                                                             <  2:909735/149‑1 (MQ=255)
                               cgGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCttt                                                                                                                       <  1:1332395/150‑1 (MQ=255)
                               cgGCTGCTGCAGTAGTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGCCCGAGGCCCttt                                                                                                                       <  2:282745/150‑1 (MQ=255)
                                ggCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGCGCAACGGCTGGGCGAGTGACCGAGGCCCTTTg                                                                                                                      >  2:3211707/1‑150 (MQ=255)
                                  ctgctgCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCTTTGagc                                                                                                                   <  2:2141109/150‑1 (MQ=255)
                                   tgctgcAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagc                                                                                                                   >  1:4359878/1‑150 (MQ=255)
                                   tgctgcAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagc                                                                                                                   <  1:2537787/150‑1 (MQ=255)
                                     ctgcAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagcca                                                                                                                 <  1:2748220/150‑1 (MQ=255)
                                        cAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccagcc                                                                                                              >  2:2050302/1‑150 (MQ=255)
                                        cAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccagcc                                                                                                              >  1:623434/1‑150 (MQ=255)
                                        cAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGAGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGagccagcc                                                                                                              >  1:4303987/1‑150 (MQ=255)
                                           tATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAg                                                                                                           >  2:4514628/1‑150 (MQ=255)
                                           tATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAg                                                                                                           >  2:1028910/1‑150 (MQ=255)
                                           tATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCACACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAg                                                                                                           >  1:557743/1‑150 (MQ=255)
                                              tCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGgcg                                                                                                        >  2:4314916/1‑150 (MQ=255)
                                                cGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCt                                                                                                      >  2:3759978/1‑150 (MQ=255)
                                                   ccTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGa                                                                                                   <  2:2084529/150‑1 (MQ=255)
                                                   acTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGATTGTCTCGGTGTCGCACTACAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGa                                                                                                   <  2:2319959/149‑1 (MQ=255)
                                                     tGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGc                                                                                                 >  2:1579455/1‑150 (MQ=255)
                                                      gTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCa                                                                                                <  2:1613426/150‑1 (MQ=255)
                                                       tCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCaa                                                                                               <  2:663479/150‑1 (MQ=255)
                                                        cGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAAt                                                                                              >  2:3213712/1‑150 (MQ=255)
                                                           tGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAAtctc                                                                                           <  1:185344/150‑1 (MQ=255)
                                                            gATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCa                                                                                          <  2:1674541/150‑1 (MQ=255)
                                                            gATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCa                                                                                          <  1:4170798/150‑1 (MQ=255)
                                                             aTGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAg                                                                                         <  1:356353/150‑1 (MQ=255)
                                                              tGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAgg                                                                                        <  2:3047832/150‑1 (MQ=255)
                                                              tGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAgg                                                                                        <  1:2373122/150‑1 (MQ=255)
                                                                    tGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGGCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAg                                                                                  <  2:784923/150‑1 (MQ=255)
                                                                    tGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAg                                                                                  <  2:1235888/150‑1 (MQ=255)
                                                                          tGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGa                                                                            >  2:2751954/1‑150 (MQ=255)
                                                                          tGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGa                                                                            <  1:2178101/150‑1 (MQ=255)
                                                                                gCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGa                                                                      <  1:2648737/150‑1 (MQ=255)
                                                                                        gTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGc                                                              >  1:3167285/1‑150 (MQ=255)
                                                                                        gTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGc                                                              <  1:4259111/150‑1 (MQ=255)
                                                                                         tctcGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAg                                                                                         <  2:90113/122‑1 (MQ=255)
                                                                                         tctcGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCAGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAg                                                                                         >  1:90113/1‑122 (MQ=255)
                                                                                           tcGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAac                                                           >  1:1558391/1‑150 (MQ=255)
                                                                                           tcGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAac                                                           >  2:807741/1‑150 (MQ=255)
                                                                                            cGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAaca                                                          <  1:4207666/150‑1 (MQ=255)
                                                                                             ggTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAacac                                                         <  2:3878573/150‑1 (MQ=255)
                                                                                                gtCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGc                                                      <  2:4116713/150‑1 (MQ=255)
                                                                                                 tCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCa                                                     >  2:4022408/1‑150 (MQ=255)
                                                                                                   gCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATc                                                   >  2:4119646/1‑150 (MQ=255)
                                                                                                    cACTTCAGCAAGTGCTATGCGGCGTATTTCGGCGACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGTTGTACCGCAACACCGCATca                                                  <  2:3173548/150‑1 (MQ=255)
                                                                                                       ttCAGCAAGTGCTATGCGGCGTATTTCGGCGACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAg                                               <  2:3504503/150‑1 (MQ=255)
                                                                                                        tCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGa                                              <  1:576475/150‑1 (MQ=255)
                                                                                                        tCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGa                                              <  2:2713090/150‑1 (MQ=255)
                                                                                                         cAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGAt                                             >  1:3199262/1‑150 (MQ=255)
                                                                                                            cAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCACTCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCAt                                          >  2:2322397/1‑150 (MQ=255)
                                                                                                            cAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTGCCGCAACACCGCATCACAGATCAt                                          <  1:1084099/150‑1 (MQ=255)
                                                                                                            cAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCAt                                          <  2:1473763/150‑1 (MQ=255)
                                                                                                             aaGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATg                                         <  2:37087/150‑1 (MQ=255)
                                                                                                             aaGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATg                                         >  2:1228665/1‑150 (MQ=255)
                                                                                                                  cTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCttt                                    <  1:2064233/150‑1 (MQ=255)
                                                                                                                   tATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTtg                                   >  1:83354/1‑150 (MQ=255)
                                                                                                                    aTGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTtgt                                  <  1:1038907/150‑1 (MQ=255)
                                                                                                                     tGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTtgtg                                 <  1:2358352/150‑1 (MQ=255)
                                                                                                                      gcggcgTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTgg                                <  2:2218452/150‑1 (MQ=255)
                                                                                                                      gcggcgTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTgg                                <  2:2173939/150‑1 (MQ=255)
                                                                                                                       cggcgTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGgc                               >  2:3761850/1‑150 (MQ=255)
                                                                                                                          cgTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCa                            <  2:3692724/150‑1 (MQ=255)
                                                                                                                          cgTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCa                            >  2:2695663/1‑150 (MQ=255)
                                                                                                                              tttCGGCTACCCGCCTTCGCGCGAGCTACGGCTGGGCGAGTTACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCATGTACGAGGCTCGAGCTGGATGTACCGCAACTCCGCATCACAGATCATGGCTTTGTGGCGCAGCtt                        <  2:422449/150‑1 (MQ=255)
                                                                                                                              tttCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCtt                        >  1:368238/1‑150 (MQ=255)
                                                                                                                              tttCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCtt                        >  1:4468042/1‑150 (MQ=255)
                                                                                                                               ttCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCttt                       <  1:1967925/150‑1 (MQ=255)
                                                                                                                                  ggCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTAtt                    >  2:1222975/1‑150 (MQ=255)
                                                                                                                                    cTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTg                  <  2:2509071/150‑1 (MQ=255)
                                                                                                                                     gcccagccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGc                 >  2:1898901/6‑150 (MQ=255)
                                                                                                                                      aCCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCt                <  2:1993664/150‑1 (MQ=255)
                                                                                                                                       cccGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGgcg                              >  2:2414233/1‑135 (MQ=255)
                                                                                                                                       cccGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtc               <  2:340451/150‑1 (MQ=255)
                                                                                                                                       ccagccgTTGCTCGCGCGACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtc               >  1:421498/4‑150 (MQ=255)
                                                                                                                                       ccagccgTTGCTCGCGCGACGGCGGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtc               >  1:750205/4‑150 (MQ=255)
                                                                                                                                          gccgTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtcatc            >  1:3958779/1‑150 (MQ=255)
                                                                                                                                          gccgTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtcatc            >  1:1476428/1‑150 (MQ=255)
                                                                                                                                           ccgTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCtcatca           <  1:3306558/150‑1 (MQ=255)
                                                                                                                                              tCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGAc        <  1:1013951/150‑1 (MQ=255)
                                                                                                                                              tCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGAc        <  1:1783901/150‑1 (MQ=255)
                                                                                                                                               cgcgcgAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACa       <  2:2598453/150‑1 (MQ=255)
                                                                                                                                                    gAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAGGTCg  >  2:4403039/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TTACATGGAGCTGCGCATCACCAAGAGCCGCCGGCTGCTGCAGTATTCCGACCTGTCGGTGATGGAAGTGGCGGTGGCCTGCGGGTTTGTCTCGGTGTCGCACTTCAGCAAGTGCTATGCGGCGTATTTCGGCTACCCGCCGTCGCGCGAGCAACGGCTGGGCGAGTGACCGAGGCCCTTTGAGCCAGCCCAGGCGCTCGAGCAATCTCAGGTACGAGGCTCGAGCTGGATGTACCGCAACACCGCATCACAGATCATGGCTTTGTGGCGCAGCTTTATTTGCTCATCAGACAGGTCG  >  NC_002947/4000120‑4000417

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: