Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,699 G→T 7.2% P41Q (CCG→CAG)  parA ← chromosome partition protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,6990GT 7.2% 400.4 / 14.1 125P41Q (CCG→CAG) parAchromosome partition protein
Reads supporting (aligned to +/- strand):  ref base G (59/57);  new base T (4/5);  total (63/62)
Fisher's exact test for biased strand distribution p-value = 7.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.53e-01

Only 100 of 128 total aligned reads displayed.

GGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTG  >  NC_002947/1556‑1848
                                                                                                                                               |                                                                                                                                                     
 gTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                <  1:625878/150‑1 (MQ=255)
 gTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                <  1:1413030/150‑1 (MQ=255)
 gTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGATCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                <  2:944517/150‑1 (MQ=255)
  ttGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGGCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAgg                                                                                                                                               <  1:1968370/150‑1 (MQ=255)
     gCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGTGGGCTTGCCCTGCGGATCGAgggcca                                                                                                                                           >  2:2105254/1‑146 (MQ=255)
       cGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAt                                                                                                                                          <  2:2270276/150‑1 (MQ=255)
            gcaACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgca                                                                                                                                     <  2:189030/150‑1 (MQ=255)
             caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCGTGGTGGCGTTGCCCTGCTGATCGACCTCGATccgcag                                                                                                                                    <  1:2198075/150‑1 (MQ=255)
             caACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATcagcag                                                                                                                                    <  2:1886198/150‑1 (MQ=255)
                cTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGac                                                                                                                                 >  1:531005/1‑150 (MQ=255)
                 tGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAACAGGaca                                                                                                                                >  1:103368/1‑150 (MQ=255)
                   gAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACacg                                                                                                                              >  2:351363/1‑150 (MQ=255)
                     aGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACacgac                                                                                                                            <  2:1620327/150‑1 (MQ=255)
                       ccTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGc                                                                                                                          <  2:1959845/150‑1 (MQ=255)
                         tCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGa                                                                                                                                                     >  1:485818/1‑121 (MQ=255)
                         tCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGa                                                                                                                                                     <  2:485818/121‑1 (MQ=255)
                           cGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTgg                                                                                                                      <  1:1629223/150‑1 (MQ=255)
                                tCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGca                                                                                                                 >  2:1407983/1‑150 (MQ=255)
                                tCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGca                                                                                                                 >  2:2003694/1‑150 (MQ=255)
                                  ggAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                <  1:2407465/117‑1 (MQ=255)
                                  ggAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAg                                                                                                                                                >  2:2407465/1‑117 (MQ=255)
                                      tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                        >  1:403438/1‑137 (MQ=255)
                                      tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCtt                                                                                                                        <  2:403438/137‑1 (MQ=255)
                                      tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTg                                                                                                                       >  1:447673/1‑138 (MQ=255)
                                      tagtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTg                                                                                                                       <  2:447673/138‑1 (MQ=255)
                                        gtGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCg                                                                                                                                                       <  1:920645/104‑1 (MQ=255)
                                            aTAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATCCGCAgg                                                                                                                                   <  2:441431/120‑1 (MQ=255)
                                             tAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcgg                                                                                                    <  1:1213415/150‑1 (MQ=255)
                                             tAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcgg                                                                                                    >  1:1299998/1‑150 (MQ=255)
                                             tAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggcgg                                                                                                    >  2:1253111/1‑150 (MQ=255)
                                                      ccAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGTAGCCAGCGAGGCGGCGAGATTGa                                                                                           <  2:1502441/150‑1 (MQ=255)
                                                      ccAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGa                                                                                           >  1:1842029/1‑150 (MQ=255)
                                                       cAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGAt                                                                                          >  1:2136943/1‑150 (MQ=255)
                                                        aaGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATa                                                                                         <  1:1306652/150‑1 (MQ=255)
                                                           tCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                      >  2:1038708/1‑150 (MQ=255)
                                                           tCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAg                                                                                      >  1:111306/1‑150 (MQ=255)
                                                             gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgcag                                                                                                                                    >  1:2276146/1‑102 (MQ=255)
                                                             gCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCTGATCGACCTCGATccgcag                                                                                                                                    <  2:2276146/102‑1 (MQ=255)
                                                              cACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtg                                                                                   >  1:472139/1‑150 (MQ=255)
                                                              cACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAgg                                                                                     >  2:120062/1‑148 (MQ=255)
                                                              cACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACAACGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGGGACAggtg                                                                                   >  1:102754/1‑150 (MQ=255)
                                                               aCTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgt                                                                                  <  2:654376/150‑2 (MQ=255)
                                                                   cccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGt                                                                              >  2:779186/1‑150 (MQ=255)
                                                                   cccGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGt                                                                              >  1:862215/1‑150 (MQ=255)
                                                                     cGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCt                                                                            >  1:68235/1‑150 (MQ=255)
                                                                           gCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTg                                                                          >  1:1922451/1‑146 (MQ=255)
                                                                            cAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtg                                                                                   >  2:1363311/1‑136 (MQ=255)
                                                                            cAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTCCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtg                                                                                   <  1:1363311/136‑1 (MQ=255)
                                                                            cAGGTCGTAGACAGAATGCTCGATCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTGGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCcac                                                                     <  1:2300499/150‑1 (MQ=255)
                                                                                   tAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCt                                                              <  2:862215/150‑1 (MQ=255)
                                                                                    aGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAggtgg                                                                                  >  1:752946/1‑129 (MQ=255)
                                                                                      aCCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTc                                                           <  2:472139/150‑1 (MQ=255)
                                                                                        cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAt                                                                                              <  2:1522785/113‑1 (MQ=255)
                                                                                        cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGAt                                                                                              >  1:1522785/1‑113 (MQ=255)
                                                                                        cGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTg                                                         <  1:2368681/150‑1 (MQ=255)
                                                                                         gAATGCTCGAGCTCTTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTgg                                                        <  1:430178/150‑1 (MQ=255)
                                                                                         gAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGtt                                                                                <  1:1700729/126‑1 (MQ=255)
                                                                                         gAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGtt                                                                                >  2:1700729/1‑126 (MQ=255)
                                                                                          aaTGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                             <  2:1985810/128‑1 (MQ=255)
                                                                                          aaTGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTc                                                                             >  1:1985810/1‑128 (MQ=255)
                                                                                           aTGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGtt                                                      <  2:2105602/150‑1 (MQ=255)
                                                                                               tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGa                                                                                                                                           >  1:188571/1‑61 (MQ=255)
                                                                                               tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGa                                                                                                                                           <  2:188571/61‑1 (MQ=255)
                                                                                               tCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgcg                                                  >  1:1693897/1‑150 (MQ=255)
                                                                                               tCGAGCGCGGGCTTGGCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGAGCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTcgcg                                                  <  2:921461/150‑1 (MQ=255)
                                                                                                cGAGCTCGCGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCTCTGCGGATCGAGGCCGATCACCAGGACATGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACTGGTGGTTGTCTTGCCCACACCCCCTTTCTGGTTCGCGa                                                 <  2:1780481/150‑1 (MQ=255)
                                                                                                   gCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagccag                                                                                                            <  2:1773357/88‑1 (MQ=255)
                                                                                                   gCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGcagccag                                                                                                            >  1:1773357/1‑88 (MQ=255)
                                                                                                    cTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATa                                                                                         <  2:1158822/106‑1 (MQ=255)
                                                                                                    cTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATa                                                                                         >  1:1158822/1‑106 (MQ=255)
                                                                                                    cTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGc                                             <  2:2136943/150‑1 (MQ=255)
                                                                                                     tCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCg                                            >  2:720560/1‑150 (MQ=255)
                                                                                                       gTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGaa                                          <  1:720560/150‑1 (MQ=255)
                                                                                                         gCTTTTCCAAACCGCTGCCGATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCGGACCACGACGCTTCATCGCAGCCAGCGAGGCGGCGAGATTGATACAGGAGGATGTCTGGCCCAGCCCACCTTTCTGGTTCGCGGTTGCGAAtc                                        >  2:765121/1‑149 (MQ=255)
                                                                                                          cTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAc                                       <  1:2004625/150‑1 (MQ=255)
                                                                                                           ttGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAggc                                                                                                      <  1:130112/86‑1 (MQ=255)
                                                                                                           ttGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAcc                                      >  2:1688552/1‑150 (MQ=255)
                                                                                                             gTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCtt                                    >  2:1357827/1‑150 (MQ=255)
                                                                                                             gTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCtt                                    >  2:2281683/1‑150 (MQ=255)
                                                                                                              tCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTa                                   <  1:1357827/150‑1 (MQ=255)
                                                                                                              tCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTa                                   <  1:1476993/150‑1 (MQ=255)
                                                                                                                 acacCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCGGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGcc                                >  1:221383/1‑150 (MQ=255)
                                                                                                                          cccATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCg                       >  2:1786478/1‑150 (MQ=255)
                                                                                                                            cATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAg                     >  1:1313379/1‑150 (MQ=255)
                                                                                                                                  ggCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTcc               <  2:103368/150‑1 (MQ=255)
                                                                                                                                    cGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCa             <  2:1693897/150‑1 (MQ=255)
                                                                                                                                         tccTGCTGATCGACCTCGATCCGCAGGGCACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGTTGGTTGTCTTGcc                                                                        <  2:1418937/85‑1 (MQ=25)
                                                                                                                                         tccTGCTGATCGACCTCGATCCGCAGGGCACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGcc                                                                        >  1:1418937/2‑86 (MQ=25)
                                                                                                                                          ccTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATg       <  2:111306/150‑1 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTg                                                                          >  2:2088443/1‑82 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgtt                          >  2:2353698/1‑130 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCAttgtt                          <  1:2353698/130‑1 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGc      >  2:1857315/1‑150 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGc      <  1:948633/150‑1 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTACGAGTGTTCCCAGTCATGc      >  2:1586332/1‑150 (MQ=255)
                                                                                                                                           cTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTATc                                                                             <  2:572389/79‑1 (MQ=255)
                                                                                                                                            tGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAc                                       >  2:761775/1‑116 (MQ=255)
                                                                                                                                            tGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATAc                                       <  1:761775/116‑1 (MQ=255)
                                                                                                                                             gCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGcgt    >  1:375571/1‑148 (MQ=255)
                                                                                                                                               ggATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTg  <  2:927216/150‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                     
GGTTGGCCGGCAGCAACTGGAAGCCTCCGTGCTCGGAGTAGTGCATAGCCTGGGCCAAGTCGCACTCCCCGATGAGCAGGTCGTAGACCGAATGCTCGAGCTCGTGCTTGTCCACACCGCTGCCCATGGTGGCGTTGCCCTGCGGATCGAGGTCGATCAGCAGGACACGACGCTTGGTCGCAGCCAGCGAGGCGGCGAGATTGATACAGGTGGTTGTCTTGCCCACACCACCTTTCTGGTTCGCGATTGCGAATACCTTAGCCATTGTTGCGAGTGTTCCCAGTCATGCCTTG  >  NC_002947/1556‑1848

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: