Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,345,571:1 +A 5.4% coding (249/1587 nt) ybiT ← ABC‑F family ATPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,345,5711.A 5.4% 1214.1 / 58.9 279coding (249/1587 nt)ybiTABC‑F family ATPase
Reads supporting (aligned to +/- strand):  ref base . (136/128);  new base A (3/12);  total (139/140)
Fisher's exact test for biased strand distribution p-value = 3.05e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 284 total aligned reads displayed.

ACTCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCC  >  NC_002947/3345423‑3345720
                                                                                                                                                     |                                                                                                                                                     
 cTCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTc                                                                                                                                                     <  2:2688623/150‑1 (MQ=255)
  tCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCg                                                                                                                                                    >  1:582343/1‑150 (MQ=255)
           tAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTg                                                                                                                                           >  2:4624107/1‑150 (MQ=255)
           gaGCCGTCCATGTCGGCGAATTCGGGTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTg                                                                                                                                           <  2:579152/149‑1 (MQ=255)
            aGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCCCAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGa                                                                                                                                          <  2:4827313/150‑1 (MQ=255)
            aGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGa                                                                                                                                          <  2:3575293/150‑1 (MQ=255)
            aGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGa                                                                                                                                          <  1:1773190/150‑1 (MQ=255)
               cGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGatca                                                                                                                                       <  2:1924687/150‑4 (MQ=255)
                gTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCGCCGTGTCGATCACGGTGAATTc                                                                                                                                      >  2:2024261/1‑150 (MQ=255)
                 tCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTcc                                                                                                                                     >  1:5165516/1‑150 (MQ=255)
                 tCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTcc                                                                                                                                     >  1:5108384/1‑150 (MQ=255)
                 tCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTcc                                                                                                                                     >  1:3071280/1‑150 (MQ=255)
                  ccATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGATCGACACGGTGatcatgg                                                                                                                                    <  1:4222524/150‑7 (MQ=255)
                   cATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTc                                                                                                                                   <  2:4711043/150‑1 (MQ=255)
                      ttCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTa                                                                                                                                <  1:2352150/150‑1 (MQ=255)
                       tCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAg                                                                                                                               <  1:4936814/150‑1 (MQ=255)
                        cGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAgg                                                                                                                              <  1:1681643/150‑1 (MQ=255)
                            gAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGaa                                                                                                                          >  2:5404369/1‑150 (MQ=255)
                               ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                       <  2:5060330/150‑1 (MQ=255)
                               ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                       >  2:132949/1‑150 (MQ=255)
                               ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                       >  2:2925990/1‑150 (MQ=255)
                               ttCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTg                                                                                                                       >  2:2581386/1‑150 (MQ=255)
                                 cGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTa                                                                                                                                >  1:632987/1‑139 (MQ=255)
                                 cGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGt                                                                                                                     >  2:1240851/1‑150 (MQ=255)
                                  ggTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTc                                                                                                                    >  2:1499268/1‑150 (MQ=255)
                                    aaacgaGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCt                                                                                                                  <  1:3162268/145‑1 (MQ=255)
                                         aGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGgcgc                                                                                                             <  2:3443459/150‑1 (MQ=255)
                                         aGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGgcgc                                                                                                             >  1:2143069/1‑150 (MQ=255)
                                         aGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGgcgc                                                                                                             >  1:4889383/1‑150 (MQ=255)
                                          gCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCa                                                                                                            <  1:3018575/150‑1 (MQ=255)
                                           cTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAg                                                                                                           >  1:4362132/1‑150 (MQ=255)
                                            tCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGc                                                                                                          <  1:554767/150‑1 (MQ=255)
                                              ggCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCtt                                                                                                        <  2:5015683/150‑1 (MQ=255)
                                                cGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGc                                                                                                      <  2:554608/150‑1 (MQ=255)
                                                 gACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGcc                                                                                                     >  2:4243301/1‑150 (MQ=255)
                                                    ggCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCAATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAg                                                                                                  >  2:3921417/1‑150 (MQ=255)
                                                      ccATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGc                                                                                                >  1:2607045/1‑150 (MQ=255)
                                                      ccATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGc                                                                                                <  2:3289491/150‑1 (MQ=255)
                                                       cATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCg                                                                                               >  1:136635/1‑150 (MQ=255)
                                                          gCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGt                                                                                            >  1:4045165/1‑150 (MQ=255)
                                                           ccGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTa                                                                                           <  2:3778448/150‑1 (MQ=255)
                                                              tcttctTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTc                                                                                        >  2:4494552/1‑150 (MQ=255)
                                                                ttcttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCgg                                                                                      >  2:1418381/1‑150 (MQ=255)
                                                                   ttcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                   <  1:1002588/150‑1 (MQ=255)
                                                                   ctcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                   <  1:35749/149‑1 (MQ=255)
                                                                   ctcGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTc                                                                                   <  1:36295/149‑1 (MQ=255)
                                                                     cGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGa                                                                                 >  1:4673573/1‑150 (MQ=255)
                                                                      ggTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAg                                                                                >  2:93196/1‑150 (MQ=255)
                                                                      ggTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAg                                                                                >  2:610338/1‑150 (MQ=255)
                                                                       gTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGc                                                                               >  1:3888761/1‑150 (MQ=255)
                                                                       gTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGc                                                                               >  2:4455101/1‑150 (MQ=255)
                                                                         cATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCAt                                                                             >  2:2304503/1‑150 (MQ=255)
                                                                           tttCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGa                                                                           >  2:524536/1‑150 (MQ=255)
                                                                           tttCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGa                                                                           >  1:1830136/1‑150 (MQ=255)
                                                                            gtCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGGGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGAc                                                                          <  2:5198780/149‑1 (MQ=255)
                                                                                 gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                     <  2:3901507/150‑1 (MQ=255)
                                                                                 gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                     <  2:3901995/150‑1 (MQ=255)
                                                                                 gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                     <  2:1801336/150‑1 (MQ=255)
                                                                                 gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                     <  2:4343164/150‑1 (MQ=255)
                                                                                 gCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGc                                                                     >  2:2648519/1‑150 (MQ=255)
                                                                                      gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                >  2:3035231/1‑150 (MQ=255)
                                                                                      gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                >  2:309/1‑150 (MQ=255)
                                                                                      gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                >  1:4286970/1‑150 (MQ=255)
                                                                                      gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                <  1:4479291/150‑1 (MQ=255)
                                                                                      gAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCg                                                                >  1:924794/1‑150 (MQ=255)
                                                                                        gTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGa                                                              <  1:84222/150‑1 (MQ=255)
                                                                                         tAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGAc                                                             <  2:3139069/150‑1 (MQ=255)
                                                                                           gATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACgg                                                           >  2:656181/1‑150 (MQ=255)
                                                                                                 gTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAg                                                     <  1:3790113/150‑1 (MQ=255)
                                                                                                     cGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCa                                                 >  1:1121368/1‑150 (MQ=255)
                                                                                                        tCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAccg                                              >  1:3258531/1‑150 (MQ=255)
                                                                                                         cGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAccgc                                             <  2:2438316/150‑1 (MQ=255)
                                                                                                          ggCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCATGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCAccgcc                                            >  1:5036229/1‑150 (MQ=255)
                                                                                                            cTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGa                                          <  2:2797311/150‑1 (MQ=255)
                                                                                                             tttGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAg                                         >  1:2506362/1‑150 (MQ=255)
                                                                                                               tGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGa                                       <  2:2235306/150‑1 (MQ=255)
                                                                                                               tGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGa                                       <  1:2372223/150‑1 (MQ=255)
                                                                                                                   tttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtttt                                   >  1:1207835/1‑150 (MQ=255)
                                                                                                                   tttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtttt                                   >  1:310390/1‑150 (MQ=255)
                                                                                                                   tttCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTAGTTGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGAtttt                                   >  2:353948/1‑150 (MQ=255)
                                                                                                                     tCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCa                                 <  2:5108384/150‑1 (MQ=255)
                                                                                                                      ccACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCAt                                >  2:3674299/1‑150 (MQ=255)
                                                                                                                       cacaGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATg                               >  2:3376941/1‑150 (MQ=255)
                                                                                                                        acaGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGa                              >  1:2244939/1‑150 (MQ=255)
                                                                                                                          aGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAg                            <  1:1058550/150‑1 (MQ=255)
                                                                                                                          aGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAg                            >  2:3223440/1‑150 (MQ=255)
                                                                                                                          aGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAg                            >  2:4889507/1‑150 (MQ=255)
                                                                                                                           gCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAgg                           >  2:757193/1‑150 (MQ=255)
                                                                                                                             tGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTc                         >  2:3194096/1‑150 (MQ=255)
                                                                                                                             tGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTc                         >  2:1088620/1‑150 (MQ=255)
                                                                                                                                    ggCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACAGCCGAGGATTTTCATGAAGGTCGACTTGc                  <  2:825654/150‑1 (MQ=255)
                                                                                                                                       ccATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGc               >  1:1571749/1‑150 (MQ=255)
                                                                                                                                         aTGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCaa             <  1:4095342/150‑1 (MQ=255)
                                                                                                                                           gATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAAcc           >  2:2314695/1‑150 (MQ=255)
                                                                                                                                            ttCACCGTGATCGACACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAGGGTCGACTTGCCGCAACCg          >  1:3877396/2‑150 (MQ=255)
                                                                                                                                             tCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGt         >  1:1096476/1‑150 (MQ=255)
                                                                                                                                              cACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGtt        <  2:4124092/150‑1 (MQ=255)
                                                                                                                                              cACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGtt        >  2:4069222/1‑150 (MQ=255)
                                                                                                                                               aCCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTg       <  2:1147369/150‑1 (MQ=255)
                                                                                                                                                 cGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGgc     <  1:2713066/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTCGGCGGTGTAGCCGTCCATTTCGGCGAATTCGGTTTCCAGCTCGGCGACGGCCATGCCGTCTTCTTCGGTCATTTCCGGCAGCGAGTAGATACGGTCACGCTCGGCTTTGACTTTCCACAGCTGGCCGTGGCCCATGATCACCGTGTCGATCACGGTGAATTCCTCGTAGGCGAACTGGTCCTGGCGCAGCTTGCCCAGGCGGGTATTCGGCTCGAGCATGACCTGGCCTGCGGACGGCTCCAGGTCACCGCCGAGGATTTTCATGAAGGTCGACTTGCCGCAACCGTTGGCGCC  >  NC_002947/3345423‑3345720

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: