Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,068 G→T 13.0% R171R (CGC→CGA PP_4277 ← GntR family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,866,0680GT13.0% 228.1 / 17.6 83R171R (CGC→CGAPP_4277GntR family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (35/37);  new base T (7/4);  total (42/41)
Fisher's exact test for biased strand distribution p-value = 5.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

GTCAGCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTGGTC  >  NC_002947/4865920‑4866217
                                                                                                                                                     |                                                                                                                                                     
gTCAGCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGgc                                                                                                                                                      >  2:929340/1‑150 (MQ=255)
gTCAGCCGCTGGCGAGTGCCGCCAGCACCAGTTCGCAGCTGCGCGGGGCTTAGGCGGTGAGCACCTCGCGAATGCTCGGCAGGTCACGCGGGACCACGGGTTCTAGCAAGCGGTCGAACAGGTCTAGATAGTCGATGTTGTTGGCCTGgc                                                                                                                                                      <  2:114113/150‑1 (MQ=255)
  cAGCCTCTGGCCAGTCCTGCCAGCACCGGATCGGAGCTGTGCTGGCCATAGGCGGTTAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCCCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGgcgc                                                                                                                                                    <  2:1128694/150‑1 (MQ=255)
    gCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTg                                                                                                                                                  <  2:688303/150‑1 (MQ=255)
     ccTCTGGCCAGTGCCGCCAGCGCCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGTGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTgg                                                                                                                                                 >  2:344914/1‑150 (MQ=255)
          ggCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAg                                                                                                                                            >  1:1356203/1‑150 (MQ=255)
             cAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCAGTGAGCACCGCGCGAATGCGCGGCAGGTCACGCCAGACCCCATACTCGAGCAAGCAGGCGAACAGGTCGAGCTAGTCGATCATGTGGGTCTGGCTCTCCTCCAGccc                                                                                                                                         >  2:1056216/1‑148 (MQ=255)
                tGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAg                                                                                                                                      <  2:1197393/150‑1 (MQ=255)
                 gcggccAGCACCAGTTCGCACCTCCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGCCCACGGCTTCGAACAAGCGGACGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAgg                                                                                                                                     <  1:433627/147‑1 (MQ=255)
                 gccgccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAgg                                                                                                                                     <  1:36967/150‑1 (MQ=255)
                 gccgccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGACCTGTCGCTGGACCAGCGCCAgg                                                                                                                                     <  1:1593493/150‑1 (MQ=255)
                     ccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAAt                                                                                                                                 <  2:957522/150‑1 (MQ=255)
                     ccAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGAGGCCGAACAGGTCGACCTCCCCGATCGATGGGGCCTGGCGCTGGTCCAGCGACAGGTAAt                                                                                                                                 >  2:1771141/1‑150 (MQ=255)
                      cAGCACCCGTTCGCAGCTGCGCTGGCCATAGGCGGTGACCACCTCGGGCATAACCGGCAGGTCACGCGAGCCCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATAATGTTGGCCTGTCGCTTGTGCAGCGagagtttatt                                                                                                                                >  2:942977/1‑144 (MQ=255)
                      cAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATa                                                                                                                                >  2:446378/1‑150 (MQ=255)
                       aGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAg                                                                                                                               >  2:1773/1‑150 (MQ=255)
                         cACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGc                                                                                                                             <  2:1841250/150‑1 (MQ=255)
                            tcatggccAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACg                                                                                                                          >  1:1623953/8‑150 (MQ=255)
                             aGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACgg                                                                                                                         >  1:2156931/1‑150 (MQ=255)
                                    aGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCGGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATgg                                                                                                                  >  1:1929593/1‑150 (MQ=255)
                                          gcTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACCGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCt                                                                                                            >  1:418883/1‑150 (MQ=255)
                                            tGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGc                                                                                                          >  1:950754/1‑150 (MQ=255)
                                            tGGACATAGGCGGTGAGCAGCTCGCGAGTGCGCGGCAGGTCACGCGAGACCACGACTGCGAGCAAGCGGGCGAACAGGTCGAGATAGTGGATCATGTTGGCCGGGCGATGGTCCAGCGACAGGTAATAGGCGCGGCTCATGGCCGGCTGc                                                                                                          >  2:166512/1‑150 (MQ=255)
                                              gCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAg                                                                                                        <  2:1052221/150‑1 (MQ=255)
                                                 aTAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCgg                                                                                                              >  2:943445/1‑141 (MQ=255)
                                                 aTAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCgg                                                                                                              <  1:943445/141‑1 (MQ=255)
                                                     gCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTcc                                                                                                 >  2:2096405/1‑150 (MQ=255)
                                                      cGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCa                                                                                                <  2:798033/150‑1 (MQ=255)
                                                         tGAGCACCTCGTGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACgg                                                                                             <  1:344914/150‑1 (MQ=255)
                                                              aCCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTcc                                                                                        >  2:251201/1‑150 (MQ=255)
                                                               ccTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCt                                                                                       <  2:418883/150‑1 (MQ=255)
                                                                cTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAgg                                                                                                                                     <  1:1182446/103‑1 (MQ=255)
                                                                cTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAgg                                                                                                                                     >  2:1182446/1‑103 (MQ=255)
                                                                 tCGCGAATGCGCGGCAGGTCACGCGAGGCCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGACCTGTCGCTGGACCAGCgcca                                                                                                                                       >  2:1293648/1‑97 (MQ=255)
                                                                 tCGCGAATGCGCGGCAGGTCACGCGAGGCCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGACCTGTCGCTGGACCAGCGCCa                                                                                                                                       <  1:1293648/100‑1 (MQ=255)
                                                                 tCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGGGGGTGAACACCCAGAGAATCTCTACCAATAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTATTCCGCTGGTGCCTGc                                                                                     >  1:782028/1‑150 (MQ=255)
                                                                  cgcgAATGCGCGGCAGGTCACGCGAGACCACGGCGTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCa                                                                                    <  1:2312290/150‑1 (MQ=255)
                                                                       aTGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAAc                                                                               <  1:446378/150‑1 (MQ=255)
                                                                        tGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCgcca                                                                                                                                       >  2:1310089/1‑90 (MQ=255)
                                                                        tGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGTCGCTGGACCAGCGCCa                                                                                                                                       <  1:1310089/93‑1 (MQ=255)
                                                                         gcgcGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACgg                                                                             <  1:242540/150‑1 (MQ=255)
                                                                          cgcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACggg                                                                            <  1:708815/150‑1 (MQ=255)
                                                                          cgcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACggg                                                                            <  1:1773/150‑1 (MQ=255)
                                                                          cgcgGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACggg                                                                            <  1:88452/150‑1 (MQ=255)
                                                                                  gTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGc                                                                    <  1:553308/150‑1 (MQ=255)
                                                                                   tCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATa                                                                                                                                <  2:460586/89‑1 (MQ=255)
                                                                                   tCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATa                                                                                                                                >  1:460586/1‑89 (MQ=255)
                                                                                    cACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGGCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTGTACTACGGTTTCCTGCAGGAACGGGTTGTCGGCGa                                                                  <  1:1771141/150‑1 (MQ=255)
                                                                                        cgAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATCTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACgg                                                              >  1:173288/1‑150 (MQ=255)
                                                                                         gagaCCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGAACGGCTCATGGCCGGCTGCAGGTTTAGATCTGGGTCGTGCAGGAACGGGTTGTCGGCGAACGGa                                                             >  2:194809/1‑150 (MQ=255)
                                                                                              cACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATg                                                                                                                   >  1:1074223/1‑91 (MQ=255)
                                                                                              cACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATg                                                                                                                   <  2:1074223/91‑1 (MQ=255)
                                                                                               aCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGt                                                                                                                                    <  2:1893183/73‑1 (MQ=255)
                                                                                               aCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGt                                                                                                                                    >  1:1893183/1‑73 (MQ=255)
                                                                                               aCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAg                                                          <  2:335731/147‑1 (MQ=255)
                                                                                               aCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAg                                                          >  1:335731/1‑147 (MQ=255)
                                                                                               aCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGcc                                                       <  1:2096405/150‑1 (MQ=255)
                                                                                               aCGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGcc                                                       >  2:107396/1‑150 (MQ=255)
                                                                                                   cTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCcgcac                                                   <  2:908421/150‑1 (MQ=255)
                                                                                                   cTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCcgcac                                                   <  2:1314411/150‑1 (MQ=255)
                                                                                                    ttCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCcgcacg                                                  >  1:181997/1‑150 (MQ=255)
                                                                                                    ttCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCcgcacg                                                  >  1:2079343/1‑150 (MQ=255)
                                                                                                       gAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCAt                                               >  1:124237/1‑150 (MQ=255)
                                                                                                          cAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGAc                                            <  2:1536484/150‑1 (MQ=255)
                                                                                                                     aCAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGt                                 >  2:1334386/1‑150 (MQ=255)
                                                                                                                         gTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCGAAACTGTCGGc                             >  1:822786/1‑150 (MQ=255)
                                                                                                                            gagaTAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCAc                                                                                               >  1:867376/1‑81 (MQ=255)
                                                                                                                            gagaTAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCAc                                                                                               <  2:867376/81‑1 (MQ=255)
                                                                                                                               aTAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGaa                       >  2:2073664/1‑150 (MQ=255)
                                                                                                                                 aGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAgg                     >  2:2074704/1‑150 (MQ=255)
                                                                                                                                  gTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTc                                                                       <  2:1659684/99‑1 (MQ=255)
                                                                                                                                  gTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTc                                                                       >  1:1659684/1‑99 (MQ=255)
                                                                                                                                   tCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTc                   >  1:2328535/1‑150 (MQ=255)
                                                                                                                                     gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACgg                                                                             >  2:378696/1‑90 (MQ=255)
                                                                                                                                     gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACgg                                                                             <  1:378696/90‑1 (MQ=255)
                                                                                                                                     gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCtt                 <  1:251201/150‑1 (MQ=255)
                                                                                                                                     gATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCCAAAAACTGTGGGCGAGGAATGTcgt                 >  2:1342942/1‑148 (MQ=255)
                                                                                                                                       tCATGTTTGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGTATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGa               <  1:929340/150‑1 (MQ=255)
                                                                                                                                        cATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGAt              <  2:2156931/150‑1 (MQ=255)
                                                                                                                                           ctattaccTGTCGCTGGACCAGCGCCAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGt            >  1:408083/7‑150 (MQ=255)
                                                                                                                                           ctattaccTGTCGCTGGACCAGCGCCAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGt            >  1:1692857/7‑150 (MQ=255)
                                                                                                                                            gTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCCGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTc           <  1:976929/150‑1 (MQ=255)
                                                                                                                                                gCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGc       <  2:124237/150‑1 (MQ=255)
                                                                                                                                                  cTGTCGCTGGACCAGCGCCAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTg     <  2:1692857/150‑1 (MQ=255)
                                                                                                                                                  cTGTCGCTGGACCAGCGCCAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTg     >  2:1961614/1‑150 (MQ=255)
                                                                                                                                                   tGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTgg    <  2:1929593/150‑1 (MQ=255)
                                                                                                                                                     gcgcTGGTCCAGCGACAGGTAATAGGAACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCACAAAGCTGTCGGCGACGAAGGTCTTGATGTCGAGCTGGTc  <  1:194809/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTCAGCCTCTGGCCAGTGCCGCCAGCACCAGTTCGCAGCTGCGCTGGCCATAGGCGGTGAGCACCTCGCGAATGCGCGGCAGGTCACGCGAGACCACGGCTTCGAGCAAGCGGGCGAACAGGTCGAGATAGTCGATCATGTTGGCCTGGCGCTGGTCCAGCGACAGGTAATAGGCACGGCTCATGGCCGGCTGCAGGTTTTCCACGGTTTCCTGCAGGAACGGGTTGTCGGCGAACGGATAGGCCGCACGCATGACCGCAAAACTGTCGGCGACGAAGGTCTTGATGTCGAGCTGGTC  >  NC_002947/4865920‑4866217

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: