Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 384,794 C→G 5.4% A423A (GCC→GCG PP_0320 → methyl‑accepting chemotaxis transducer

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947384,7940CG 5.4% 945.3 / 30.3 280A423A (GCC→GCGPP_0320methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (138/127);  new base G (12/3);  total (150/130)
Fisher's exact test for biased strand distribution p-value = 5.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 286 total aligned reads displayed.

TCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGTCG  >  NC_002947/384645‑384942
                                                                                                                                                     |                                                                                                                                                    
tCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGcc                                                                                                                                                      >  1:5425067/1‑150 (MQ=255)
tCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGcc                                                                                                                                                      >  2:1490876/1‑150 (MQ=255)
  cGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCgg                                                                                                                                                    >  2:1392209/1‑150 (MQ=255)
  cGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCgg                                                                                                                                                    >  2:8128765/1‑150 (MQ=255)
      cGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGAt                                                                                                                                                <  1:4184097/150‑1 (MQ=255)
         cTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCg                                                                                                                                             >  2:3815246/1‑150 (MQ=255)
              aCGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCa                                                                                                                                        <  2:2645734/150‑1 (MQ=255)
                 ggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGc                                                                                                                                     <  1:693307/150‑1 (MQ=255)
                 ggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGc                                                                                                                                     >  1:5448819/1‑150 (MQ=255)
                 ggCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGc                                                                                                                                     <  2:7202945/150‑1 (MQ=255)
                    cGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTc                                                                                                                                  <  1:4407813/150‑1 (MQ=255)
                    cGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTc                                                                                                                                  >  1:8170218/1‑150 (MQ=255)
                       ggAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCa                                                                                                                               <  1:6059883/150‑1 (MQ=255)
                           cAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCg                                                                                                                           <  1:5379567/150‑1 (MQ=255)
                           cAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCg                                                                                                                           >  2:5956194/1‑150 (MQ=255)
                           cAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCg                                                                                                                           >  2:291705/1‑150 (MQ=255)
                             ggTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGc                                                                                                                         >  1:3654971/1‑150 (MQ=255)
                             ggTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGc                                                                                                                         >  2:7735697/1‑150 (MQ=255)
                              gTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCa                                                                                                                        <  2:3364607/150‑1 (MQ=255)
                                 aCCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacga                                                                                                                     <  2:7560679/150‑1 (MQ=255)
                                 aCCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacga                                                                                                                     >  2:637624/1‑150 (MQ=255)
                                  cccAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgag                                                                                                                    >  2:4427667/1‑150 (MQ=255)
                                   ccAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgaga                                                                                                                   >  2:8165929/1‑150 (MQ=255)
                                    cAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCacgagac                                                                                                                  >  1:8567090/1‑150 (MQ=255)
                                        tggtggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAgg                                                                                                              >  1:595389/1‑150 (MQ=255)
                                         ggtggACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGc                                                                                                             <  1:3700657/150‑1 (MQ=255)
                                               caacaCCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCa                                                                                                       <  2:3503484/150‑1 (MQ=255)
                                                aacaCCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAg                                                                                                      <  2:126454/150‑1 (MQ=255)
                                                 acacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAgg                                                                                                     >  1:5191740/1‑150 (MQ=255)
                                                  cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGc                                                                                                    <  2:6984915/150‑1 (MQ=255)
                                                  cacCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGc                                                                                                    >  2:3256231/1‑150 (MQ=255)
                                                   acCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGATATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGcc                                                                                                   >  2:6210147/1‑150 (MQ=255)
                                                   acCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGcc                                                                                                   >  2:2031208/1‑150 (MQ=255)
                                                    ccgccgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGccc                                                                                                  <  1:8723828/150‑1 (MQ=255)
                                                     cgccgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCcg                                                                                                 >  2:2567559/1‑150 (MQ=255)
                                                       ccgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCcgcg                                                                                               >  2:5157678/1‑150 (MQ=255)
                                                        cgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGa                                                                                              <  1:6267901/150‑1 (MQ=255)
                                                        cgAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGa                                                                                              <  1:6414671/150‑1 (MQ=255)
                                                         gAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAg                                                                                             <  2:1299279/150‑1 (MQ=255)
                                                           gcgcTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggt                                                                                           <  1:7069665/150‑1 (MQ=255)
                                                           gcacTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggt                                                                                           >  2:8848340/4‑150 (MQ=255)
                                                            cgcTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtg                                                                                          <  1:57262/150‑1 (MQ=255)
                                                              cTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAggtggt                                                                                        >  1:6192589/1‑150 (MQ=255)
                                                               tCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGgtggtg                                                                                       >  1:9088479/1‑150 (MQ=255)
                                                               tCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGgtggtg                                                                                       >  1:7651428/1‑150 (MQ=255)
                                                                 cGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCGGTACACGAGATGGGCCTGACCGTGCAGGACAGCGCCCGCACTGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCa                                                                                     >  2:7414852/1‑150 (MQ=255)
                                                                  gggCGCCAGCAGGAGATGACCGTAATTGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCTCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCTGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCaa                                                                                    <  2:9030215/150‑1 (MQ=255)
                                                                     cgcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGa                                                                                 <  1:2231503/150‑1 (MQ=255)
                                                                     cgcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGa                                                                                 <  1:2926354/150‑1 (MQ=255)
                                                                      gcCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGAt                                                                                <  2:1157049/150‑1 (MQ=255)
                                                                              gagaTGACCGCACTGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCACTGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGt                                                                        <  2:8734060/150‑1 (MQ=255)
                                                                              gagaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGt                                                                        >  2:5719680/1‑150 (MQ=255)
                                                                               agaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTg                                                                       <  1:4452984/150‑1 (MQ=255)
                                                                                gaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTgg                                                                      >  1:2084396/1‑150 (MQ=255)
                                                                                gaTGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTgg                                                                      <  1:8001124/150‑1 (MQ=255)
                                                                                   gACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCAt                                                                   >  2:502643/1‑150 (MQ=255)
                                                                                   gACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCAt                                                                   >  2:6960051/1‑150 (MQ=255)
                                                                                   gACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCAt                                                                   >  1:1202381/1‑150 (MQ=255)
                                                                                     ccGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCg                                                                 >  1:2739641/1‑150 (MQ=255)
                                                                                      cGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGa                                                                >  1:3191339/1‑150 (MQ=255)
                                                                                      cGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGa                                                                <  1:8862470/150‑1 (MQ=255)
                                                                                       gAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAg                                                               >  1:7470762/1‑150 (MQ=255)
                                                                                       gAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAg                                                               <  2:7283388/150‑1 (MQ=255)
                                                                                       gAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAg                                                               >  2:5139687/1‑150 (MQ=255)
                                                                                             gTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATg                                                         >  2:5254067/1‑150 (MQ=255)
                                                                                                 ccACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCgg                                                     <  1:1135264/150‑1 (MQ=255)
                                                                                                  cACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCggg                                                    <  1:4800953/150‑1 (MQ=255)
                                                                                                   aCTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGc                                                   <  1:6675164/150‑1 (MQ=255)
                                                                                                             acGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCa                                         >  2:2061499/1‑150 (MQ=255)
                                                                                                             acGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCa                                         >  1:5611613/1‑150 (MQ=255)
                                                                                                             acGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCa                                         >  1:1482131/1‑150 (MQ=255)
                                                                                                               gagaTGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAg                                       <  2:5425067/150‑1 (MQ=255)
                                                                                                               gagaTGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAg                                       <  1:8007837/150‑1 (MQ=255)
                                                                                                                 gaTGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAggc                                     >  1:4566776/1‑150 (MQ=255)
                                                                                                                   tGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAggcgg                                   >  1:3972468/1‑150 (MQ=255)
                                                                                                                      gCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCg                                <  1:1771373/150‑1 (MQ=255)
                                                                                                                         tGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATg                             >  1:7010043/1‑150 (MQ=255)
                                                                                                                          gACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGc                            <  1:2337288/150‑1 (MQ=255)
                                                                                                                          gACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGc                            >  2:3142516/1‑150 (MQ=255)
                                                                                                                            ccGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCgg                          <  1:5081321/150‑1 (MQ=255)
                                                                                                                            ccGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCgg                          <  2:8416697/150‑1 (MQ=255)
                                                                                                                               tGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTca                       <  2:4587717/150‑1 (MQ=255)
                                                                                                                                gCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcag                      >  2:676502/1‑150 (MQ=255)
                                                                                                                                  aGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagcc                    >  2:895212/1‑150 (MQ=255)
                                                                                                                                    gACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagccag                  >  1:4640831/1‑150 (MQ=255)
                                                                                                                                    cgcATCGCCGGCATTGCGGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagccag                  >  2:5782087/3‑150 (MQ=255)
                                                                                                                                     gcATCGCCGGCATTGCGGGCGATGTGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTcagccagc                 >  1:7503628/2‑150 (MQ=255)
                                                                                                                                      cATCGCCGGCATTGCGGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCt                >  1:8775949/1‑150 (MQ=255)
                                                                                                                                      cATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCt                >  2:4243744/1‑150 (MQ=255)
                                                                                                                                        tCGCCGGCATTGCGGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTgg              >  2:5901048/1‑150 (MQ=255)
                                                                                                                                        tCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTgg              >  2:2271064/1‑150 (MQ=255)
                                                                                                                                              gCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCcgacg        >  2:1629496/1‑150 (MQ=255)
                                                                                                                                              gCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCcgacg        >  2:3039191/1‑150 (MQ=255)
                                                                                                                                               cAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCcgacga       >  2:1767600/1‑150 (MQ=255)
                                                                                                                                                aaggCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGaa      <  2:4788169/147‑1 (MQ=255)
                                                                                                                                                aaTGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGGTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGaa      <  2:3988028/150‑1 (MQ=255)
                                                                                                                                                  tGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGt    >  2:6260368/1‑150 (MQ=255)
                                                                                                                                                  tGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGt    <  2:5461768/150‑1 (MQ=255)
                                                                                                                                                    ccGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGTTGGCCGACGAAGTCg  >  2:8222706/1‑150 (MQ=255)
                                                                                                                                                    ccGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGTCg  >  2:5626287/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
TCCGAGCGGCTGCGACGGGCCGTGGAGCAGGTAACCCAGGTGGTGGACAACACCGCCGAGCGCTCCGGGCGCCAGCAGGAGATGACCGAAATGGTCGCCACTGCTGTACACGAGATGGGCCTGACCGTGCAGGACATCGCCCGCAATGCCGGCGATGCGGCCCAGGCGTCGCAGTCGGCACGAGACGAGGCGTTGCAGGCCCGCGAGGTGGTGCAACGATCCATTCGTGGCATCGAGGGCATGTCGGGCGACATCGGCAAGGCGGCCGATGCGGTCAGCCAGCTGGCCGACGAAGTCG  >  NC_002947/384645‑384942

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: