Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,726,630 A→C 5.6% intergenic (‑40/+24) PP_4178 ← / ← htpG dienelactone hydrolase family protein/chaperone protein HtpG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,726,6300AC 5.6% 1044.0 / 33.9 305intergenic (‑40/+24)PP_4178/htpGdienelactone hydrolase family protein/chaperone protein HtpG
Reads supporting (aligned to +/- strand):  ref base A (156/131);  new base C (8/9);  total (164/141)
Fisher's exact test for biased strand distribution p-value = 6.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

Only 100 of 306 total aligned reads displayed.

CCATCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAACGGT  >  NC_002947/4726482‑4726778
                                                                                                                                                    |                                                                                                                                                    
  aTCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGc                                                                                                                                                   <  2:6752100/150‑1 (MQ=255)
  aTCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACCGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGc                                                                                                                                                   >  2:1578408/1‑150 (MQ=255)
    cACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCgg                                                                                                                                                 >  1:2234319/1‑150 (MQ=255)
        aGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCtt                                                                                                                                             <  2:4466662/150‑1 (MQ=255)
        aGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCtt                                                                                                                                             <  1:3847092/150‑1 (MQ=255)
          cccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCTTCGGCGGGCtttt                                                                                                                                           <  1:808902/150‑2 (MQ=255)
          cccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTg                                                                                                                                           >  1:8282127/1‑150 (MQ=255)
           ccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGc                                                                                                                                          >  1:5097432/1‑150 (MQ=255)
           ccTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGc                                                                                                                                          <  2:6428787/150‑1 (MQ=255)
            cTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCt                                                                                                                                         <  2:7632014/150‑1 (MQ=255)
              gggCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCtgt                                                                                                                                       <  1:3639528/150‑1 (MQ=255)
                 cggcCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCa                                                                                                                                    >  2:2965624/1‑150 (MQ=255)
                 ccgcCCAGGGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCa                                                                                                                                    >  1:7714750/3‑150 (MQ=255)
                   gcCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATc                                                                                                                                  >  1:1791555/1‑150 (MQ=255)
                     ccAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCaa                                                                                                                                >  1:3430676/1‑150 (MQ=255)
                      cAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAAc                                                                                                                               <  1:5708920/150‑1 (MQ=255)
                          gcACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCaa                                                                                                                           <  1:2747563/150‑1 (MQ=255)
                          gcACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCaa                                                                                                                           <  2:5577388/150‑1 (MQ=255)
                           cACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAACACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAg                                                                                                                          <  2:1056292/150‑1 (MQ=255)
                            aCCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGc                                                                                                                         <  2:4427564/150‑1 (MQ=255)
                            aCCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGc                                                                                                                         <  2:6985403/150‑1 (MQ=255)
                              ccGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAg                                                                                                                       >  2:8545215/1‑150 (MQ=255)
                                   tCATAGACCAAACGGCTTACGTAGGTCTTGCCAGACAGGTGATAAACCAGCCATACGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGc                                                                                                                  >  2:3732802/1‑150 (MQ=255)
                                     aTAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                                >  2:7804648/1‑150 (MQ=255)
                                     aTAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTc                                                                                                                >  1:988303/1‑150 (MQ=255)
                                      tAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCg                                                                                                               >  2:3459834/1‑150 (MQ=255)
                                       aGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGa                                                                                                              >  1:8732563/1‑150 (MQ=255)
                                        gACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGAc                                                                                                             <  2:5505472/150‑1 (MQ=255)
                                        gACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGAc                                                                                                             >  2:3098794/1‑150 (MQ=255)
                                          ccAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCa                                                                                                           <  2:4988379/150‑1 (MQ=255)
                                          ccAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCa                                                                                                           <  2:204954/150‑1 (MQ=255)
                                           cAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAg                                                                                                          >  2:3852908/1‑150 (MQ=255)
                                           cAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAg                                                                                                          >  2:835028/1‑150 (MQ=255)
                                            aaaCGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGt                                                                                                         <  2:2594961/150‑1 (MQ=255)
                                               cGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGc                                                                                                      <  2:6878331/150‑1 (MQ=255)
                                                ggCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCt                                                                                                     <  2:4035555/150‑1 (MQ=255)
                                                ggCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCt                                                                                                     >  2:250593/1‑150 (MQ=255)
                                                 gCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCtt                                                                                                    <  1:7054246/150‑1 (MQ=255)
                                                   tttCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGt                                                                                                  >  1:2975756/1‑150 (MQ=255)
                                                       gTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAg                                                                                              >  1:6619526/1‑150 (MQ=255)
                                                       gTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAg                                                                                              <  1:6848007/150‑1 (MQ=255)
                                                       gTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGAACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCTACCAGTAGCTTGTTCAg                                                                                              >  1:8603546/1‑150 (MQ=255)
                                                        tAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAgg                                                                                             >  1:3339903/1‑150 (MQ=255)
                                                        tAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAgg                                                                                             >  1:4558742/1‑150 (MQ=255)
                                                         aGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                            <  1:234755/150‑1 (MQ=255)
                                                         aGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTACATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGc                                                                                            <  1:6249719/150‑1 (MQ=255)
                                                           gTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcga                                                                                          >  2:5345295/1‑150 (MQ=255)
                                                           gTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcga                                                                                          >  2:3968806/1‑150 (MQ=255)
                                                              ttGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcgacga                                                                                       >  1:7451220/1‑150 (MQ=255)
                                                              atGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGcgacga                                                                                       >  1:6184942/2‑150 (MQ=255)
                                                                gCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAAc                                                                                     >  2:1878267/1‑150 (MQ=255)
                                                                gCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAAc                                                                                     <  1:5253925/150‑1 (MQ=255)
                                                                        ggTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAgccg                                                                             >  1:6502029/1‑150 (MQ=255)
                                                                         gTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAgccgc                                                                            >  1:7924199/1‑150 (MQ=255)
                                                                            aTAAACCAGCTATTCGACTATTACTTTTCTCATAAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGATTGTTCAGGCGACGAACGTAAGCCGCCgg                                                                         <  1:119147/150‑1 (MQ=255)
                                                                            aTAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCCCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCgg                                                                         <  2:3772608/150‑1 (MQ=255)
                                                                                  cAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCtt                                                                   >  1:763316/1‑150 (MQ=255)
                                                                                  cAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGAAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCCACGGGTCCtt                                                                   >  1:369494/1‑150 (MQ=255)
                                                                                   cgCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGCGTCCTTc                                                                  >  2:5156215/2‑150 (MQ=255)
                                                                                   aGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTc                                                                  >  2:8215331/1‑150 (MQ=255)
                                                                                   aGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTc                                                                  <  2:6532522/150‑1 (MQ=255)
                                                                                     cGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAg                                                                <  1:1644029/150‑1 (MQ=255)
                                                                                      gATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAgg                                                               <  2:8107542/150‑1 (MQ=255)
                                                                                      gATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAgg                                                               <  1:6479626/150‑1 (MQ=255)
                                                                                        ttCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCt                                                             >  2:7462415/1‑150 (MQ=255)
                                                                                        ttCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCt                                                             <  1:5835027/150‑1 (MQ=255)
                                                                                              gaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCAcc                                                       <  2:8658762/150‑1 (MQ=255)
                                                                                              gaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCAcc                                                       >  1:3986695/1‑150 (MQ=255)
                                                                                              gaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCAcc                                                       >  2:6167666/1‑150 (MQ=255)
                                                                                              gaTTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCAcc                                                       >  2:7139200/1‑150 (MQ=255)
                                                                                              gaTTACTTTGCTCATCAAGCACTCCTGTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACCGCTCGACCAGTAGCGTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCGGCCGGCTGTCAcc                                                       >  2:8260710/1‑150 (MQ=255)
                                                                                                       gCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAg                                              >  2:6030861/1‑150 (MQ=255)
                                                                                                         tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                            >  1:3089042/1‑150 (MQ=255)
                                                                                                         tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                            >  2:4315755/1‑150 (MQ=255)
                                                                                                         tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                            >  1:6034989/1‑150 (MQ=255)
                                                                                                         tcatcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAggg                                            >  1:4115143/1‑150 (MQ=255)
                                                                                                           atcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGgcc                                          <  2:5896971/150‑1 (MQ=255)
                                                                                                           aacaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCCAAGCGGGCTTTGCTGTGCATCAACTCAAGCACACAGCACGACCAGTAGCTTGTTCAGGCGACGAACCTAAGCCCCCGGGTCCTTCAGCCTGACACCGGCCGCCAGCgcc                                          >  1:1857508/3‑150 (MQ=255)
                                                                                                            tcaAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGgccg                                         >  2:6673736/1‑150 (MQ=255)
                                                                                                             caAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGgccgc                                        <  1:2677723/150‑1 (MQ=255)
                                                                                                               aGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCt                                      >  2:8193293/1‑150 (MQ=255)
                                                                                                                 cACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGa                                    >  2:8610056/1‑150 (MQ=255)
                                                                                                                  aCTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGAt                                   <  2:172679/150‑1 (MQ=255)
                                                                                                                       tttGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAAc                              >  1:4786845/1‑150 (MQ=255)
                                                                                                                       tttGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAAc                              >  2:4736345/1‑150 (MQ=255)
                                                                                                                        ttGGCACTCGTCATAAAAAAGCCCGCCTAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACa                             <  1:1012523/150‑1 (MQ=255)
                                                                                                                        ttGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACa                             <  1:8379327/150‑1 (MQ=255)
                                                                                                                                cGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGc                     <  2:7714750/150‑1 (MQ=255)
                                                                                                                                  tCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGa                   >  2:2228655/1‑150 (MQ=255)
                                                                                                                                   cATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGAc                  >  2:2343891/1‑150 (MQ=255)
                                                                                                                                    aTAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACa                 >  1:489401/1‑150 (MQ=255)
                                                                                                                                     tAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAg                >  2:3604614/1‑150 (MQ=255)
                                                                                                                                     tAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAg                >  1:4942567/1‑150 (MQ=255)
                                                                                                                                     cagcaaaGCCCGCTTCGGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAg                >  2:7352547/5‑150 (MQ=255)
                                                                                                                                      aaaaaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGc               >  1:8153447/1‑150 (MQ=255)
                                                                                                                                       gcaaaGCCCGCTTCGGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCt              >  1:546536/3‑150 (MQ=255)
                                                                                                                                         aaaGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCg            >  1:1948521/1‑150 (MQ=255)
                                                                                                                                              ccgccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGaa       <  1:6726365/150‑1 (MQ=255)
                                                                                                                                              ccgccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGaa       >  2:5785691/1‑150 (MQ=255)
                                                                                                                                                 ccgAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAACg    >  1:246354/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CCATCACCAGCCCTGGGCGGCCCAGCGCACCCGGCTCATAGACCAAACGGCTTTCGTAGGTCTTGCCAGACAGGTGATAAACCAGCGATTCGACGATTACTTTGCTCATCAAGCACTCCTTTGGCACTCGTCATAAAAAAGCCCGCCGAAGCGGGCTTTGCTGTGCATCAACTCAAGCAGACAGCTCGACCAGTAGCTTGTTCAGGCGACGAACGTAAGCCGCCGGGTCCTTCAGGCTGTCACCGGCCGCCAGGGCCGCCTGATCGAACAGGATGTGCGACAGCTCGGCGAAACGGT  >  NC_002947/4726482‑4726778

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: