Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,891,173 G→C 9.2% intergenic (+58/+65) 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,1730GC 9.2% 536.0 / 30.8 168intergenic (+58/+65)PP_4295/PP_4296TetR family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (75/77);  new base C (6/10);  total (81/87)
Fisher's exact test for biased strand distribution p-value = 4.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 169 total aligned reads displayed.

AAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCG  >  NC_002947/4891026‑4891322
                                                                                                                                                   |                                                                                                                                                     
aaGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGaa                                                                                                                                                     <  2:4654590/150‑1 (MQ=255)
 aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACTTGCCTGAAg                                                                                                                                                    <  2:4671292/150‑1 (MQ=255)
 aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAg                                                                                                                                                    <  2:2679239/150‑1 (MQ=255)
 aGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAg                                                                                                                                                    <  1:4431761/150‑1 (MQ=255)
  ggCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAgg                                                                                                                                                   <  2:3341765/150‑1 (MQ=255)
   gCCAAGCTGGCAGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGGATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACAATCTCGAAGGGACATTCGCGGGCAAGCCCGCTCCACCAGCATCACCAAATGCCTCACGGt                                                                                                                                                  >  2:2938353/1‑150 (MQ=255)
    ccAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGtt                                                                                                                                                 <  2:3872950/150‑1 (MQ=255)
      caGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGCATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCAACTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                               <  2:1601790/149‑1 (MQ=255)
      aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                               <  2:4320770/150‑1 (MQ=255)
      aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                               <  1:4281539/150‑1 (MQ=255)
      aaGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgt                                                                                                                                               >  2:3929935/1‑150 (MQ=255)
       aGCTGGCCGACAGTGTTTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGATTGGTTGTGAGCCTGAGGCAGCCTGAAGCATCTCGTAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGTTCTCCACATGCGTGAAGGTtgtg                                                                                                                                              <  2:3824035/150‑1 (MQ=255)
       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                              <  2:3328018/150‑1 (MQ=255)
       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                              <  2:4593385/150‑1 (MQ=255)
       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                              <  2:1666174/150‑1 (MQ=255)
       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                              <  2:1514414/150‑1 (MQ=255)
       aGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                              <  2:3790/150‑1 (MQ=255)
        gCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCGGAAGGTTGTgg                                                                                                                                             >  1:390561/1‑150 (MQ=255)
         cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                            <  1:5441252/150‑1 (MQ=255)
         cTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGa                                                                                                                                            <  1:508638/150‑1 (MQ=255)
          tGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAg                                                                                                                                           <  1:4613534/150‑1 (MQ=255)
            gCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTc                                                                                                                                         <  1:3189964/150‑1 (MQ=255)
             ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                                        <  2:1394576/150‑1 (MQ=255)
             ccGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                                        <  1:69400/150‑1 (MQ=255)
              cGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTcc                                                                                                                                        >  1:3584271/1‑149 (MQ=255)
               gACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCt                                                                                                                                      <  2:3895883/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)
                       ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGc                                                                                                                              <  1:1494836/150‑1 (MQ=18)
                       ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTtgtg                                                                                                                                              <  2:4227205/134‑1 (MQ=255)
                       ttATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGc                                                                                                                              <  2:2320382/150‑1 (MQ=255)
                         aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                            <  1:2308331/150‑1 (MQ=18)
                         aTGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                            <  1:3763251/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)
                           gACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGc                                                                                                                          <  2:2065517/150‑1 (MQ=17)
                             cgcgGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGCCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCtt                                                                                                                        <  1:4297999/150‑1 (MQ=17)
                                       gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                            <  1:4162016/136‑1 (MQ=17)
                                       gagaCCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCgg                                                                                                                            >  2:4161471/1‑136 (MQ=17)
                                             atcatTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATGTCGATGGGCCATTCGCGGGCAAGCCCGGTCCCACAGGATCAACACATGCCTGAAGGTTGTGGAGTCCGTGTGGGAGCGGGCTTGCCCGCGAATGGGGCa                                                                                                        >  2:2899392/1‑150 (MQ=39)
                                                   cGTATGGTCTTGAAAGGTGGTGAGCCTGAGGCAGCCTGAACCACCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAg                                                                                                  <  2:3538131/150‑1 (MQ=255)
                                                         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)
                                                               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)
                                                                 aGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                                    <  1:5601818/150‑1 (MQ=255)
                                                                 aGGTTGTGAGCCAGAGGAAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTa                                                                                    <  2:4835300/150‑1 (MQ=255)
                                                                    ttGTGAGCCTGAGGCCGCCTGAACAATCTCGAAGGGCCATTCGCGGGCACGCCCGCTCCCACAGGATCCCCACATGCCTGAAGGTTGTGGAGCCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGc                                                                                 <  1:2707262/150‑1 (MQ=38)
                                                                        gAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACAAGCCTGAAGGTTGTGGAGTCCCTATGGGAGCGGGCTTTCCCGCGAAAGGGCCAGAACACGCAATAGAGATTTAGGCCGcg                                                                             >  1:4813710/1‑149 (MQ=39)
                                                                                       tGACCCCTCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACCGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGcc                                                              <  2:5027479/150‑1 (MQ=255)
                                                                                              tctcGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCCCCACATGCCTGTAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACCGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcc                                                       >  2:3589163/1‑150 (MQ=255)
                                                                                                ccGACGGGCCTTGCGCGGGCACGCCCGCCCCCACAGGACCCCCACAGCCCTGACGGTGGTGGCGCCCCGGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGcccc                                                     <  1:5260095/149‑1 (MQ=255)
                                                                                                  gAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCAc                                                   >  1:4127105/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)
                                                                                                                    aaGCCCGCTCCCACAGGGACTCCACAACCTTCAGGCATGTGGTGATCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGAt                                 >  2:4036908/1‑150 (MQ=18)
                                                                                                                     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)
                                                                                                                        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)
                                                                                                                              ccACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGc                       >  1:951113/1‑150 (MQ=255)
                                                                                                                               cacaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCt                      >  1:1217515/1‑150 (MQ=255)
                                                                                                                                acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc                     >  2:3273153/1‑150 (MQ=255)
                                                                                                                                acaGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtc                     >  1:1035668/1‑150 (MQ=255)
                                                                                                                                 ccggATCACCACCTGCCGGACGGTTGGGGCGTCCCTGTGGGCGGGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcg                    <  1:5368190/148‑1 (MQ=255)
                                                                                                                                 caGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcg                    >  1:5355516/1‑150 (MQ=255)
                                                                                                                                    gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg                 >  2:5464207/1‑150 (MQ=255)
                                                                                                                                    gATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCtcgtcg                 >  1:4307345/1‑150 (MQ=255)
                                                                                                                                     aacaccacAACCCTCAAGGTTCTGCACACCCTGTGGGAGCGCCCTTCCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa                <  1:2811970/148‑1 (MQ=255)
                                                                                                                                     aTCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGa                >  2:3825092/1‑150 (MQ=255)
                                                                                                                                          cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc           >  1:1755346/1‑150 (MQ=255)
                                                                                                                                          cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc           >  1:3639106/1‑150 (MQ=255)
                                                                                                                                          cacaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATcc           >  1:4623131/1‑150 (MQ=255)
                                                                                                                                           acaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCg          >  2:3084485/1‑150 (MQ=255)
                                                                                                                                           acaTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCg          >  2:1829554/1‑150 (MQ=255)
                                                                                                                                            caTGCCTGACGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGa         <  1:2691196/150‑1 (MQ=255)
                                                                                                                                            caTGCCTGAAGGTTGTGGAGTCCCTGTGGGCGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGa         <  2:1164027/150‑1 (MQ=255)
                                                                                                                                            caTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGa         >  1:3051987/1‑150 (MQ=255)
                                                                                                                                             aTGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGAc        >  1:2234086/1‑150 (MQ=255)
                                                                                                                                              tGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACg       >  1:1944197/1‑150 (MQ=255)
                                                                                                                                                ccTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCACTAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTa     >  2:4917781/1‑150 (MQ=255)
                                                                                                                                                  tGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATc   >  2:1936444/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
AAGGCCAAGCTGGCCGACAGTGATTATGACGCGGCGGCCGAGACCATCATTCGTATGGTCTTGAAAGGTTGTGAGCCTGAGGCAGCCTGAACCATCTCGAAGGGCCATTCGCGGGCAAGCCCGCTCCCACAGGATCACCACATGCCTGAAGGTTGTGGAGTCCCTGTGGGAGCGGGCTTGCCCGCGAATGGGCCAGAACAGGCAATAGAGATTTAGGCCGCTACCCCGGCATCTGCCCGCAGCCCCACCTCCTCGATCGCACTGATCGCACACTGCTCGTCGATATCCGACGTATCG  >  NC_002947/4891026‑4891322

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: