Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,672,258 G→T 5.1% intergenic (+36/+35) nhaP → / ← PP_4975 NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,672,2580GT 5.1% 1243.4 / 37.0 371intergenic (+36/+35)nhaP/PP_4975NhaP‑type Na+(K+)/H+ antiporter/acyl‑CoA thioesterase
Reads supporting (aligned to +/- strand):  ref base G (167/184);  new base T (13/7);  total (180/191)
Fisher's exact test for biased strand distribution p-value = 1.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 376 total aligned reads displayed.

GGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACAG  >  NC_002947/5672109‑5672405
                                                                                                                                                     |                                                                                                                                                   
   gAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAc                                                                                                                                                  >  2:5941727/1‑150 (MQ=255)
    aGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAcc                                                                                                                                                 <  2:5626866/150‑1 (MQ=255)
     ggAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAccc                                                                                                                                                >  1:1753698/1‑150 (MQ=255)
     ggAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAAccc                                                                                                                                                >  1:3265961/1‑150 (MQ=255)
         cGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGa                                                                                                                                            <  1:2467653/150‑1 (MQ=255)
          gggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaa                                                                                                                                           <  1:4847379/150‑1 (MQ=255)
          gggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGaa                                                                                                                                           <  1:8008593/150‑1 (MQ=255)
          gggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGaa                                                                                                                                           <  2:3698409/150‑1 (MQ=255)
          gggACCTGCTACTGTCGATCACCTACATCGTGGTGTTGTCGTCGAGTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTGTGCGGCCCATTCGCGGGTGAACCCGCGaa                                                                                                                                           <  2:3516127/150‑1 (MQ=255)
           ggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCTCTTCGCGGGTTCACCCGCGaat                                                                                                                                          <  2:1317663/150‑2 (MQ=255)
           ggACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTTCACCCGCGaat                                                                                                                                          <  1:5834555/150‑2 (MQ=255)
            gACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAaga                                                                                                                                         >  1:6082197/1‑150 (MQ=255)
            gACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAAAGCCGGGGCCGCTTTGCGGCCCAATCGCGGGTGAACCCGCGAaga                                                                                                                                         >  1:5216555/1‑150 (MQ=255)
             aCCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAagag                                                                                                                                        <  2:5981181/150‑1 (MQ=255)
             aCCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAagag                                                                                                                                        <  2:6583952/150‑1 (MQ=255)
             aCCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAagag                                                                                                                                        <  1:4695421/150‑1 (MQ=255)
             aCCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAagag                                                                                                                                        <  2:179602/150‑1 (MQ=255)
             aCCTGCTGCTGTCGAACACCTACATCGTGGTGTTGACGTCGATACTGGTGCAGGGGTTGAGCAACGCGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAagag                                                                                                                                        <  2:2264504/150‑1 (MQ=39)
              ccTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCTCTTCGCGGGTTCACCCGCGaatggg                                                                                                                                       <  2:202383/150‑5 (MQ=255)
                tgctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGATGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCAATCGCGGGTGAACCCGCGAAGAGGcc                                                                                                                                     >  1:4571720/1‑150 (MQ=255)
                 gctgctgTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCCGGGGTTGAGCATTGGGGGGGTGGTGCGAAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCg                                                                                                                                    <  1:3788699/150‑1 (MQ=255)
                      ggTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCATTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGcaca                                                                                                                               <  2:5120811/149‑1 (MQ=255)
                       gTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAg                                                                                                                              <  2:3345457/150‑1 (MQ=255)
                       gTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAg                                                                                                                              <  2:6951858/150‑1 (MQ=255)
                          gATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGGGAACCCGCGAAGAGGCCGGCACAGGcc                                                                                                                           >  1:5477864/1‑149 (MQ=255)
                           aTCACCTACATTGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTTAGAAAGCCGGGGCCGCTTTTCGGCCCATTACCGGGGGAACCCGCGAAGAGGCCGGCCCAGGCaa                                                                                                                          >  1:7607665/1‑150 (MQ=255)
                           aTCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCaa                                                                                                                          <  2:6082197/150‑1 (MQ=255)
                                cTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAgg                                                                                                                     >  1:8598567/1‑150 (MQ=255)
                                 tACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAggg                                                                                                                    >  2:8439430/1‑150 (MQ=255)
                                  aCATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                   <  1:8791845/150‑1 (MQ=255)
                                  aCATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                   >  1:7215870/1‑150 (MQ=255)
                                  aCATCGTGGTGTGGTCGTCGATTCTGGTGCAGGGGTTGAGCATGGGGCGGGTGGGGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGt                                                                                                                   <  2:5309229/150‑1 (MQ=255)
                                    aTCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCt                                                                                                                 <  1:592707/150‑1 (MQ=255)
                                             ttgtCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTg                                                                                                        <  1:3629818/150‑1 (MQ=255)
                                                  gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGaa                                                                                                   >  2:7292120/1‑150 (MQ=255)
                                                  gtcgATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGACGGGGCCGCTTTTCGGCCCAATCGCGGGTGAACCCGCCAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGaa                                                                                                   >  1:8657623/1‑150 (MQ=255)
                                                     gATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATcc                                                                                                >  2:6713465/1‑150 (MQ=255)
                                                        tCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCggg                                                                                             >  2:6534226/1‑150 (MQ=255)
                                                        tCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCggg                                                                                             >  1:993662/1‑150 (MQ=255)
                                                         cTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGa                                                                                            >  2:6267314/1‑150 (MQ=255)
                                                          tGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGaa                                                                                           >  2:4243363/1‑150 (MQ=255)
                                                          tGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGaa                                                                                           >  1:5752493/1‑150 (MQ=255)
                                                           ggTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAAc                                                                                          <  2:3010150/150‑1 (MQ=255)
                                                                 ggggTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTcc                                                                                    >  2:3259239/1‑150 (MQ=255)
                                                                  gggTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCt                                                                                   <  1:7805290/150‑1 (MQ=255)
                                                                    gTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGc                                                                                 >  2:889013/1‑150 (MQ=255)
                                                                      tGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCAc                                                                               <  2:5688797/150‑1 (MQ=255)
                                                                          cATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTAc                                                                           <  1:405136/150‑1 (MQ=255)
                                                                           aTTGGGCGGGTGGTGCGCCAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCGTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACt                                                                          >  2:4594024/1‑150 (MQ=255)
                                                                             tGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTg                                                                        >  2:2634248/1‑150 (MQ=255)
                                                                                  gggTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGAtctc                                                                   <  2:6364193/150‑1 (MQ=255)
                                                                                   ggtggtGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                  <  2:4844945/150‑1 (MQ=255)
                                                                                   ggtggtGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCg                                                                  <  2:776361/150‑1 (MQ=255)
                                                                                    cgcgggagcAAGTTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCTCGGGTGAACCCGCGAAGAGGCCGGCAAAGGCTACCAGGGTCTATTCGACTGCCGAATCCGGTAACTGGTCCTGCACGTACTTGATCTCgg                                                                 <  1:1823597/143‑1 (MQ=255)
                                                                                       gtgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTgc                                                              <  2:217704/150‑1 (MQ=255)
                                                                                        tgCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCTGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTgcgc                                                            <  1:2326064/150‑1 (MQ=255)
                                                                                          cgcAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGcc                                                           >  1:7228291/1‑150 (MQ=255)
                                                                                            cAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCAATAGCGGGTGAACCCGCGAAGAGGCCGGCCCAGGCAACCAGGGGCTATTCGACTGCCCAATACGGGCACAGGTCCGGCACGTAATTTTTGTCGGTGCGcccc                                                         >  1:437773/1‑149 (MQ=255)
                                                                                             aaGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCAt                                                        <  1:8288817/150‑1 (MQ=255)
                                                                                              aGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTTCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATg                                                       <  1:3476985/150‑1 (MQ=255)
                                                                                               ggTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATgg                                                      >  2:7494022/1‑150 (MQ=255)
                                                                                                gTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATggg                                                     >  1:6526581/1‑150 (MQ=255)
                                                                                                 tCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGc                                                    <  1:138156/150‑1 (MQ=255)
                                                                                                      gcCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCg                                               <  1:2144716/150‑1 (MQ=255)
                                                                                                       cccAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCAATCGCGGGGGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCAGGGAACAGGGCCAGCACGTAATTGATCGCGGTGCGCCCATGGGCTGCCgg                                              >  1:8437583/1‑150 (MQ=255)
                                                                                                         cAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCa                                            >  2:7600480/1‑150 (MQ=255)
                                                                                                           gCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAg                                          <  2:2822727/150‑1 (MQ=255)
                                                                                                           gCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAg                                          <  2:793498/150‑1 (MQ=255)
                                                                                                             cTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCCACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGcc                                        >  2:4682211/1‑150 (MQ=255)
                                                                                                             cTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGcc                                        <  1:4998213/150‑1 (MQ=255)
                                                                                                              ttGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCa                                       >  2:6882498/1‑150 (MQ=255)
                                                                                                              ttGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCAGGGGCTGCCGGCCAGccc                                       >  2:1221181/1‑149 (MQ=255)
                                                                                                                   aaTGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTc                                  >  2:3447510/1‑150 (MQ=255)
                                                                                                                    gcgaatgggCCGCAAAGCGGCCCCTTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCCTCcgcg                                 >  1:811649/7‑150 (MQ=255)
                                                                                                                     tGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGc                                <  1:4155258/150‑1 (MQ=255)
                                                                                                                      gaatgggCCGCAAAGCGGCCCCTTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCCTGGGCTGCCGGCAAGCCATCCTCGcc                               >  2:8212797/5‑150 (MQ=255)
                                                                                                                      gaatgggCCGCAAAGCGGCCCCTTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGcc                               >  1:1158388/5‑150 (MQ=255)
                                                                                                                      gaatgggCCGCAAAGCGGCCCATTCGCGGGTGAACCCGCGAAGTGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGcc                               >  1:710499/5‑150 (MQ=255)
                                                                                                                        cGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCa                             <  1:5339639/150‑1 (MQ=255)
                                                                                                                        cGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCa                             <  1:5942021/150‑1 (MQ=255)
                                                                                                                         ggggCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAg                            <  2:4131060/150‑1 (MQ=255)
                                                                                                                              cGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGa                       >  2:1697419/1‑150 (MQ=255)
                                                                                                                               gCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc                      >  1:1243391/1‑150 (MQ=255)
                                                                                                                               gCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc                      >  1:5652567/1‑150 (MQ=255)
                                                                                                                               ctgtgccGGCCTCTTCGCGGGTTCACCCGCGAATGGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGAc                      >  1:6907772/7‑150 (MQ=38)
                                                                                                                                  ttGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTTGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGaa                   <  2:7051061/150‑1 (MQ=255)
                                                                                                                                  ttGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGaa                   <  2:2507613/150‑1 (MQ=255)
                                                                                                                                     cGGCCTCTTCGCGGGTTCACCCGCGAATGGGCCGGCACAGGCAACCAGGGGCGATTCGACGGCCGAATCCGGGAACGGGTCCTGCACGTACTTGATCTCGGTGCGCCCCTGGGCTGCCGGCAAGCCCTCCTCGCCCCGGGGGACGAacac                >  2:1591454/1‑150 (MQ=255)
                                                                                                                                       gCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCa              >  1:1699940/1‑150 (MQ=255)
                                                                                                                                          tcttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCt           >  2:551849/3‑150 (MQ=255)
                                                                                                                                           aTTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCtt          >  1:509825/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTTCACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTg         >  1:2920375/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTg         >  2:38184/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGAGCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGGTGACGAACACCAGCTTg         >  2:8109221/1‑150 (MQ=255)
                                                                                                                                             tCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGt        <  2:5605353/150‑1 (MQ=255)
                                                                                                                                             tCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGt        <  2:8598567/150‑1 (MQ=255)
                                                                                                                                             tCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGt        >  1:2672827/1‑150 (MQ=255)
                                                                                                                                              cgcgGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTc       <  1:6267314/150‑1 (MQ=255)
                                                                                                                                                  ggTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACa   <  1:6562951/150‑1 (MQ=255)
                                                                                                                                                   gTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACaa  >  2:1955918/1‑149 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GGCGAGGAACGGGACCTGCTGCTGTCGATCACCTACATCGTGGTGTTGTCGTCGATTCTGGTGCAGGGGTTGAGCATTGGGCGGGTGGTGCGCAAGGTCAGCGCCCAGCCTTGAGAATGCCGGGGCCGCTTTGCGGCCCATTCGCGGGTGAACCCGCGAAGAGGCCGGCACAGGCAACCAGGGTCTATTCGACTGCCGAATCCGGGAACTGGTCCTGCACGTACTTGATCTCGGTGCGCCCATGGGCTGCCGGCAAGCCATCCTCGCCCAGGTTGACGAACACCATCTTGTCGACAG  >  NC_002947/5672109‑5672405

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: