Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,055,789 A→G 7.8% E549G (GAA→GGA)  spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,055,7890AG 7.8% 631.1 / 10.4 216E549G (GAA→GGA) spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (63/123);  new base G (15/2);  total (91/125)
Fisher's exact test for biased strand distribution p-value = 1.51e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 219 total aligned reads displayed.

AGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAAC  >  NC_002947/6055641‑6055937
                                                                                                                                                    |                                                                                                                                                     
agcagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGaa                                                                                                                                                      <  1:3897424/150‑1 (MQ=255)
agcagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGaa                                                                                                                                                      <  1:2678307/150‑1 (MQ=255)
 gcagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAAc                                                                                                                                                     <  1:7166552/150‑1 (MQ=255)
  cagcCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGGCGCCTGCTGTGGGCCGAAGGGGaaaa                                                                                                                                                    >  2:3327224/1‑148 (MQ=255)
           gaAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGgc                                                                                                                                           <  1:1112539/150‑1 (MQ=255)
             aaaTCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGGCAGCTGTTCGCCTTCGGCCGACAGCAGGCGgcgc                                                                                                                                         <  1:2890694/150‑1 (MQ=37)
               aTCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCa                                                                                                                                       <  2:973999/150‑1 (MQ=255)
                  ccccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCGGCTGTTGGGCGTAGGCGGAAAGCTGGACGCGGCAGaa                                                                                                                                    >  2:2158712/1‑150 (MQ=255)
                  ccccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCGGTAGGCCCAAGGCGAAAACCAGCCGGCGccccaa                                                                                                                                    >  1:2416041/1‑146 (MQ=255)
                   cccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAg                                                                                                                                   <  1:4044623/150‑1 (MQ=255)
                    ccAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCGGCGGTCGGCCGCAGGCGCACAGCAGCCCGCGcccccagc                                                                                                                                  >  1:2030770/1‑144 (MQ=255)
                     cAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCTGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGc                                                                                                                                 <  1:6473329/150‑1 (MQ=255)
                     cAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGc                                                                                                                                 <  1:2042382/150‑1 (MQ=255)
                      aGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGcc                                                                                                                                <  2:2887381/150‑1 (MQ=255)
                         aaCGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCAc                                                                                                                             <  1:1154060/150‑1 (MQ=255)
                          aCGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACt                                                                                                                            <  1:5732097/150‑1 (MQ=255)
                          aCGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACt                                                                                                                            <  1:5543733/150‑1 (MQ=255)
                           cGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCCCCGGCTGTCGGGCGACGGCGCACACCGGGCGGCGGCCCAAGGcccccag                                                                                                                           >  2:3347093/1‑146 (MQ=255)
                             cATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCAGCTGTCGGCCCCACGCGAACAGCGGCCGGCGCCCGAACGCCCACAGGc                                                                                                                         >  2:5251248/1‑150 (MQ=255)
                              aTCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCg                                                                                                                        <  1:3345081/150‑1 (MQ=255)
                              aTCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCAGCGATGGGCCGAAGGCGAACAGCTGACGGCGCCACAATGCCCAcagtcg                                                                                                                        >  2:778044/1‑145 (MQ=255)
                                ccAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCTGCTGGCGGCCGAAGGGGAACAGCTGCCGGCGCCCCAAGGCCCACAGGCGAt                                                                                                                      >  2:2253224/1‑150 (MQ=255)
                                ccAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCCGTCGGCCGCAGGCGCACACCAGCCCGCGCCAGAAGGCCCACAGGCGAt                                                                                                                      >  2:7162614/1‑150 (MQ=255)
                                ccAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCAGCAGGCGGCCGCACGCGCACACCAGCCGGCGCCCGCACGCCCCCAGGCgga                                                                                                                      >  2:3424210/1‑148 (MQ=255)
                                 cAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCTGCGGTCGGCCCAAGGCGAACAGCAGCCCGCGCCCCAAGGCCCACAGGCGATc                                                                                                                     >  1:3655019/1‑150 (MQ=255)
                                    tCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCcgc                                                                                                                  <  1:4994145/150‑1 (MQ=255)
                                      tATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAAc                                                                                                                                                     >  1:897154/1‑113 (MQ=255)
                                       aTTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCAGCAGTCGGCCCGAGGGGAAAAGCAGCACGCGCCCGAATTGCCACATGCCATCCGCGGc                                                                                                               >  1:608628/1‑150 (MQ=255)
                                          ctcGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCAcc                                                                                                            <  1:6615767/150‑1 (MQ=255)
                                           tcGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCCCGGTCGGCCGCAGGCGCACACCAGCCCGCGCCCCAAGGCCCCCACGCGCTCCCCGGccccc                                                                                                           >  2:4979352/1‑149 (MQ=255)
                                             gCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGaa                                                                                                         <  2:1560449/150‑1 (MQ=255)
                                                  gTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCt                                                                                                    <  1:6823510/150‑1 (MQ=255)
                                                      cGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCCCCAGCGGTCGGCCGAAGGCGGAAACCAGCCGGCGCCCCAAAGCCCCCAGGCGATACGCGGCCCCCAAGGCCTGGTg                                                                                                >  1:6165900/1‑150 (MQ=255)
                                                      cGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCGGTAGGCCCCACGCGAACACCAGCCGGCGCCCCAAAGCCCACAGGCGATCCCCGGCACCCCAGGCCAGGTg                                                                                                >  1:6166157/1‑150 (MQ=255)
                                                       gCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGGCTGCTGTGGGGCGAAGGCGTACAGGTTGGGGGGCCAGTTGGCGCACTTGGGTTAAGCGGGACCGAAGTCCTTGttg                                                                                               >  2:6654970/1‑148 (MQ=255)
                                                       gCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGc                                                                                               <  1:4870169/150‑1 (MQ=255)
                                                       gCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGc                                                                                               <  1:4825319/150‑1 (MQ=255)
                                                        ccTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCGGTCGGCCCCACGCCGCCCCCGGCCCGCGCCCCAAGGCCCCCAGGCGATCCGCCGCCCCCCAGGCCCGGTggg                                                                                              >  2:4409946/1‑148 (MQ=255)
                                                         cTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTc                                                                                             <  2:2207148/150‑1 (MQ=255)
                                                          tGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCa                                                                                            <  1:4101887/150‑1 (MQ=255)
                                                          tGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCa                                                                                            <  2:7280449/150‑1 (MQ=255)
                                                          tGGAGCTCATAGAAGACCGGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCa                                                                                            <  2:2940017/150‑1 (MQ=255)
                                                            gAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCGGCGGTCGGCCGAAGGCGGACAGCAGCCCGCCCCCGAAGGCCCACAGGCGCGCCCCGGCCCCCCAGGCCGGGTGGGCAGc                                                                                          >  2:4142672/1‑150 (MQ=255)
                                                                tCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCCCCCGCGGGCCGCCCCAGGCGAACAGCAGCCCGCGCCACAAGGCCCACAGGCGAAACGCCGCCCCCAACGCCAGGGGCGCAGcaaac                                                                                      >  2:3307879/1‑146 (MQ=255)
                                                                         aCCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGtgct                                                                             <  2:2802727/150‑1 (MQ=255)
                                                                          ccTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGtgctg                                                                            <  1:1342759/150‑1 (MQ=255)
                                                                          ccTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCTGTCGGCCGGAGGCGAACACCTGCCGGCGCCCGAAGGCCCACAGGCGATCCGCGGCCCCCAAGGCCAGGTGCTCAGCTACCCCAAGttctg                                                                            >  2:4249168/1‑150 (MQ=255)
                                                                           cTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGtgctgc                                                                           <  1:1124322/150‑1 (MQ=255)
                                                                               tCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGc                                                                       <  2:6006089/150‑1 (MQ=255)
                                                                                      cATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATccc                                                               <  2:2467890/150‑1 (MQ=255)
                                                                                      cATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATccc                                                               <  1:5448383/150‑1 (MQ=255)
                                                                                       aTCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGGGCCGCCCGCGGTCGGGCGCAGGCGGACACCTGCCGGCGCCCCAAGGCCCCCAGGCGCTCCCCGGCCCCCCAGGCCGGGGGCTCAGCGAACCCAACTGCGGCACGccccacccg                                                               >  2:5084035/1‑143 (MQ=255)
                                                                                       aTCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCg                                                              <  1:3179074/150‑1 (MQ=255)
                                                                                        tCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCgg                                                             <  1:357444/150‑1 (MQ=255)
                                                                                        tCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCgg                                                             <  2:6377759/150‑1 (MQ=255)
                                                                                        tCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCgg                                                             <  2:5500540/150‑1 (MQ=255)
                                                                                              tGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCGGGGGGCCGAAGGCGAACAGCGGCCGGCGCCCGAACGCCCCCGGGCGCTCCGCGGCCCCCCACGCCGGGTGCGCCGCGTCCCCCAGGGCTGCACCCCCCTCCCggggggcc                                                       >  2:809544/1‑145 (MQ=255)
                                                                                               gggCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCGGGCGGCCGGAGGCGGACACCAGCCCGCGCCCCAAGGCCCACGGGCGCTCCCCGGCACCCCAGGCCGGGTGGTCAGCAACCCCCAGTGCTGCACCCCGCACCCCgggggccc                                                       >  2:5359666/1‑145 (MQ=255)
                                                                                                 gCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGa                                                    <  2:1542365/150‑1 (MQ=255)
                                                                                                 gCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGa                                                    <  1:7400469/150‑1 (MQ=255)
                                                                                                 gCAACCGCATGGCCTACGTGGTCGCGCGCCGCCGGCTGTGGGGCGAAGGCGAACCGCTGGCGGCGCCAGCAGGCCCCCGGGCGCGCCGCGGCCCCCAAGGCCAGGGGCGCAGCGACGCCCAGGGCGGCCCGCCCAGCCCCGGTGTccccc                                                    >  2:6729060/1‑148 (MQ=255)
                                                                                                   aaCCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGAtt                                                  <  2:1639559/150‑1 (MQ=255)
                                                                                                   aaCCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGAtt                                                  <  2:552821/150‑1 (MQ=255)
                                                                                                   aaCCGCATGGCCTACGTGGTCGCACGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCg                                                              <  1:3321295/138‑1 (MQ=255)
                                                                                                    aCCGCATGGCCTACGTGGTCGCGGGCCCCCCGCGGTGGGCCGCCGGCGGACAGCGGGCGGCGGCACCAGGCCCACAGGCGGTGCCCCGCACCCCAGGGCCGGTGCTCAGCGAAGCCCACTGCGGGACGCCCCGCCCCGGTGGccccaatt                                                 >  2:1919209/1‑145 (MQ=255)
                                                                                                      cGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCAGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTc                                               <  1:5159431/150‑1 (MQ=255)
                                                                                                           ggCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCa                                          <  2:3250030/150‑1 (MQ=255)
                                                                                                                 cGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTc                                    <  2:2056600/150‑1 (MQ=255)
                                                                                                                    ggTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggc                                 <  2:386534/150‑1 (MQ=255)
                                                                                                                      tCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccg                               <  1:4634858/150‑1 (MQ=255)
                                                                                                                      tCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccg                               >  2:5487721/1‑150 (MQ=255)
                                                                                                                        gcgcgcCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTctgccggc                             >  1:5936674/1‑150 (MQ=255)
                                                                                                                        gcgcgcCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTcggccggc                             >  1:3361386/1‑150 (MQ=255)
                                                                                                                          gcgcCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCaa                           <  2:5920067/150‑1 (MQ=255)
                                                                                                                            gccgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCACgg                         >  2:4784231/1‑150 (MQ=255)
                                                                                                                            gccgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAgg                         <  1:4902317/150‑1 (MQ=255)
                                                                                                                             ccgccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAggg                        >  1:2025195/1‑150 (MQ=255)
                                                                                                                               gccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCa                      >  2:4739067/1‑150 (MQ=255)
                                                                                                                                ccTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCAAGAAGGCCCACTGGCGATCCTCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCAt                     >  2:1164897/1‑150 (MQ=255)
                                                                                                                                 ctgctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACGGGCGCTCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATg                    <  2:5909877/150‑1 (MQ=255)
                                                                                                                                  tgctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATgg                   >  1:1368998/1‑150 (MQ=255)
                                                                                                                                   gctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACAGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGt                  <  2:3060552/150‑1 (MQ=255)
                                                                                                                                    ctgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTc                 <  2:448770/150‑1 (MQ=255)
                                                                                                                                     tgTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCg                >  2:2714609/1‑150 (MQ=255)
                                                                                                                                      gTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTTGGCCGGCAAGGGCATGGTCGt               >  1:3015937/1‑150 (MQ=255)
                                                                                                                                      gTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGt               <  1:6867577/150‑1 (MQ=255)
                                                                                                                                       tCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTg              >  1:1914265/1‑150 (MQ=255)
                                                                                                                                         ggCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCa            >  2:739026/1‑150 (MQ=255)
                                                                                                                                         ggCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCa            >  2:1872909/1‑150 (MQ=255)
                                                                                                                                          gCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCAc           <  2:7056480/150‑1 (MQ=255)
                                                                                                                                           ccGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCAcc          >  1:7045614/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt         <  1:1004509/150‑1 (MQ=255)
                                                                                                                                            cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt         >  2:6140404/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt         >  1:6539866/1‑150 (MQ=255)
                                                                                                                                            cGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCt         >  2:5165682/1‑150 (MQ=255)
                                                                                                                                             gAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTg        >  1:523505/1‑150 (MQ=255)
                                                                                                                                                ggCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGgag     >  2:5549368/1‑150 (MQ=255)
                                                                                                                                                 gCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGt                                                           <  1:6771957/95‑1 (MQ=255)
                                                                                                                                                   gAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAAc  >  2:6928425/1‑150 (MQ=255)
                                                                                                                                                   gAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAAc  >  2:1501060/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
AGCAGCCTGGAGAAAATCCCCCAGGAACGCATCCAGTCTATTCTCGCCGAGTACCGCCTGGAGCTCATAGAAGACCTGCTCGAAGACATCGGCCTGGGCAACCGCATGGCCTACGTGGTCGCGCGCCGCCTGCTGTCGGCCGAAGGCGAACAGCTGCCGGCGCCAGAAGGCCCACTGGCGATCCGCGGCACCGAAGGCCTGGTGCTCAGCTACGCCAAGTGCTGCACGCCGATCCCGGGTGACCCGATTGTCGGCCACCTGTCGGCCGGCAAGGGCATGGTCGTGCACCTGGAGAAC  >  NC_002947/6055641‑6055937

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