Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,160 T→A 9.8% intergenic (+45/+78) PP_4295 → / ← PP_4296 TetR family transcriptional regulator/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,891,1600TA 9.8% 469.3 / 29.2 153intergenic (+45/+78)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (60/77);  new base A (5/10);  total (65/88)
Fisher's exact test for biased strand distribution p-value = 5.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.37e-01

Only 100 of 154 total aligned reads displayed.

CTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGAT  >  NC_002947/4891012‑4891309
                                                                                                                                                    |                                                                                                                                                     
cTGCGGTGACCGGTAAGTACATGCTGGCCGACAGTGATTATGACGCGGATGCAGAGAACATCATTCTTATGTTCTTTAATGGTTGTGATCCTGAGGCAGCATGTACCATCTAGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGTTc                                                                                                                                                      <  2:3670147/150‑1 (MQ=255)
   cGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATcacc                                                                                                                                                   <  1:4786995/150‑1 (MQ=255)
           ggTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCt                                                                                                                                           <  2:434770/150‑1 (MQ=255)
              aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                        <  2:4654590/150‑1 (MQ=255)
              aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGACTGAACCATCTCGAAGGGCCATTCGCCGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                        <  2:4019132/150‑1 (MQ=255)
               aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAg                                                                                                                                       <  1:4431761/150‑1 (MQ=255)
                ggCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAgg                                                                                                                                      <  2:3341765/150‑1 (MQ=255)
                  ccAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGtt                                                                                                                                    <  2:3872950/150‑1 (MQ=255)
                  ccAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGtt                                                                                                                                    <  1:1558258/150‑1 (MQ=255)
                    caGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGCATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCAACTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  <  2:1601790/149‑1 (MQ=255)
                    aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  >  2:3929935/1‑150 (MQ=255)
                    aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  <  2:4320770/150‑1 (MQ=255)
                    aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                  <  1:4281539/150‑1 (MQ=255)
                     aGCTGGCCGACAGTGTTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGATTGGTTGTGAGCCTGAGGCAGCCTGAAGCATCTCGTAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCACATGCGTGAAGGTtgtg                                                                                                                                 <  2:3824035/150‑1 (MQ=255)
                     aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:4593385/150‑1 (MQ=255)
                     aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:4657739/150‑1 (MQ=255)
                     aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:1666134/150‑1 (MQ=255)
                       cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                               <  1:508638/150‑1 (MQ=255)
                       cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                               <  1:5441249/150‑1 (MQ=255)
                        tGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAg                                                                                                                              <  1:4613534/150‑1 (MQ=255)
                          gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAATCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTc                                                                                                                            >  2:574924/1‑150 (MQ=255)
                          gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTc                                                                                                                            <  1:3189964/150‑1 (MQ=255)
                           ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                           <  1:69400/150‑1 (MQ=255)
                           ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                           <  2:1394576/150‑1 (MQ=255)
                            cGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                           <  2:3584271/149‑1 (MQ=255)
                            cGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                           >  1:3584271/1‑149 (MQ=255)
                               cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                                       <  2:3306616/150‑1 (MQ=255)
                                aGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgtg                                                                                                                      <  1:221975/150‑1 (MQ=255)
                                 gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  2:5104069/150‑1 (MQ=255)
                                 gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  1:1317116/150‑1 (MQ=255)
                                 gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCACTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  2:918417/150‑1 (MQ=255)
                                 gTGATTATGACGCGGCGGCCGAGACCATCAGTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                     <  2:4321479/150‑1 (MQ=255)
                                  tGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACCTGCCTGAAGGTTGTGGAGTCCCTGTggg                                                                                                                    >  2:5007229/1‑150 (MQ=255)
                                    aTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAg                                                                                                                  <  2:586635/150‑1 (MQ=18)
                                    aTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAg                                                                                                                  <  1:3841912/150‑1 (MQ=255)
                                     ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGc                                                                                                                 <  1:1494836/150‑1 (MQ=18)
                                     ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 <  2:4227205/134‑1 (MQ=255)
                                     ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                 >  1:4227205/1‑134 (MQ=255)
                                     ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                 <  2:2320382/150‑1 (MQ=255)
                                       aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:3763251/150‑1 (MQ=18)
                                       aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:2308331/150‑1 (MQ=18)
                                       aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCgg                                                                                                               <  1:1662750/150‑1 (MQ=255)
                                       aTGACGCGGCGGCCGAGAACATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:2940701/150‑1 (MQ=18)
                                         gACGCGGCTGCCGAGACCATAATTAGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCTAATGGCCATTCGCTGGCAAGCCCTCTCCCACAGGATCAACACATGCCTGAAGGTTTTGGAGTCCCTGTGGGAGCGGGc                                                                                                             <  2:1735419/150‑1 (MQ=255)
                                         gACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGc                                                                                                             <  2:2065517/150‑1 (MQ=17)
                                              ggcggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGcc                                                                                                        <  2:3813785/150‑1 (MQ=255)
                                                     gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               >  2:4161471/1‑136 (MQ=17)
                                                     gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                               <  1:4161471/136‑1 (MQ=17)
                                                     gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAAt                                                                                                 >  2:3016717/1‑150 (MQ=255)
                                                         ccATCATTCGTCTGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGc                                                                                             <  2:4157397/150‑1 (MQ=255)
                                                           atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATGTCGATGGGCCATTCGCGGGCAAGCCCGGTCCCACAGGATCAACACATGCCTGAAGGTTGTGGAGTCCGTGTGGGAGCGGGCTTGCCCGCGAATGGGGCa                                                                                           >  2:2899392/1‑150 (MQ=39)
                                                            tcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGTAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAg                                                                                          >  2:1398077/1‑150 (MQ=255)
                                                                 cGTATGGTCTTGAAAGGTGGTGAGCCTGAGGCAGCCTGAACCACCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAg                                                                                     <  2:3538131/150‑1 (MQ=255)
                                                                  gtttGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACCTGCCTGAAGGTTGTGGTGTCCCTGTGGGAGGGGGCTTGCCCGCGAATGGGCCAGAACAgg                                                                                    <  2:5162648/147‑1 (MQ=39)
                                                                     tGGTCTTGAAAGGTTGTGTGCCTGAGGCAGCCTGAACCATCTCTAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCCCTTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGGTTGCCCGCGAATGGGCCAGAACAGGCaa                                                                                 <  2:1965184/150‑1 (MQ=39)
                                                                       gTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATa                                                                               <  2:3148952/150‑1 (MQ=255)
                                                                        tCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATag                                                                              >  2:3679941/1‑150 (MQ=255)
                                                                         cTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGCCTGAAGGTTGTGGAGTCCCGGGGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATaga                                                                             >  2:2486102/1‑150 (MQ=255)
                                                                          ttGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACAGGATCACCACATGCCAGCAGGTTGTGGGGTGCCAGTGGGAGCGGGCGTGCCCGCGAAAGGGCCAGAACAGGCAAAagag                                                                            >  1:44697/1‑150 (MQ=37)
                                                                               aGGTTTTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACATGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                       <  1:5601815/150‑1 (MQ=255)
                                                                                  ttGTGAGCCTGAGGCCGCCTGAACAATCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCCCCACATGCCTGAAGGTTGTGGAGCCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc                                                                    <  1:2707262/150‑1 (MQ=38)
                                                                                  tggtgCGCCTGCGGCAGCCTGAACCATCTCGACGGGCCATGCGCGGGCACGCCCGCGCCCCCCGGATCCCCACCTGCCTGACGGGTGTGGCGTCCCTGTGGGAGCGGGCTGGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc                                                                    <  2:1135695/148‑1 (MQ=255)
                                                                                      gAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGCCTGAAGGTTGTGGAGTCCCTATGGGAGCGGGCTTTCCCGCGAAAGGGCCAGAACACGCAATAGAGATTTAGGCCGcg                                                                >  1:4813710/1‑149 (MQ=39)
                                                                                                     tGACCCCTCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGcc                                                 <  2:5027479/150‑1 (MQ=255)
                                                                                                                gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                      >  1:4127105/1‑150 (MQ=255)
                                                                                                                gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAGGCCTGAAGGTTGTGGAGTCCCTGTGGGAGTGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                      <  2:2924855/150‑1 (MQ=255)
                                                                                                                    ggCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAcctcc                                  >  2:5248050/1‑150 (MQ=255)
                                                                                                                    ggCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCAGTGGGAGCGGGCTTGCCCGCGAATGGGCCCGAACAGGCACTAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCaccacc                                  >  2:5247612/1‑150 (MQ=255)
                                                                                                                      ccATTCGCGGGCAAGCCCGCTCCCACAGGGACGACACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACctcctc                                >  2:3785388/1‑150 (MQ=17)
                                                                                                                        aTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGCCAGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACACGCAATAGAGATTTAGGCCGCAACCCCCGCATCTGCCCGCAGCCCCaccaccacca                              >  1:4426383/1‑146 (MQ=38)
                                                                                                                              gggCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAAAAGAGAATTATGCCGCTACCCCGGCATCTGACCGCAGCCCCACCTCCACGATCGCAc                        >  1:4504149/1‑150 (MQ=39)
                                                                                                                                 cAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGa                     <  2:1313308/150‑1 (MQ=255)
                                                                                                                                   aGCCCGCTCCCCCTGGATCACCACAGGCCTGATGGTTGTGGTGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATc                   <  2:4698736/150‑1 (MQ=255)
                                                                                                                                   aGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATc                   <  2:4674747/150‑1 (MQ=255)
                                                                                                                                   aGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATc                   <  1:3679941/150‑1 (MQ=255)
                                                                                                                                    gCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCg                  >  1:3395685/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGc                 >  1:1833223/1‑150 (MQ=255)
                                                                                                                                     cccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGc                 >  1:1313007/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                >  1:2766609/1‑150 (MQ=18)
                                                                                                                                      ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                >  1:5556559/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                >  1:331340/1‑150 (MQ=255)
                                                                                                                                      ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                <  2:3149887/150‑1 (MQ=255)
                                                                                                                                        gCTCCCACCGTATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCgagga              <  2:44697/150‑5 (MQ=255)
                                                                                                                                        gCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTCCCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca              >  2:3003929/1‑150 (MQ=255)
                                                                                                                                        gCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca              >  1:4930921/1‑150 (MQ=255)
                                                                                                                                           cccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTg           <  2:4127105/150‑1 (MQ=255)
                                                                                                                                            ccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGc          >  2:3727829/1‑150 (MQ=255)
                                                                                                                                             cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt         >  1:5317852/1‑150 (MQ=255)
                                                                                                                                             cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt         >  1:1217515/1‑150 (MQ=255)
                                                                                                                                              acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  2:5327437/1‑150 (MQ=255)
                                                                                                                                              acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  1:1035668/1‑150 (MQ=255)
                                                                                                                                              acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc        >  2:3273153/1‑150 (MQ=255)
                                                                                                                                               caGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcg       >  1:5355516/1‑150 (MQ=255)
                                                                                                                                                aGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgt      >  1:3044566/1‑150 (MQ=255)
                                                                                                                                                 ggATCACCACTTGCCTGACGGTTGTGGCGGCCCGGGGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtc     <  1:609575/150‑1 (MQ=255)
                                                                                                                                                  gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCTCACTGCtcgtcg    >  2:4864851/1‑150 (MQ=255)
                                                                                                                                                  gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg    >  1:4307196/1‑150 (MQ=255)
                                                                                                                                                  gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg    >  2:5464207/1‑150 (MQ=255)
                                                                                                                                                   aTCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa   >  2:3825092/1‑150 (MQ=255)
                                                                                                                                                    tCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGat  >  1:4335102/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CTGCGGTGACCGGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGAT  >  NC_002947/4891012‑4891309

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