Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,054,447 A→T 6.1% K102* (AAG→TAG) ‡ spoT → bifunctional (p)ppGpp synthase/hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029476,054,4470AT 6.1% 938.7 / 34.2 278K102* (AAG→TAG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (147/114);  new base T (11/6);  total (158/120)
Fisher's exact test for biased strand distribution p-value = 6.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 281 total aligned reads displayed.

GGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTCT  >  NC_002947/6054299‑6054596
                                                                                                                                                    |                                                                                                                                                     
  ccAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGc                                                                                                                                                    <  2:2704772/150‑1 (MQ=255)
   cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                   >  2:1984896/1‑150 (MQ=255)
   cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                   >  1:3659191/1‑150 (MQ=255)
      cATCCTCGCCGACATGCACCTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGAc                                                                                                                                                <  2:5059968/150‑1 (MQ=255)
       aTCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGAcc                                                                                                                                               <  2:13959/150‑1 (MQ=255)
        tCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAccc                                                                                                                                              <  1:3846663/150‑1 (MQ=255)
          cTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAcccca                                                                                                                                            <  2:1216146/150‑3 (MQ=255)
           tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCatc                                                                                                                                           <  1:1511215/150‑3 (MQ=255)
           tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGa                                                                                                                                           <  1:1149646/150‑1 (MQ=255)
                  cATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAAttt                                                                                                                                    <  1:1856695/150‑1 (MQ=255)
                   aTGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTc                                                                                                                                   <  1:1665643/150‑1 (MQ=255)
                   aTGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTc                                                                                                                                   <  1:404247/150‑1 (MQ=255)
                    tGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCg                                                                                                                                  <  2:5017776/150‑1 (MQ=255)
                     gCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                 >  2:3475514/1‑150 (MQ=255)
                     gCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                 <  1:2485595/150‑1 (MQ=255)
                     gCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                 <  1:3984823/150‑1 (MQ=255)
                     gCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                 <  1:2322822/150‑1 (MQ=255)
                      cacaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgag                                                                                                                                <  1:3210227/150‑1 (MQ=255)
                       acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                               <  2:904904/150‑1 (MQ=255)
                       acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                               >  2:3611631/1‑150 (MQ=255)
                        caTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                              <  2:64194/150‑1 (MQ=255)
                            gACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAg                                                                                                                          >  2:4125344/1‑150 (MQ=255)
                             aCCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAgg                                                                                                                         <  2:156777/150‑1 (MQ=255)
                             aCCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAgg                                                                                                                         <  2:5590564/150‑1 (MQ=255)
                              ccATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGc                                                                                                                        >  2:4439289/1‑150 (MQ=255)
                                 tCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGa                                                                                                                     >  1:3113221/1‑150 (MQ=255)
                                   agagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAgg                                                                                                                   >  2:2515062/1‑150 (MQ=255)
                                    gagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACCCCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTGGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgc                                                                                                                  <  2:4029897/150‑1 (MQ=255)
                                    gagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgc                                                                                                                  >  1:4972276/1‑150 (MQ=255)
                                        cTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAg                                                                                                              >  2:3921994/1‑150 (MQ=255)
                                         tGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGAGGGGCTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGTGGCGCAgg                                                                                                             >  1:639779/1‑150 (MQ=255)
                                          gATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGc                                                                                                            >  1:3693499/1‑150 (MQ=255)
                                          gATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGc                                                                                                            >  1:720195/1‑150 (MQ=255)
                                                ggcCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGaaaa                                                                                                      >  1:4137293/1‑150 (MQ=255)
                                                 gcCATGCTGCACGAGGTGATCGACGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGAGGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAAc                                                                                                     <  2:2892464/150‑1 (MQ=255)
                                                      gctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGACGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCa                                                                                                <  2:4627300/150‑1 (MQ=255)
                                                      gctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCa                                                                                                >  2:2172469/1‑150 (MQ=255)
                                                      gctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCa                                                                                                >  2:5002787/1‑150 (MQ=255)
                                                      gctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCa                                                                                                <  1:5127383/150‑1 (MQ=255)
                                                       ctgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCag                                                                                               <  1:1293403/150‑1 (MQ=255)
                                                        tgcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCaga                                                                                              >  1:2902458/1‑150 (MQ=255)
                                                          cACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGACAACTTCCagaag                                                                                            >  2:4851384/1‑150 (MQ=255)
                                                             gATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatg                                                                                         <  2:673853/150‑1 (MQ=255)
                                                             gATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatg                                                                                         >  1:1033072/1‑150 (MQ=255)
                                                             gATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACAGTGGCAGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatg                                                                                         >  2:3157645/1‑150 (MQ=255)
                                                              aTGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatgg                                                                                        <  2:5134330/150‑1 (MQ=255)
                                                               tgtgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggc                                                                                       >  2:378127/1‑150 (MQ=255)
                                                                gtgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggcc                                                                                      >  1:3500450/1‑150 (MQ=255)
                                                                 tgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggcca                                                                                     <  1:5295741/150‑1 (MQ=255)
                                                                 tgATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggcca                                                                                     >  2:1027909/1‑150 (MQ=255)
                                                                  gATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccat                                                                                    <  2:5301481/150‑1 (MQ=255)
                                                                  gATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccat                                                                                    >  2:4314734/1‑150 (MQ=255)
                                                                    tCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatgg                                                                                  >  1:1104084/1‑150 (MQ=255)
                                                                     cGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggc                                                                                 >  2:442122/1‑150 (MQ=255)
                                                                         gACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatg                                                                             <  1:3067748/150‑1 (MQ=255)
                                                                         gACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatg                                                                             <  1:2385205/150‑1 (MQ=255)
                                                                         gACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatg                                                                             <  1:421669/150‑1 (MQ=255)
                                                                          acacCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatgg                                                                            <  1:1718544/150‑1 (MQ=255)
                                                                            acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggcc                                                                          >  1:3439093/1‑150 (MQ=255)
                                                                            acCGGCATCGCCAAGGAAGCCCACAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatcttc                                                                          <  1:5186365/150‑5 (MQ=255)
                                                                              cGGCATCGCCACGGACGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCcg                                                                        <  1:4721878/150‑1 (MQ=255)
                                                                              cGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggccat                                                                        <  1:1771749/150‑3 (MQ=255)
                                                                                gCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCcgcg                                                                      >  2:2837545/1‑150 (MQ=255)
                                                                                 cATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAACACTTCCAGAAGATGGCCATGGCCATGGCCCGCGa                                                                     >  2:2954698/1‑150 (MQ=255)
                                                                                          ggAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGt                                                            >  1:5110068/1‑150 (MQ=255)
                                                                                           gAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATTGCCCGCGATATCCGCGTg                                                           >  1:3323510/1‑150 (MQ=255)
                                                                                           gAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTg                                                           >  2:5508143/1‑150 (MQ=255)
                                                                                            aaGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGa                                                          >  1:3618334/1‑150 (MQ=255)
                                                                                               cccTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATcc                                                       >  1:1397810/1‑150 (MQ=255)
                                                                                                 cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                                     <  2:393401/150‑1 (MQ=255)
                                                                                                    agccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGAAGGCCATCGCCATGGCCCGCGATATACGCGTGATCCTGGTc                                                  >  2:2276319/1‑150 (MQ=255)
                                                                                                      ccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTcc                                                                                                 >  2:5372595/1‑101 (MQ=255)
                                                                                                      ccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTcc                                                                                                 <  1:5372595/101‑1 (MQ=255)
                                                                                                      ccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCaa                                                >  1:2958877/1‑150 (MQ=255)
                                                                                                         gcagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCt                                             <  1:3570157/150‑1 (MQ=255)
                                                                                                            gTTTGGCGAGACCGTGGCCTAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGc                                          <  1:247199/150‑1 (MQ=255)
                                                                                                            gTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGc                                          <  2:1452444/150‑1 (MQ=255)
                                                                                                            gTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGc                                          <  1:553882/150‑1 (MQ=255)
                                                                                                             tttGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGcc                                         >  1:4527300/1‑150 (MQ=255)
                                                                                                             tttGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGcc                                         <  2:1778020/150‑1 (MQ=255)
                                                                                                             tttGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGcc                                         >  1:2863524/1‑150 (MQ=255)
                                                                                                                 gCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGAcc                                     >  1:1947582/1‑150 (MQ=255)
                                                                                                                              cGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACAt                        >  1:1002780/1‑150 (MQ=255)
                                                                                                                               gAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATg                       >  1:3250860/1‑150 (MQ=255)
                                                                                                                                taTTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATgc                      <  1:4164958/149‑1 (MQ=255)
                                                                                                                                     gTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCAcc                 >  2:1039468/1‑150 (MQ=255)
                                                                                                                                     gTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCAcc                 >  2:3207641/1‑150 (MQ=255)
                                                                                                                                      tCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCAccc                >  1:4693062/1‑150 (MQ=255)
                                                                                                                                         tctgggTCAGCTTGCTGACCCCGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTgg             >  1:894067/4‑150 (MQ=255)
                                                                                                                                         tctgggTCAGCTTGCTGACCCCGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTgg             >  2:1169460/4‑150 (MQ=255)
                                                                                                                                         tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTgg             >  1:1679212/4‑150 (MQ=255)
                                                                                                                                         tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTgg             >  1:4350019/4‑150 (MQ=255)
                                                                                                                                            gggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAg          >  1:2111143/1‑150 (MQ=255)
                                                                                                                                             ggTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGt         >  2:3411131/1‑150 (MQ=255)
                                                                                                                                             ggTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGt         >  2:720039/1‑150 (MQ=255)
                                                                                                                                               tCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGc       <  1:4801054/150‑1 (MQ=255)
                                                                                                                                               tCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGc       >  1:748317/1‑150 (MQ=255)
                                                                                                                                                cAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCt      <  1:4439289/150‑1 (MQ=255)
                                                                                                                                                 aGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTg     >  2:3697766/1‑150 (MQ=255)
                                                                                                                                                    aaGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTCt  <  1:3403069/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GGCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTCT  >  NC_002947/6054299‑6054596

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