Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,986,810 G→T 6.0% D96E (GAC→GAA ycaC‑II ← putative hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,986,8100GT 6.0% 709.1 / 24.4 215D96E (GAC→GAAycaC‑IIputative hydrolase
Reads supporting (aligned to +/- strand):  ref base G (108/93);  new base T (7/6);  total (115/100)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 216 total aligned reads displayed.

CGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAG  >  NC_002947/5986662‑5986959
                                                                                                                                                    |                                                                                                                                                     
 gTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTc                                                                                                                                                     >  1:3881301/1‑150 (MQ=255)
    gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTc                                                                                                                                                  <  2:4342026/150‑1 (MQ=255)
    gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTc                                                                                                                                                  <  2:1364768/150‑1 (MQ=255)
    gAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTc                                                                                                                                                  <  1:135958/150‑1 (MQ=255)
        gTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTg                                                                                                                                              >  1:3859422/1‑150 (MQ=255)
           cccGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTcc                                                                                                                                           >  1:800770/1‑150 (MQ=255)
            ccGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTccc                                                                                                                                          >  1:494540/1‑150 (MQ=255)
              gAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAg                                                                                                                                        <  2:4933437/150‑1 (MQ=255)
               aGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAgg                                                                                                                                       >  2:1141955/1‑150 (MQ=255)
                   aTCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGtt                                                                                                                                   <  2:5172131/150‑1 (MQ=255)
                     cGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGa                                                                                                                                 <  1:3109682/150‑1 (MQ=255)
                      ggTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGAt                                                                                                                                >  1:3223587/1‑150 (MQ=255)
                         caccacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTg                                                                                                                             <  2:2837302/150‑1 (MQ=255)
                            cacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGcc                                                                                                                          >  2:3505222/1‑150 (MQ=255)
                            cacGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGcc                                                                                                                          >  1:2553945/1‑150 (MQ=255)
                                    cTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACggg                                                                                                                  <  1:2846689/150‑1 (MQ=255)
                                     ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGc                                                                                                                 >  2:1310776/1‑150 (MQ=255)
                                     ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGc                                                                                                                 >  1:2947662/1‑150 (MQ=255)
                                     ttCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGc                                                                                                                 >  1:2680471/1‑150 (MQ=255)
                                      tCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCa                                                                                                                <  1:1169504/150‑1 (MQ=255)
                                        aaaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAATTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAAt                                                                                                              >  1:5436074/1‑150 (MQ=255)
                                         aaaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATg                                                                                                             >  1:5247927/1‑150 (MQ=255)
                                          aaCCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGt                                                                                                            >  1:1066887/1‑150 (MQ=255)
                                            ccTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTAc                                                                                                          >  1:134835/1‑150 (MQ=255)
                                              ttCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACgg                                                                                                        >  2:2256594/1‑150 (MQ=255)
                                                   ccAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCAt                                                                                                   >  1:2425263/1‑150 (MQ=255)
                                                      ttcccGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCATTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCg                                                                                                >  1:1132674/4‑150 (MQ=255)
                                                      gcgcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCg                                                                                                <  2:1872624/150‑1 (MQ=255)
                                                       cgcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCgg                                                                                               <  2:3560261/150‑1 (MQ=255)
                                                        gcCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCggg                                                                                              <  1:3974005/150‑1 (MQ=255)
                                                         ccGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGa                                                                                             <  2:2553945/150‑1 (MQ=255)
                                                          cGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGaa                                                                                            <  1:4551965/150‑1 (MQ=255)
                                                          cGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGaa                                                                                            <  1:2215330/150‑1 (MQ=255)
                                                           gACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAAc                                                                                           >  2:215592/1‑150 (MQ=255)
                                                                ggTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACAtctc                                                                                      >  2:701785/1‑150 (MQ=255)
                                                                 gTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATctct                                                                                     >  1:556841/1‑150 (MQ=255)
                                                                  tCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCtt                                                                                    <  2:2402619/150‑1 (MQ=255)
                                                                     ggAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCt                                                                                                                                                    >  2:2218431/1‑83 (MQ=255)
                                                                      gAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAg                                                                                >  1:2416132/1‑150 (MQ=255)
                                                                           ccACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCgg                                                                           <  2:816164/150‑1 (MQ=255)
                                                                           ccACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCgg                                                                           >  2:334945/1‑150 (MQ=255)
                                                                             aCGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGgcgc                                                                                                     <  1:1439267/122‑1 (MQ=255)
                                                                             aCGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGgcgc                                                                                                     >  2:1439267/1‑122 (MQ=255)
                                                                               gCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCAc                                                                       <  1:265148/150‑1 (MQ=255)
                                                                               gCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCAc                                                                       <  2:755475/150‑1 (MQ=255)
                                                                                cacacaTCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCAcc                                                                      <  2:556841/150‑1 (MQ=255)
                                                                                       cGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGc                                                               >  1:2079222/1‑150 (MQ=255)
                                                                                         gTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCa                                                             >  1:1904466/1‑150 (MQ=255)
                                                                                            accacACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCgg                                                                                               <  1:171194/113‑1 (MQ=255)
                                                                                             ccacACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCATCAGCGGGCCATTCg                                                         >  1:4436592/1‑150 (MQ=255)
                                                                                             ccacACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCg                                                         >  1:2464487/1‑150 (MQ=255)
                                                                                                 acCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTg                                                                                                                                  <  2:1717637/73‑1 (MQ=255)
                                                                                                   cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAg                                                                                                                                        <  2:3572604/65‑1 (MQ=255)
                                                                                                   cGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAg                                                                                                                                        >  1:3572604/1‑65 (MQ=255)
                                                                                                     gCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGt                                                                                                                                                >  1:1448015/1‑55 (MQ=255)
                                                                                                     gCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGc                                                 >  2:4685383/1‑150 (MQ=255)
                                                                                                       gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                               <  2:1033172/150‑1 (MQ=255)
                                                                                                       gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                               >  1:730852/1‑150 (MQ=255)
                                                                                                       gatgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                               <  2:81144/150‑1 (MQ=255)
                                                                                                       gatgatCAGTTGCTTGCGGCCCGTGGCCTTGAGTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTc                                               <  2:2997968/150‑1 (MQ=255)
                                                                                                         tgatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGa                                             >  1:1664602/1‑150 (MQ=255)
                                                                                                          gatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                            >  2:5036570/1‑150 (MQ=255)
                                                                                                          gatCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGaa                                            >  2:2398533/1‑150 (MQ=255)
                                                                                                             cAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCggg                                                       >  1:1238435/1‑136 (MQ=255)
                                                                                                             cAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCggg                                                       <  2:1238435/136‑1 (MQ=255)
                                                                                                                ttgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                 <  1:2257499/107‑1 (MQ=255)
                                                                                                                ttgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCa                                                                                 >  2:2257499/1‑107 (MQ=255)
                                                                                                                 tgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                     <  2:1447211/150‑1 (MQ=255)
                                                                                                                 tgcttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                     <  1:3099415/150‑1 (MQ=255)
                                                                                                                 tgcttgcGGCCCGTGGCCGTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCtggtg                                     <  2:3924575/150‑1 (MQ=255)
                                                                                                                    ttgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGgtggtg                                  >  1:3300048/1‑150 (MQ=255)
                                                                                                                     tgcGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGa                                 <  1:1028736/150‑1 (MQ=255)
                                                                                                                     tcttgCCGGTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTggtggtca                                 <  2:1677880/146‑3 (MQ=2)
                                                                                                                        ggCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGa                              >  2:3260243/1‑150 (MQ=255)
                                                                                                                           cggTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAgaacag                           >  1:233465/3‑145 (MQ=2)
                                                                                                                           ccGTGGCCTTGATTGCCTTGACGAAGTGTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGGTCGAAGCTGGTGGTGAGGATgg                           >  2:519439/1‑150 (MQ=255)
                                                                                                                           ccGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCgg                                                                                                                        >  2:4187983/1‑57 (MQ=255)
                                                                                                                           ccGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCgg                                                                                                                        <  1:4187983/57‑1 (MQ=255)
                                                                                                                            ggTGGCCTTCACTGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAgaacagt                          <  1:988327/149‑7 (MQ=2)
                                                                                                                            cGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGt                          <  1:1568358/150‑1 (MQ=255)
                                                                                                                             gTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTc                         <  2:2707846/150‑1 (MQ=255)
                                                                                                                                  cTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTATGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                    >  2:3358551/1‑150 (MQ=255)
                                                                                                                                  cTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                    >  2:2869980/1‑150 (MQ=255)
                                                                                                                                  cTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAg                    >  2:1195928/1‑150 (MQ=255)
                                                                                                                                    tGATTGCCTTGACGAAGTCTACGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGc                  <  2:1530726/150‑1 (MQ=255)
                                                                                                                                     gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc                 >  1:3366681/1‑150 (MQ=255)
                                                                                                                                     gATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGcc                 >  1:4461532/1‑150 (MQ=255)
                                                                                                                                      aTTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCa                <  2:5323263/150‑1 (MQ=255)
                                                                                                                                       ttGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaa               >  1:4654211/1‑150 (MQ=255)
                                                                                                                                       ctGCCTTGACGAATTCCTCATTGTCCCAGGCGTTGATCTGGCCGGGGCGAGCGATATAAGGCGCATCCGGGAATGTTTCTTTCAGCTCCGGCACCAATGGGCCATTCGGGCCGTTTTCAAAACTGGTGGTCAGAACAGTCGGCAgcttga               >  2:3979398/2‑145 (MQ=2)
                                                                                                                                        tGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCaaa              <  2:2947662/150‑1 (MQ=255)
                                                                                                                                          ccTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGa            >  1:2285181/1‑150 (MQ=255)
                                                                                                                                           cTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATCCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGaa           >  2:4787440/1‑150 (MQ=255)
                                                                                                                                               aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg       >  1:3657046/1‑150 (MQ=255)
                                                                                                                                               aCGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTg       <  2:4428936/150‑1 (MQ=255)
                                                                                                                                                 gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc     <  1:690387/150‑1 (MQ=255)
                                                                                                                                                 gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc     >  2:916448/1‑150 (MQ=255)
                                                                                                                                                 gAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGc     >  1:484980/1‑150 (MQ=255)
                                                                                                                                                  aaGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGcc    <  1:1824828/150‑1 (MQ=255)
                                                                                                                                                    gTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAg  >  2:319881/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CGTTGAAGGTGCCCGAGGCATCGGTCACCACGAACACTTCAAAACCTTCCGCCAGCGCCGACAGGGTCGGGAAGGCCACGCACACATCGGTCACCACACCGGCGATGATCAGTTGCTTGCGGCCCGTGGCCTTGATTGCCTTGACGAAGTCTTCGTTGTCCCAGGCGTTGATCTGGCCGGGACGGGCAATGTACGGCGCATCCGGGAACATCTCTTTCAGCTCGGGCACCAGCGGGCCATTCGGGCCTTGCTCGAAGCTGGTGGTGAGGATGGTCGGCAGGCCAAAGAACTTGGCCAG  >  NC_002947/5986662‑5986959

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