Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,473,762 A→T 5.0% L130Q (CTG→CAG)  PP_3086 ← RNA polymerase sigma‑70 factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,473,7620AT 5.0% 321.4 / 5.5 99L130Q (CTG→CAG) PP_3086RNA polymerase sigma‑70 factor
Reads supporting (aligned to +/- strand):  ref base A (53/41);  new base T (3/2);  total (56/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.73e-01

Only 100 of 106 total aligned reads displayed.

CCCAGGATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTAGGCCCGCTCGACAT  >  NC_002947/3473613‑3473909
                                                                                                                                                     |                                                                                                                                                   
cccAGGATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGTAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGAcagca                                                                                                                                                     <  1:814407/150‑1 (MQ=255)
cccAGGATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGAcagca                                                                                                                                                     >  1:2108017/1‑150 (MQ=255)
    ggATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGaa                                                                                                                                                 <  1:1956371/150‑1 (MQ=255)
     gATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAac                                                                                                                                                >  1:780923/1‑150 (MQ=255)
     gATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAac                                                                                                                                                >  2:2146498/1‑150 (MQ=255)
       tCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAacac                                                                                                                                              <  1:1199963/150‑1 (MQ=255)
           atgcacAATGCGCTGCACCATCACTCGAGGCTCAGGCAGGCAAGTTATGCCTTACGCATGTCGCGCTTGATCGTGATCAGCTAAANCTGCAGNTGTTCGGNGATTNGGGCGTAGGTCAGGNCGTCCAGCTGTGACAGCAGGAACACCTGa                                                                                                                                          <  1:509876/145‑1 (MQ=255)
                cAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTAc                                                                                                                                     >  1:678896/1‑150 (MQ=255)
                        tGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGa                                                                                                                                          >  2:1111326/1‑137 (MQ=255)
                        tGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGa                                                                                                                                          <  1:1111326/137‑1 (MQ=255)
                         gCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCgg                                                                                                                            <  1:1811907/150‑1 (MQ=255)
                           aCCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGAAGGAAAACCTGACGTACCAGCGCCGTCa                                                                                                                          >  2:1960220/1‑150 (MQ=255)
                             cATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATg                                                                                                                        >  1:1794001/1‑150 (MQ=255)
                             cATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGCGACAGCAGGAACACCTGACGTACCCGCGCCGGCATg                                                                                                                        >  1:1128765/1‑150 (MQ=255)
                              aTCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCt                                                                                                                                            <  2:2252272/129‑1 (MQ=255)
                              aTCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCt                                                                                                                                            >  1:2252272/1‑129 (MQ=255)
                              aTCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGt                                                                                                                       >  2:2235365/1‑150 (MQ=255)
                              aTCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGt                                                                                                                       >  2:2391394/1‑150 (MQ=255)
                                 aCTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTcgc                                                                                                                    >  2:2328082/1‑150 (MQ=255)
                                 aCTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCGGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTcgc                                                                                                                    >  1:242298/1‑150 (MQ=255)
                                  cTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTa                                                                                                                                      <  2:1770265/131‑1 (MQ=255)
                                  cTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTa                                                                                                                                      >  1:1770265/1‑131 (MQ=255)
                                  cTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTcgcg                                                                                                                   >  1:2129572/1‑150 (MQ=255)
                                  cTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTcgcg                                                                                                                   >  2:621125/1‑150 (MQ=255)
                                    cGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCgca                                                                                                                 <  1:2146498/150‑1 (MQ=255)
                                      aTGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCgcagc                                                                                                               >  1:1676871/1‑150 (MQ=255)
                                             ggcaggcaAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGt                                                                                                        >  1:2174058/1‑150 (MQ=255)
                                                   caAGGAATGCCGTACGCATGTCGCGCTTGATCGTGAACAGCGAAAACTGCAGCTGTTCGGAGAAGTGGGAGTAGATCAGGACGTACAGCTGAGACAGCAGGAGCACCTGACGTACGCGCGCCGGCATGTCGGGCGGCATGGCGTCGATGc                                                                                                  >  2:520146/1‑150 (MQ=255)
                                                        aaggCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTAccc                                                                                                                                   <  1:257174/109‑1 (MQ=255)
                                                        aaTGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTAccc                                                                                                                                   >  2:257174/1‑112 (MQ=255)
                                                            ccGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAggg                                                                                         <  2:2174058/150‑1 (MQ=255)
                                                              gTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGtt                                                                                       >  2:705154/1‑150 (MQ=255)
                                                              gTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGACTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGtt                                                                                       <  2:2277022/150‑1 (MQ=255)
                                                               tACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGGTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGttt                                                                                      <  1:509214/150‑1 (MQ=255)
                                                                aCGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTc                                                                                     <  1:1590227/150‑1 (MQ=255)
                                                                  gCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCaa                                                                                   <  2:383461/150‑1 (MQ=255)
                                                                  gCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCaa                                                                                   <  1:528211/150‑1 (MQ=255)
                                                                      gTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGAt                                                                               <  2:481550/150‑1 (MQ=255)
                                                                       tCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATc                                                                              >  2:384989/1‑150 (MQ=255)
                                                                         gcgcTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATcag                                                                            >  1:1071536/1‑150 (MQ=255)
                                                                             ttGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATcagccag                                                                        <  2:2293762/150‑1 (MQ=255)
                                                                             ttGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATcagccag                                                                        <  2:1676871/150‑1 (MQ=255)
                                                                                 tCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCaa                                                                                   >  2:2392813/1‑135 (MQ=255)
                                                                                 tCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCaa                                                                                   <  1:2392813/135‑1 (MQ=255)
                                                                                  cGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCgcgc                                                                                                                               >  2:2281206/1‑90 (MQ=255)
                                                                                  cGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCgcgc                                                                                                                               <  1:2281206/90‑1 (MQ=255)
                                                                                   gTGATCAGCGAAGCCTGCAGCGGTTCGGCGATTTGGGCGTCGGTCAGGCCGTCCAGCTGTGACCGCAGGAACACCTGACGTACAGGCGCCAGCCTGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGtt                                                                  <  2:577099/150‑1 (MQ=255)
                                                                                   gTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGtt                                                                  <  2:139784/150‑1 (MQ=255)
                                                                                     gATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTc                                                                >  2:1154768/1‑150 (MQ=255)
                                                                                     gATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTc                                                                >  2:1610458/1‑150 (MQ=255)
                                                                                        cAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCgg                                                             >  1:2341965/1‑150 (MQ=255)
                                                                                         aGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGc                                                            >  2:354278/1‑150 (MQ=255)
                                                                                          gCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCgg                                                             >  1:1659207/1‑148 (MQ=255)
                                                                                          gCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCgg                                                             <  2:1659207/148‑1 (MQ=255)
                                                                                           cGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGa                                                          >  2:880897/1‑150 (MQ=255)
                                                                                            gAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGaa                                                         <  2:68539/150‑1 (MQ=255)
                                                                                            gAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGaa                                                         <  2:68550/150‑1 (MQ=255)
                                                                                            gAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGCCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGaa                                                         <  1:1118116/150‑1 (MQ=255)
                                                                                             aaaCCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAg                                                        >  1:1694174/1‑150 (MQ=255)
                                                                                              aaCCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAgg                                                       <  1:1246138/150‑1 (MQ=255)
                                                                                                      gCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCgg                                               <  2:2129572/150‑1 (MQ=255)
                                                                                                       cTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCggg                                              >  2:1877708/1‑150 (MQ=255)
                                                                                                          ttCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTg                                           >  1:749146/1‑150 (MQ=255)
                                                                                                          ttCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTg                                           >  2:265472/1‑150 (MQ=255)
                                                                                                            cGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCg                                         <  1:996487/150‑1 (MQ=255)
                                                                                                             gtacgtttGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCgg                                        >  1:683929/6‑150 (MQ=255)
                                                                                                             ggCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCgg                                        >  2:486917/1‑150 (MQ=255)
                                                                                                              gCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGc                                       >  2:58194/1‑150 (MQ=255)
                                                                                                                gATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGCTGCGGcca                                     >  1:491842/1‑148 (MQ=255)
                                                                                                                    tGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATgg                                 <  1:1610458/150‑1 (MQ=255)
                                                                                                                    tGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATgg                                 <  1:1877708/150‑1 (MQ=255)
                                                                                                                     gggCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATggg                                >  2:1587112/1‑150 (MQ=255)
                                                                                                                         gTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACa                                                                                            <  1:55768/86‑1 (MQ=255)
                                                                                                                         gTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACa                                                                                            >  2:55768/1‑86 (MQ=255)
                                                                                                                         gTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGTCACCTCGGGTTGCGGCAAATGGGCGAt                            <  2:1414840/150‑1 (MQ=255)
                                                                                                                         gTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGAt                            <  1:2350419/150‑1 (MQ=255)
                                                                                                                            ggTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGc                         <  1:67474/150‑1 (MQ=255)
                                                                                                                            ggTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGc                         >  2:67474/1‑150 (MQ=255)
                                                                                                                             gTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACTTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGcc                        >  2:633450/1‑150 (MQ=255)
                                                                                                                             gTCAGGCCGGCCAGCTGTGACAGCAGGAACACCTGACGTGGCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGGGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGcc                        <  2:242298/150‑1 (MQ=255)
                                                                                                                              tCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACATGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCt                       >  1:1547517/1‑150 (MQ=255)
                                                                                                                              tCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCt                       >  2:2282486/1‑150 (MQ=255)
                                                                                                                              tCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCt                       >  2:858302/1‑150 (MQ=255)
                                                                                                                                 ggCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCa                    >  1:1458404/1‑150 (MQ=255)
                                                                                                                                  gCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAg                   <  2:2108017/150‑1 (MQ=255)
                                                                                                                                     ttcctgCTGTCACAGCTGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGc                                          >  2:1136310/6‑124 (MQ=255)
                                                                                                                                     ttcctgCTGTCACAGCTGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGc                                          <  1:1136310/119‑1 (MQ=255)
                                                                                                                                     ttcatgCTGTCACAGCTGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTa                >  1:220950/6‑150 (MQ=255)
                                                                                                                                      tCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATg                           <  1:950669/138‑1 (MQ=255)
                                                                                                                                      tCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATg                           >  2:950669/1‑138 (MQ=255)
                                                                                                                                      tCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGAGATGGCCTCCAGGTAg               >  2:1217198/1‑150 (MQ=255)
                                                                                                                                       cctgCTGTCACAGCTGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTAgg              >  2:393023/4‑150 (MQ=255)
                                                                                                                                         aGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCTTCCAGGTAGGcc            <  2:1404277/150‑1 (MQ=255)
                                                                                                                                              tgACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTAGGCCCGCTc       >  1:614011/1‑150 (MQ=255)
                                                                                                                                                aCAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTAGGCCCGCTCGa     >  1:227377/1‑150 (MQ=255)
                                                                                                                                                  agcagGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTc                                                 >  2:2230458/1‑104 (MQ=255)
                                                                                                                                                  agcagGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTc                                                 <  1:2230458/104‑1 (MQ=255)
                                                                                                                                                  agcagGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGAt                            <  1:921996/125‑1 (MQ=255)
                                                                                                                                                  agcagGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGAt                            >  2:921996/1‑125 (MQ=255)
                                                                                                                                                   gcagGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTAGGCCCGCTCGACAt  >  1:1959104/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
CCCAGGATCAGCGGGTCAATGCGCTGCACCATCACTCGATGCTCAGGCAGGCAAGGAATGCCGTACGCATGTCGCGCTTGATCGTGATCAGCGAAACCTGCAGCTGTTCGGCGATTTGGGCGTAGGTCAGGCCGTCCAGCTGTGACAGCAGGAACACCTGACGTACCCGCGCCGGCATGTCGCGCAGCATGGCGTCGATGCGGTACAGGGTTTCAAGGATCAGCCAGCGCGTTTCCGGCGAAGGCACCTCGGGTTGCGGCAAATGGGCGATGGCCTCCAGGTAGGCCCGCTCGACAT  >  NC_002947/3473613‑3473909

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