Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 6,049,498 A→T 6.8% 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,049,4980AT 6.8% 273.7 / 9.2 88K102* (AAG→TAG) ‡spoTbifunctional (p)ppGpp synthase/hydrolase
Reads supporting (aligned to +/- strand):  ref base A (34/48);  new base T (4/2);  total (38/50)
Fisher's exact test for biased strand distribution p-value = 3.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.18e-01

GCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTC  >  NC_002947/6049351‑6049646
                                                                                                                                                   |                                                                                                                                                    
gCCACCATCCCCGCCGACATGCACATGGGCCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAg                                                                                                                                                    <  1:1347728/150‑1 (MQ=255)
 ccAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCTGCCCGCAGCTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGc                                                                                                                                                   <  1:904431/150‑1 (MQ=255)
  cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                  >  2:1222393/1‑150 (MQ=255)
  cAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCt                                                                                                                                                  >  1:1503332/1‑150 (MQ=255)
    gCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGa                                                                                                                                                <  2:1282464/150‑1 (MQ=255)
    gCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGa                                                                                                                                                <  2:1812627/150‑1 (MQ=255)
    gCAGCCTCGCCGACATGCACATGGACCATCAGAGCCGGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGa                                                                                                                                                <  2:276398/150‑1 (MQ=255)
     cATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGAc                                                                                                                                               >  1:1500938/1‑150 (MQ=255)
      aTCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGAcc                                                                                                                                              <  1:1895512/150‑1 (MQ=255)
        ccTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCa                                                                                                                                            <  2:1670816/150‑1 (MQ=255)
         cTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCTTGCTGAcccca                                                                                                                                           <  1:1274134/150‑3 (MQ=255)
          tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGa                                                                                                                                          <  1:1551102/150‑1 (MQ=255)
          tCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATAGCCGAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACGGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGa                                                                                                                                          >  1:601564/1‑150 (MQ=255)
           cGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGAt                                                                                                                                         <  1:64818/150‑1 (MQ=255)
           cGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGAt                                                                                                                                         <  1:820175/150‑1 (MQ=255)
                    gCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                <  1:1622023/150‑1 (MQ=255)
                     cacaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTGTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgag                                                                                                                               >  2:822938/1‑150 (MQ=255)
                      acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                              >  1:874502/1‑150 (MQ=255)
                      acaTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                              <  1:828121/150‑1 (MQ=255)
                       caTGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGAc                                                                                                                             <  1:685862/150‑1 (MQ=255)
                           gACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAg                                                                                                                         <  1:1357318/150‑1 (MQ=255)
                           gACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAg                                                                                                                         >  2:557820/1‑150 (MQ=255)
                             ccATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGc                                                                                                                       <  2:1503332/150‑1 (MQ=255)
                              cATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGcc                                                                                                                      >  2:1008765/1‑150 (MQ=255)
                               aTCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCg                                                                                                                     >  1:830107/1‑150 (MQ=255)
                                tCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGa                                                                                                                    >  1:1001317/1‑150 (MQ=255)
                                tCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAg                                                                                                                   <  2:894370/150‑1 (MQ=255)
                                   gagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgc                                                                                                                 <  1:1033946/150‑1 (MQ=255)
                                   gagCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGgc                                                                                                                 >  2:380355/1‑150 (MQ=255)
                                      ccTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCa                                                                                                              <  1:735842/150‑1 (MQ=255)
                                       cTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAg                                                                                                             <  2:1583411/150‑1 (MQ=255)
                                        tGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAgg                                                                                                            >  2:756797/1‑150 (MQ=255)
                                         gagggcggcCATGCTGCACGATGTGAGCGAAGACACCGGGATCGCCAAGTAATCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGAGCCAGATGAATCCCGAGACCAAGGCCGAGGCGCAGGc                                                                                                           <  2:1552343/147‑1 (MQ=255)
                                             gaggcCATGCTGCAAGATGTGATAGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGaa                                                                                                       <  2:874502/148‑1 (MQ=255)
                                                gcCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAAc                                                                                                    >  2:532384/1‑150 (MQ=255)
                                                        gcACGATGTGATCGAAGACATCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                            <  1:1896982/150‑1 (MQ=255)
                                                        gcACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTAACCCGACCGTGGCGGCACCGGTCGAGGGGTTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                            >  1:231918/1‑150 (MQ=255)
                                                          aCGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaaga                                                                                          <  1:984287/150‑1 (MQ=255)
                                                          aCGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaaga                                                                                          <  2:1369069/150‑1 (MQ=255)
                                                           cGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGat                                                                                         <  1:3016/150‑1 (MQ=255)
                                                             aTGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatgg                                                                                       <  2:956826/150‑1 (MQ=255)
                                                                 gATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccat                                                                                   <  2:1500938/150‑1 (MQ=255)
                                                                  aTCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatg                                                                                  <  2:335004/150‑1 (MQ=255)
                                                                  aTCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCCAAAGCGGAAAACTTCCAGAAGatggccatg                                                                                  <  1:50733/150‑1 (MQ=255)
                                                                    cGAAGACACCGGCATCGCCAAGGAAGCCATCAGGCATTAGTTTATCGACTCCACGTCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggc                                                                                >  1:239674/1‑150 (MQ=255)
                                                                     gAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGCATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTc                                                                                                 <  1:1076468/132‑1 (MQ=255)
                                                                     gAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGTATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggcc                                                                               >  1:498450/1‑150 (MQ=255)
                                                                     gAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTc                                                                                                 >  2:1076468/1‑132 (MQ=255)
                                                                           acCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCTAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggccatggccatggcc                                                                         <  1:214399/150‑1 (MQ=255)
                                                                                 aTCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                <  2:1399639/89‑1 (MQ=255)
                                                                                 aTCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGGCAGCAAGCTGACCCAGATGAATTTCga                                                                                                                                >  1:1399639/1‑89 (MQ=255)
                                                                                    gCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                            >  1:299919/1‑122 (MQ=255)
                                                                                    gCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCagaa                                                                                            <  2:299919/122‑1 (MQ=255)
                                                                                      cAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcg                                                              >  1:174189/1‑150 (MQ=255)
                                                                                       aaGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcgc                                                             <  2:174189/150‑1 (MQ=255)
                                                                                        aGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCcgcg                                                            <  2:290982/150‑1 (MQ=255)
                                                                                           aaGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGa                                                         <  1:530370/150‑1 (MQ=255)
                                                                                               ccTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCt                                                     >  2:423738/1‑150 (MQ=255)
                                                                                                cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                                    <  2:479955/150‑1 (MQ=255)
                                                                                                cTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTg                                                    >  1:555450/1‑150 (MQ=255)
                                                                                                 tCAGCCAGCAGTTTGGCGGGACCGTGGCCGAATTGCTCGATGGGGTCAGCAAGCTGCCCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGCAGATGGCCATGGCCCTGGCCCGCGATATCCGCGTGATCCTgg                                                   <  2:563816/150‑1 (MQ=255)
                                                                                                 tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggcca                                                                                    >  2:1555492/1‑117 (MQ=255)
                                                                                                 tCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGatggcca                                                                                    <  1:1555492/117‑1 (MQ=255)
                                                                                                  cagccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGt                                                  <  2:1785651/150‑1 (MQ=255)
                                                                                                   agccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAGAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTc                                                 <  1:380355/150‑1 (MQ=255)
                                                                                                   agccagcAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTc                                                 >  2:1180761/1‑150 (MQ=255)
                                                                                                       agcagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGc                                             >  2:1865553/1‑150 (MQ=255)
                                                                                                       agcagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGc                                             >  2:1840290/1‑150 (MQ=255)
                                                                                                        gcagTTTAGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCt                                            <  1:557820/150‑1 (MQ=255)
                                                                                                         cagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACGCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTg                                           <  2:1001317/150‑1 (MQ=255)
                                                                                                         cagTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTg                                           >  1:1160241/1‑150 (MQ=255)
                                                                                                                   agaCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCAGAAAACTTCCAGCATATGTCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGcc                                 <  1:234381/150‑1 (MQ=255)
                                                                                                                    gaCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCt                                <  2:239674/150‑1 (MQ=255)
                                                                                                                     aCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                              >  2:1897084/1‑55 (MQ=255)
                                                                                                                     aCCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCgaga                                                                                                                              <  1:1897084/55‑1 (MQ=255)
                                                                                                                       catGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGca                             >  2:940603/3‑150 (MQ=255)
                                                                                                                           gCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCacaac                         >  2:1524590/1‑150 (MQ=255)
                                                                                                                              gAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCACGTGATCCTGGTCAAGCTGGCCGACCCCCCGCACAACATg                      >  2:265441/1‑150 (MQ=255)
                                                                                                                                        tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTgg            >  2:656607/4‑150 (MQ=255)
                                                                                                                                        tctgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTgg            >  2:317322/4‑150 (MQ=255)
                                                                                                                                          tgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGaa          <  1:656607/149‑1 (MQ=255)
                                                                                                                                          tgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGaa          >  2:1702912/2‑150 (MQ=255)
                                                                                                                                          tgggTCAGCTTGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGaa          >  2:1112402/2‑150 (MQ=255)
                                                                                                                                             gccAGCAAGCTGACCCAGCTGAATTTAGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGTGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCAAAACATGCGCACCCTGGAAGTg       >  2:963158/3‑150 (MQ=255)
                                                                                                                                              tCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGc      >  1:1208005/1‑150 (MQ=255)
                                                                                                                                                  cAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCacaac                         <  1:1264140/127‑1 (MQ=255)
                                                                                                                                                  cAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCacaac                         >  2:1264140/1‑127 (MQ=255)
                                                                                                                                                  cAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTc  <  2:555450/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
GCCAGCATCCTCGCCGACATGCACATGGACCATCAGAGCCTGATGGCGGCCATGCTGCACGATGTGATCGAAGACACCGGCATCGCCAAGGAAGCCCTCAGCCAGCAGTTTGGCGAGACCGTGGCCGAATTGGTCGATGGGGTCAGCAAGCTGACCCAGATGAATTTCGAGACCAAGGCCGAGGCGCAGGCGGAAAACTTCCAGAAGATGGCCATGGCCATGGCCCGCGATATCCGCGTGATCCTGGTCAAGCTGGCCGACCGCCTGCACAACATGCGCACCCTGGAAGTGCTGTC  >  NC_002947/6049351‑6049646

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: