Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,171 C→T 8.4% intergenic (+56/+67) 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,1710CT 8.4% 495.3 / 22.4 167intergenic (+56/+67)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (74/78);  new base T (5/9);  total (79/88)
Fisher's exact test for biased strand distribution p-value = 4.11e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 168 total aligned reads displayed.

GGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTA  >  NC_002947/4891023‑4891319
                                                                                                                                                    |                                                                                                                                                    
ggTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCt                                                                                                                                                     <  2:434770/150‑1 (MQ=255)
   aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                                  <  2:4654590/150‑1 (MQ=255)
   aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGACTGAACCATCTCGAAGGGCCATTCGCCGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                                  <  2:4019132/150‑1 (MQ=255)
    aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACTTGCCTGAAg                                                                                                                                                 <  2:4671292/150‑1 (MQ=255)
    aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAg                                                                                                                                                 <  2:2679239/150‑1 (MQ=255)
     ggCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAgg                                                                                                                                                <  2:3341765/150‑1 (MQ=255)
         aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                            <  1:4281539/150‑1 (MQ=255)
         aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                            >  2:3929935/1‑150 (MQ=255)
         aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                            <  2:4320770/150‑1 (MQ=255)
          aGCTGGCCGACAGTGTTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGATTGGTTGTGAGCCTGAGGCAGCCTGAAGCATCTCGTAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCACATGCGTGAAGGTtgtg                                                                                                                                           <  2:3824035/150‑1 (MQ=255)
          aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                           <  2:1514414/150‑1 (MQ=255)
          aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                           <  2:1666134/150‑1 (MQ=255)
          aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                           <  2:3790/150‑1 (MQ=255)
          aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                           <  2:3824073/150‑1 (MQ=255)
            cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                         <  1:5441249/150‑1 (MQ=255)
               gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTc                                                                                                                                      <  1:3189964/150‑1 (MQ=255)
               gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTc                                                                                                                                      >  2:4629592/1‑150 (MQ=255)
                ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                                     <  2:1394576/150‑1 (MQ=255)
                 cGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                                     <  2:3584271/149‑1 (MQ=255)
                  gACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCt                                                                                                                                   <  2:3895883/150‑1 (MQ=255)
                    cAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCtgt                                                                                                                                 <  2:3306616/150‑1 (MQ=255)
                      gTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTgg                                                                                                                               <  2:5104069/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)
                       tGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCGGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGACGGTTGTGGAGTCCCTGgggg                                                                                                                              >  2:5220624/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)
                          ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                           <  2:2320382/150‑1 (MQ=255)
                            aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                         <  1:3763251/150‑1 (MQ=18)
                            aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                         <  1:2308331/150‑1 (MQ=18)
                            aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCgg                                                                                                                         <  1:1662651/150‑1 (MQ=255)
                            aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGGTGTGGAGTCCCTGTGGGAGCgg                                                                                                                         >  2:5260095/1‑150 (MQ=255)
                            aTGACGCGGCGGCCGAGAACATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                         <  1:2940701/150‑1 (MQ=18)
                              gACGCGGCTGCCGAGACCATAATTAGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCTAATGGCCATTCGCTGGCAAGCCCTCTCCCACAGGATCAACACATGCCTGAAGGTTTTGGAGTCCCTGTGGGAGCGGGc                                                                                                                       <  2:1735419/150‑1 (MQ=255)
                              gACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGc                                                                                                                       <  2:2065517/150‑1 (MQ=17)
                                cgcgGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCtt                                                                                                                     <  1:4297999/150‑1 (MQ=17)
                                      ggcCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCcgc                                                                                                               <  2:133574/150‑1 (MQ=255)
                                          gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                         <  1:4162016/136‑1 (MQ=17)
                                          gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                         >  2:4162016/1‑136 (MQ=17)
                                          gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                         >  2:4161471/1‑136 (MQ=17)
                                                atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATGTCGATGGGCCATTCGCGGGCAAGCCCGGTCCCACAGGATCAACACATGCCTGAAGGTTGTGGAGTCCGTGTGGGAGCGGGCTTGCCCGCGAATGGGGCa                                                                                                     >  2:2899392/1‑150 (MQ=39)
                                                      cGTATGGTCTTGAAAGGTGGTGAGCCTGAGGCAGCCTGAACCACCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAg                                                                                               <  2:3538131/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)
                                                                  aaaGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGAtt                                                                                   >  1:1327723/1‑150 (MQ=255)
                                                                  aaaGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGAtt                                                                                   >  2:5557344/1‑150 (MQ=255)
                                                                   aaGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCGCCCACAGGATCACCACAAGCCGGAAGGTTGTGGAGTCCCGGTGGGCGCGGGCTTGCCCGCGAATGGGCCCGAACAGGCCATAGAGAttt                                                                                  >  1:2630657/1‑150 (MQ=38)
                                                                    aGGTTTTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACATGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                                 <  1:5601815/150‑1 (MQ=255)
                                                                    aGGTTGTGAGCCAGAGGAAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                                 <  2:4835300/150‑1 (MQ=255)
                                                                       ttGTGAGCCTGAGGCCGCCTGAACAATCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCCCCACATGCCTGAAGGTTGTGGAGCCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc                                                                              <  1:2707262/150‑1 (MQ=38)
                                                                                          tGACCCCTCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGcc                                                           <  2:5027479/150‑1 (MQ=255)
                                                                                                   tcGACGGGCCATACGCGGGCACGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGCCCCTGTGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcccc                                                  <  1:3469601/150‑1 (MQ=39)
                                                                                                     gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                                >  1:4127105/1‑150 (MQ=255)
                                                                                                     gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAGGCCTGAAGGTTGTGGAGTCCCTGTGGGAGTGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                                <  2:2924855/150‑1 (MQ=255)
                                                                                                         ggCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCAGTGGGAGCGGGCTTGCCCGCGAATGGGCCCGAACAGGCACTAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCaccacc                                            >  2:5247612/1‑150 (MQ=255)
                                                                                                             aTTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGCCAGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACACGCAATAGAGATTTAGGCCGCAACCCCCGCATCTGCCCGCAGCCCCaccaccacca                                        >  1:4426383/1‑146 (MQ=38)
                                                                                                                   gggCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAAAAGAGAATTATGCCGCTACCCCGGCATCTGACCGCAGCCCCACCTCCACGATCGCAc                                  >  1:4504149/1‑150 (MQ=39)
                                                                                                                        aGCCCGCTCCCCCTGGATCACCACAGGCCTGATGGTTGTGGTGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATc                             <  2:4698736/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)
                                                                                                                           ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGTCAATAGAGATTTAGGCCGCTACCCCGGCATCTTCCCGCAGCCCCACATCCTCGATCTCACTGATCGca                          >  1:4429635/1‑150 (MQ=255)
                                                                                                                           ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                          <  2:3149887/150‑1 (MQ=255)
                                                                                                                           ccGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGca                          >  1:5556559/1‑150 (MQ=255)
                                                                                                                             gCTCCCACCGTATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCgagga                        <  2:44697/150‑5 (MQ=255)
                                                                                                                             gCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTCCCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca                        >  2:3003929/1‑150 (MQ=255)
                                                                                                                             gCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGcaca                        >  1:4930921/1‑150 (MQ=255)
                                                                                                                               tCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACt                      >  1:4852984/1‑150 (MQ=255)
                                                                                                                                cccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTg                     <  2:4127105/150‑1 (MQ=255)
                                                                                                                                 ccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGc                    >  1:951113/1‑150 (MQ=255)
                                                                                                                                 ccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGc                    >  2:3727829/1‑150 (MQ=255)
                                                                                                                                  cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt                   >  1:1217515/1‑150 (MQ=255)
                                                                                                                                  cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt                   >  1:5317852/1‑150 (MQ=255)
                                                                                                                                   acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc                  >  2:3273153/1‑150 (MQ=255)
                                                                                                                                   acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc                  >  2:3123316/1‑150 (MQ=255)
                                                                                                                                   acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc                  >  1:1035668/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)
                                                                                                                                       gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg              >  1:4307345/1‑150 (MQ=255)
                                                                                                                                       gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg              >  2:5464207/1‑150 (MQ=255)
                                                                                                                                        aacaccacAACCCTCAAGGTTCTGCACACCCTGTGGGAGCGCCCTTCCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa             <  1:2811970/148‑1 (MQ=255)
                                                                                                                                            ccacATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATc         >  2:3853408/1‑150 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        >  1:4596401/1‑150 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        >  1:4623131/1‑150 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        >  1:1755346/1‑150 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        >  1:837815/1‑150 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        <  2:3737987/150‑1 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        >  1:3639106/1‑150 (MQ=255)
                                                                                                                                             cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc        >  1:5182036/1‑150 (MQ=255)
                                                                                                                                              acaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCg       >  2:1829554/1‑150 (MQ=255)
                                                                                                                                              acaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCg       >  2:3084485/1‑150 (MQ=255)
                                                                                                                                               caTGCCTGACGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGa      <  1:2691196/150‑1 (MQ=255)
                                                                                                                                               caTGCCTGAAGGTTGTGGAGTCCCTGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGa      <  2:1164027/150‑1 (MQ=255)
                                                                                                                                                aTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGAc     >  1:2234086/1‑150 (MQ=255)
                                                                                                                                                 tGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACg    >  1:1944197/1‑150 (MQ=255)
                                                                                                                                                   ccTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCACTAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTa  >  2:4917781/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGTAAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTA  >  NC_002947/4891023‑4891319

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: