Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,499 A→T 7.6% K102M (AAG→ATG) ‡ 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,4990AT 7.6% 376.3 / 16.5 118K102M (AAG→ATG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (49/60);  new base T (5/4);  total (54/64)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

Only 100 of 119 total aligned reads displayed.

CCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTCTG  >  NC_002947/6049352‑6049648
                                                                                                                                                   |                                                                                                                                                     
ccAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGc                                                                                                                                                     <  2:1305430/150‑1 (MQ=255)
 cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                    >  1:427738/1‑150 (MQ=255)
 cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                    >  1:114412/1‑150 (MQ=255)
 cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCACGCAGTTTGGCGAGACCGTGTCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                    >  1:145767/1‑150 (MQ=255)
   gCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGATACCGGTATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGa                                                                                                                                                  <  2:2373503/150‑1 (MQ=255)
   gCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGa                                                                                                                                                  <  1:1522249/150‑1 (MQ=255)
     aTCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAcc                                                                                                                                                <  2:2087017/150‑1 (MQ=255)
      tCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAccc                                                                                                                                               >  1:1583646/1‑150 (MQ=255)
         tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAccccat                                                                                                                                            >  2:660854/1‑147 (MQ=255)
         tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGa                                                                                                                                            >  1:1694191/1‑150 (MQ=255)
         tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGa                                                                                                                                            <  2:853540/150‑1 (MQ=255)
          cGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAccccatc                                                                                                                                           <  1:660854/150‑5 (MQ=255)
              gACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGCCAGCAAGCTGACCCAGATGAAt                                                                                                                                       >  1:1011001/1‑150 (MQ=255)
              gACATGAACATGGACCATCCGAGCGTGATGGCGCTCACGCTGCACGACGTGATCGAAGTCAGCGGCAACGGCAAGGAAGCCCTCAGGCAGCAGGTTGGCGAGACCGTGGCCGAATTGGTCGCTGAGGTCAGCAAGCGGACCCAGCTGGAt                                                                                                                                       <  2:659059/150‑1 (MQ=255)
                cATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAAttt                                                                                                                                     >  1:2313643/1‑150 (MQ=255)
                 aTGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTc                                                                                                                                    <  2:2136020/150‑1 (MQ=255)
                   gCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                  <  1:130445/150‑1 (MQ=255)
                     acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                                <  1:1240095/150‑1 (MQ=255)
                       aTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACATCGGCATCGCCAAGGAAGCCCTCAGCACGAAAATGGGTGAGCGCGGCGCCGTGGGGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAcc                                                                                                                              >  1:2090355/1‑150 (MQ=255)
                         ggACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATTGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCaa                                                                                                                            >  1:267268/1‑150 (MQ=255)
                         ggACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCaa                                                                                                                            <  2:1949505/150‑1 (MQ=255)
                          gACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGCTGGGGCCAGCAAGCTGACCCAGCTGAATTTCGAGACCAAg                                                                                                                           >  2:1774103/1‑150 (MQ=255)
                           aCCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAgg                                                                                                                          <  2:2290040/150‑1 (MQ=255)
                           aCCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAgg                                                                                                                          <  1:702517/150‑1 (MQ=255)
                              aTCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCatc                                                                                                                                            >  2:1262111/1‑127 (MQ=255)
                              aTCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCg                                                                                                                       >  1:91388/1‑150 (MQ=255)
                              aTCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCg                                                                                                                       <  2:2192950/150‑1 (MQ=255)
                                  gagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgc                                                                                                                   >  1:1084637/1‑150 (MQ=255)
                                  gagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAAACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgc                                                                                                                   <  1:1785272/150‑1 (MQ=255)
                                     ccTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCa                                                                                                                <  1:1461233/150‑1 (MQ=255)
                                        gATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGc                                                                                                             >  1:374243/1‑150 (MQ=255)
                                            gcggcCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGaa                                                                                                         >  2:1562529/1‑150 (MQ=255)
                                               gcCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAAc                                                                                                      <  2:1694191/150‑1 (MQ=255)
                                                 cATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                                <  1:1280377/122‑1 (MQ=255)
                                                 cATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                                >  2:1280377/1‑122 (MQ=255)
                                                 cATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAgg                                                                                                                          >  1:972066/1‑128 (MQ=255)
                                                 cATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCCGCCggctagaac                                                                                                    >  2:1118380/1‑143 (MQ=255)
                                                  aTGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTc                                                                                                   >  1:1737569/1‑150 (MQ=255)
                                                  aTGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCAGGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTc                                                                                                   <  2:417023/150‑1 (MQ=255)
                                                   tgctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTcc                                                                                                  <  1:502109/150‑1 (MQ=255)
                                                   tgctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTcc                                                                                                  <  2:2266951/150‑1 (MQ=255)
                                                    ggtgcATGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCATCAGCCAGGAGTTTGGCGAGACCGTGGCCGACTTGGTCGATGGGGTCAGCAAGCTGACCCAGATGACTTTGGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCa                                                                                                 <  2:1011001/148‑1 (MQ=255)
                                                       gcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                              >  1:192404/1‑150 (MQ=255)
                                                       gcACGATGTGATCGAAGACACAGGCATCGCCAAGGAAGAACTCAGCCAGGAGGTTGGCGAGACCGTGGCCGAATGGGTCGATTGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                              <  2:267268/150‑1 (MQ=255)
                                                         aCGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaaga                                                                                            <  1:603702/150‑1 (MQ=255)
                                                            aTGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatgg                                                                                         <  2:1737569/150‑1 (MQ=255)
                                                             tgtgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggc                                                                                        <  2:192404/150‑1 (MQ=255)
                                                             tgtgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggc                                                                                        >  1:2166127/1‑150 (MQ=255)
                                                                  tCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatgg                                                                                   >  1:1831967/1‑150 (MQ=255)
                                                                     aaGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggcca                                                                                >  1:2268150/1‑150 (MQ=255)
                                                                        ccacCGGCATCGCAAAGGCAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatgg                                                                             >  2:2331868/2‑150 (MQ=255)
                                                                        acacCGTCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatgg                                                                             <  2:374243/150‑1 (MQ=255)
                                                                        acacCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatgg                                                                             <  2:1084637/150‑1 (MQ=255)
                                                                        acacCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatgg                                                                             <  1:2353831/150‑1 (MQ=255)
                                                                         cacCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggc                                                                            <  2:2313643/150‑1 (MQ=255)
                                                                          acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggcc                                                                           >  1:585387/1‑150 (MQ=255)
                                                                            cGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCcg                                                                         >  2:1840732/1‑150 (MQ=255)
                                                                             ggCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggccatc                                                                        <  2:91388/150‑4 (MQ=255)
                                                                               cATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTc                                                                                                   >  2:1625732/1‑121 (MQ=255)
                                                                               cATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggccatctt                                                                      <  1:1774103/150‑6 (MQ=255)
                                                                                aTCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                               >  1:693862/1‑92 (MQ=255)
                                                                                aTCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                               <  2:693862/92‑1 (MQ=255)
                                                                                     cAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                               >  2:2043535/1‑87 (MQ=255)
                                                                                     cAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                               <  1:2043535/87‑1 (MQ=255)
                                                                                     cAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcg                                                                >  2:94861/1‑150 (MQ=255)
                                                                                      aaGGAAGCCCTCAGCCAGCAGATTGGCGAGAAGTTGCTGGAAATGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcgc                                                               <  2:1831967/150‑1 (MQ=255)
                                                                                            gCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACt                                                                                                     <  2:820004/106‑1 (MQ=255)
                                                                                            gCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACt                                                                                                     >  1:820004/1‑106 (MQ=255)
                                                                                            gCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATc                                                         <  1:1562529/150‑1 (MQ=255)
                                                                                            gCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATc                                                         <  2:2090355/150‑1 (MQ=255)
                                                                                              ccTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatgg                                                                                         <  1:1261703/116‑1 (MQ=255)
                                                                                               cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                                      <  1:849477/150‑1 (MQ=255)
                                                                                               cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                                      <  1:306029/150‑1 (MQ=255)
                                                                                               cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                                      >  2:849477/1‑150 (MQ=255)
                                                                                               cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGACGATGGGATCAGCAAGCTGACCCAGATGAATTTCGAGACCATAATAGTGCCGACTGCTTAAAAGTTCCAGAAGATGGGCATGGCCATGGCCCGCGGTATCCGCGTGAGCCTg                                                      >  1:549395/1‑150 (MQ=255)
                                                                                                tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTTGAGACCAAGGCCGAGGAGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTgg                                                     <  2:145767/150‑1 (MQ=255)
                                                                                                tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTgg                                                     >  2:1376016/1‑150 (MQ=255)
                                                                                                tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTgg                                                     <  2:114706/150‑1 (MQ=255)
                                                                                                tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTgg                                                     <  2:114412/150‑1 (MQ=255)
                                                                                                tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTgg                                                     <  1:1103720/150‑1 (MQ=255)
                                                                                                  agccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTc                                                   >  2:971896/1‑150 (MQ=255)
                                                                                                     cagcagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGt                                                             >  2:2119142/1‑137 (MQ=255)
                                                                                                     cagcagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAg                                                >  1:1310416/1‑150 (MQ=255)
                                                                                                         agTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCCACACCAAGGCCGAGGGGCAGGCGGAAAACTTCCAGCAGCTGGCCCTGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTgg                                            <  1:1376016/150‑1 (MQ=255)
                                                                                                           tttGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggc                                                                            >  1:36846/1‑116 (MQ=255)
                                                                                                           tttGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggc                                                                            <  2:36846/116‑1 (MQ=255)
                                                                                                              ggCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGAc                                       <  1:971896/150‑1 (MQ=255)
                                                                                                                 gagaCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGc                                    >  2:1923909/1‑150 (MQ=255)
                                                                                                                      cGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatgg                                                                                   >  1:2344437/1‑98 (MQ=255)
                                                                                                                       gTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGcac                              >  1:1809546/1‑150 (MQ=255)
                                                                                                                         ggCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatg                                                                              <  1:1938269/100‑1 (MQ=255)
                                                                                                                            cGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGTCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGGTGGCCGACCGCCTGCACAACAt                         >  1:508503/1‑150 (MQ=255)
                                                                                                                                ttGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTcc                                                                                                  >  2:1900572/1‑73 (MQ=255)
                                                                                                                                   gTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggccatc                                                                        >  1:76738/1‑93 (MQ=255)
                                                                                                                                   gTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggccatc                                                                        <  2:76738/96‑4 (MQ=255)
                                                                                                                                    tCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGCCAAGCTGGCCGACCGCCTGCACAACATGCGCAccc                 <  1:94861/150‑1 (MQ=255)
                                                                                                                                       tctgggTCAGCTTGCTGACCGCGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTTCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCAGCCCGCCTGCACCACATGCGCACCCTgg              >  2:171602/4‑150 (MQ=255)
                                                                                                                                           ggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCACATGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGt          >  1:1366110/1‑150 (MQ=255)
                                                                                                                                                gCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGt     >  1:558657/1‑150 (MQ=255)
                                                                                                                                                   aGCTGACCCAGATGAATTTCGAGACCAAGGCCGACGCCCCGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTCTg  <  1:1933765/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTCTG  >  NC_002947/6049352‑6049648

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