Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,349 T→C 6.7% V52A (GTG→GCG) ‡ spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,049,3490TC 6.7% 448.4 / 13.9 149V52A (GTG→GCG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base T (69/70);  new base C (5/5);  total (74/75)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.90e-01

Only 100 of 150 total aligned reads displayed.

GGGTATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAG  >  NC_002947/6049200‑6049496
                                                                                                                                                     |                                                                                                                                                    
gggTATAGAAGCCTTGGCCGACCGGCTTTCGCCCTATCTTGGCCCCGAACAGGGCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGt                                                                                                                                                      <  2:1933178/150‑1 (MQ=255)
  gTATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTgg                                                                                                                                                    >  1:2290040/1‑150 (MQ=255)
  gTATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTgg                                                                                                                                                    <  1:1588004/150‑1 (MQ=255)
  gTATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTgg                                                                                                                                                    <  1:1962308/150‑1 (MQ=255)
   tataGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGc                                                                                                                                                   <  1:2049050/150‑1 (MQ=255)
     taGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCa                                                                                                                                                 <  1:1180569/150‑1 (MQ=255)
     taGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCa                                                                                                                                                 >  2:494204/1‑150 (MQ=255)
         aGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgg                                                                                                                                             <  1:239675/150‑3 (MQ=255)
         aGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgg                                                                                                                                             <  2:830321/150‑3 (MQ=255)
          gCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGNCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATc                                                                                                                                            <  1:1896885/150‑1 (MQ=255)
            cTTGGCCGAACGGCTTTCGCCACCTCAGTTCCCCGAACAAGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgggtg                                                                                                                                          <  1:1118/150‑6 (MQ=255)
             ttGGCCGAACGGCTTTCGACTTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTGCGCCGAACAGGCCCACGATTGGCAACGCCGCCGCAGTGGCGAGCCTTGCGTGACCCACCCGCTGGCCGTGGCCAGCATCCTc                                                                                                                                         <  1:2368562/150‑1 (MQ=255)
             ttGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTc                                                                                                                                         >  1:791310/1‑150 (MQ=255)
             ttGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAAGCCCAACATGCGGAACGGCGCCGCAGTGGCGAGCCCTACCTGACCCACCCGCTGGCCGTGGCCAGCATCCTc                                                                                                                                         >  2:1880960/1‑150 (MQ=255)
                gCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACGCGGCCGTGGGCCGTGGCCAGCATCCTCGcc                                                                                                                                      >  1:1811029/1‑150 (MQ=255)
                  cGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGa                                                                                                                                    >  1:2087017/1‑150 (MQ=255)
                  cGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCACGGCCAGcgg                                                                                                                                             <  2:1661952/141‑3 (MQ=255)
                   gAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGAc                                                                                                                                   >  2:1841118/1‑150 (MQ=255)
                       ggCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGc                                                                                                                               <  1:1755444/150‑1 (MQ=255)
                          tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGc                                                                                                                                       >  1:38508/1‑139 (MQ=255)
                          tttCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGcaca                                                                                                                            >  2:320179/1‑150 (MQ=255)
                               aCCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGAc                                                                                                                       >  1:2192950/1‑150 (MQ=255)
                                  tATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCAt                                                                                                                    <  1:2352121/150‑1 (MQ=255)
                                    tCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCa                                                                                                                  >  1:114706/1‑150 (MQ=255)
                                      ttGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCaga                                                                                                                <  2:2288846/150‑1 (MQ=255)
                                        ggCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGc                                                                                                              <  2:820415/150‑1 (MQ=255)
                                           cccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCa                                                                                                                     >  2:1545831/1‑140 (MQ=255)
                                           cccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCa                                                                                                                     <  1:1545831/140‑1 (MQ=255)
                                           cccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTg                                                                                                           >  2:292305/1‑150 (MQ=255)
                                            ccGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGa                                                                                                          >  1:1439822/1‑150 (MQ=255)
                                                  aGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATggcgg                                                                                                    <  2:187258/150‑1 (MQ=255)
                                                    gTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACAt                                                                                                                           >  2:974700/1‑125 (MQ=255)
                                                         ccTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtgct                                                                                             >  1:717675/1‑150 (MQ=255)
                                                           tGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCAtgctgc                                                                                           >  2:569065/1‑150 (MQ=255)
                                                                 ggCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                     <  2:1012354/150‑1 (MQ=255)
                                                                 ggCGTGCCTAGTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                     <  1:336635/150‑1 (MQ=255)
                                                                    gTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGa                                                                                  <  1:494204/150‑1 (MQ=255)
                                                                     tGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGAt                                                                                 <  1:320179/150‑1 (MQ=255)
                                                                     tGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGAt                                                                                 <  2:906545/150‑1 (MQ=255)
                                                                     tGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGAt                                                                                 <  1:1061700/150‑1 (MQ=255)
                                                                     tGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGAt                                                                                 <  1:569065/150‑1 (MQ=255)
                                                                      gCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                <  1:292305/150‑1 (MQ=255)
                                                                       ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATc                                                                                                                   <  1:1829002/114‑1 (MQ=255)
                                                                       ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATc                                                                                                                   >  2:1829002/1‑114 (MQ=255)
                                                                       ccTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGAt                                                                                 <  2:2398633/148‑1 (MQ=255)
                                                                        cTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGAt                                                                                                         >  1:978313/1‑123 (MQ=255)
                                                                        cTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGAt                                                                                                         <  2:978313/123‑1 (MQ=255)
                                                                         tATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGaa                                                                             <  1:775453/150‑1 (MQ=255)
                                                                          aTTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAg                                                                            >  1:1633843/1‑150 (MQ=255)
                                                                          aTTTCTACGCCGAACAGCCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAg                                                                            >  1:2373503/1‑150 (MQ=255)
                                                                           tttCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCTCTACGTGACCCCCCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTAATGGCGGCCATGCTGCACGATGTGATCGAAGa                                                                           <  2:546128/150‑1 (MQ=255)
                                                                              cTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACAt                                                                                                                                 >  2:942165/1‑93 (MQ=255)
                                                                              cTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACAt                                                                                                                                 <  1:942165/93‑1 (MQ=255)
                                                                               tACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGc                                                                                                                                                   >  1:547189/1‑74 (MQ=255)
                                                                               tACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCAc                                                                                         >  2:57543/1‑132 (MQ=255)
                                                                               tACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACAcc                                                                       >  1:1305430/1‑150 (MQ=255)
                                                                                 cGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCgg                                                                     >  2:1096619/1‑150 (MQ=255)
                                                                                 cGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCgg                                                                     <  2:669563/150‑1 (MQ=255)
                                                                                 cGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTAACTACGCCAGCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCgg                                                                     >  1:659059/1‑150 (MQ=255)
                                                                                     gaccAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATc                                                                 >  1:2399613/4‑150 (MQ=255)
                                                                                     gAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCTTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATc                                                                                >  1:1551070/1‑135 (MQ=255)
                                                                                         aGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGc                                                                                                              >  1:1919399/1‑101 (MQ=255)
                                                                                           gCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAg                                                           <  2:170970/150‑1 (MQ=255)
                                                                                             ccACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGAtg                                                                                     >  2:2041100/1‑122 (MQ=255)
                                                                                              cACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGac                                                        >  2:1933765/1‑149 (MQ=255)
                                                                                              cACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGaa                                                        >  1:2131189/1‑150 (MQ=255)
                                                                                              cACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGaa                                                        >  1:1940158/1‑150 (MQ=255)
                                                                                               aCGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATc                                                                                                                                            <  2:391901/65‑1 (MQ=255)
                                                                                               aCGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATc                                                                                                                                            >  1:391901/1‑65 (MQ=255)
                                                                                               aCGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAggcca                                                       <  1:1667319/150‑4 (MQ=255)
                                                                                                 gATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGcc                                                     <  2:717675/150‑1 (MQ=255)
                                                                                                 gATGGGCAACGCCGCCGCAGTGGCGAGCCCTACATGACCGAGCCGGTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGcc                                                     <  2:529430/150‑1 (MQ=255)
                                                                                                     ggCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                                 <  1:1856628/150‑1 (MQ=255)
                                                                                                      gCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAAACACCGGCATCGCCAAGGAAGCCCTcag                                                >  2:1785272/1‑150 (MQ=255)
                                                                                                       cAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagc                                               >  2:1461233/1‑150 (MQ=255)
                                                                                                        aaCGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagcc                                              >  1:2035666/1‑150 (MQ=255)
                                                                                                                  cAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGGAGCCCTCAGCCAGCAGTTTgg                                   >  2:1240095/1‑150 (MQ=255)
                                                                                                                  cAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTgg                                   <  2:729958/150‑1 (MQ=255)
                                                                                                                            cccTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTg                         <  2:2095651/150‑1 (MQ=255)
                                                                                                                               tACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGcc                      >  1:2266951/1‑150 (MQ=255)
                                                                                                                               tACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGcc                      >  2:1103720/1‑150 (MQ=255)
                                                                                                                                aCGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagccagc                                           >  2:1444232/1‑129 (MQ=255)
                                                                                                                                aCGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTcagccagc                                           <  1:1444232/129‑1 (MQ=255)
                                                                                                                                    gACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGACCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAAtt                 <  2:1633843/150‑1 (MQ=255)
                                                                                                                                     ccccacccGCGGGCCGTGGCCAGCAGCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTg                <  2:2376432/149‑1 (MQ=255)
                                                                                                                                       ccacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGt              >  1:2136020/1‑150 (MQ=255)
                                                                                                                                        cacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGc                                  >  1:1669739/1‑129 (MQ=255)
                                                                                                                                        cacccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGc                                  <  2:1669739/129‑1 (MQ=255)
                                                                                                                                         acccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCgaga                              >  1:971741/1‑132 (MQ=255)
                                                                                                                                         acccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCga                                >  2:1568207/1‑130 (MQ=255)
                                                                                                                                         acccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCg            >  2:130445/1‑150 (MQ=255)
                                                                                                                                           ccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGc                                  <  2:2279939/126‑1 (MQ=255)
                                                                                                                                           ccGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGAt          <  2:2399613/150‑1 (MQ=255)
                                                                                                                                            tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATg         >  2:2353831/2‑150 (MQ=255)
                                                                                                                                            tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATg         >  2:1522249/2‑150 (MQ=255)
                                                                                                                                            tgCTGGCCACGGCCAGCATCCTCGCCGACATGCACATGGAACATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTca                                                 >  1:2298075/2‑111 (MQ=255)
                                                                                                                                              cTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAATCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATggg       <  1:441487/150‑1 (MQ=255)
                                                                                                                                              cTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATggg       <  1:1085122/150‑1 (MQ=255)
                                                                                                                                                ggCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGt     <  1:1841118/150‑1 (MQ=255)
                                                                                                                                                   cGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACAGCGGCCTCGCCAAGGAAGACCGCGACGGCTTGTAAAGCGAGACCTTGGCCGAATGGTGTGATGGGGTCAg  >  2:635344/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGGTATAGAAGCCTTGGCCGAACGGCTTTCGACCTATCTTGGCCCCGAACAGGTCAACCTGGTTCGGCGTGCCTATTTCTACGCCGAACAGGCCCACGATGGGCAACGCCGCCGCAGTGGCGAGCCCTACGTGACCCACCCGCTGGCCGTGGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAG  >  NC_002947/6049200‑6049496

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: