Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,584,942 A→C 6.8% intergenic (+46/‑89) benC → / → benD benzoate 1,2‑dioxygenase electron transfer component/1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,584,9420AC 6.8% 341.0 / 13.7 127intergenic (+46/‑89)benC/benDbenzoate 1,2‑dioxygenase electron transfer component/1,6‑dihydroxycyclohexa‑2,4‑diene‑1‑carboxylate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (60/57);  new base C (4/5);  total (65/62)
Fisher's exact test for biased strand distribution p-value = 7.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.32e-01

Only 100 of 129 total aligned reads displayed.

TGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCG  >  NC_002947/3584795‑3585088
                                                                                                                                                    |                                                                                                                                                  
tgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGACTTACCACGTCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGc                                                                                                                                                  >  2:1189459/1‑150 (MQ=255)
tgCAGCCCGCCACGCATGGTTGTATCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGGTCCCTGATTCAAAGCAGTACGCCGACCCTGTGGGAGCGGGCGAGc                                                                                                                                                  <  2:195345/150‑1 (MQ=255)
tccGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGc                                                                                                                                                  >  2:1654228/3‑150 (MQ=255)
    gCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGAGCGTGAGGAGGGCATTACCCCGGCCAACTTAGATGACGAGAAGTTTCCGGCAGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGc                                                                                                                                              >  1:2520241/1‑150 (MQ=255)
      ccGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGCGCGACCGACGCCGCTc                                                                                                                                            >  1:2284579/1‑150 (MQ=255)
      ccGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTc                                                                                                                                            >  1:1662700/1‑150 (MQ=255)
       cGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTATTACGTGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc                                                                                                                                           <  1:53903/150‑1 (MQ=255)
       cGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc                                                                                                                                           <  2:735094/150‑1 (MQ=255)
       cGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc                                                                                                                                           <  2:1662700/150‑1 (MQ=255)
              cATGGTGGCAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTGGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                                    >  2:1446059/1‑150 (MQ=255)
              cATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                                    >  2:2607873/1‑150 (MQ=255)
              cATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGCTGTCACCCTGACGCTCCCTGATTCAAAGCACTACGACGGCCCTGTCTGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                                    >  1:1301694/1‑150 (MQ=255)
               aTGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGa                                                                                                                                   >  1:2757029/1‑150 (MQ=255)
                           gTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCtt                                                                                                                       >  2:1945608/1‑150 (MQ=255)
                             cAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTc                                                                                                                     <  1:1127307/150‑1 (MQ=255)
                                  aGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGACTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCt                                                                                                                >  2:1097007/1‑150 (MQ=255)
                                    tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTAACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGACCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                              <  1:1218376/150‑1 (MQ=255)
                                       gtgcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAg                                                                                                           >  2:1970727/1‑150 (MQ=255)
                                       gtgcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAg                                                                                                           >  1:2210715/1‑150 (MQ=255)
                                        tgcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGc                                                                                                          >  1:2170832/1‑150 (MQ=255)
                                         tcgtgAGCAGGGCATTACCCCGGCCAAGCTCTACCACGCGACGCTTCCGACAGCGGTCTCACGCTTCCTGATTCAAAGCAATACGCCGACACTCTGGCAGCGGCCGAACCCACCCCCACAGCGACCTCGCCGCCTTTCAGGCTTTTAGCa                                                                                                         >  2:2545561/2‑150 (MQ=255)
                                           gtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAag                                                                                                       <  1:1945608/150‑1 (MQ=255)
                                            tgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAaga                                                                                                      >  1:121142/1‑150 (MQ=255)
                                                cAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGt                                                                                                  >  2:2677585/1‑150 (MQ=255)
                                                cAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGt                                                                                                  >  2:1535703/1‑150 (MQ=255)
                                                 aGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                                    <  1:975477/131‑1 (MQ=255)
                                                 aGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCtt                                                                                                               <  2:2395455/136‑1 (MQ=255)
                                                   ggCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAagagag                                                                                                   >  1:349469/1‑146 (MQ=255)
                                                   ggCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAagagag                                                                                                   <  2:349469/146‑1 (MQ=255)
                                                            ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  2:1045374/150‑1 (MQ=255)
                                                            ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  1:1446059/150‑1 (MQ=255)
                                                             cGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAggg                                                                                     >  2:3070335/1‑150 (MQ=255)
                                                                ccAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                                    <  2:520150/116‑1 (MQ=255)
                                                                 cAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctc                                                                                 >  2:2814969/1‑150 (MQ=255)
                                                                   aCTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctctc                                                                               <  2:121142/150‑1 (MQ=255)
                                                                    cTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGa                                                                                                                                                    <  1:3039837/80‑1 (MQ=255)
                                                                    cTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGa                                                                                                                                                    >  2:3039837/1‑80 (MQ=255)
                                                                    cTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTTCCACAGGGTCTCTcg                                                                              >  2:1185027/1‑150 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTTCTGCCACAGGGTCTCTcgc                                                                             >  2:2681029/1‑150 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  1:1741093/141‑1 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      >  2:1741093/1‑141 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             <  2:490533/150‑1 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             <  2:1841232/150‑1 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGCGGGAGGGGCCGGCCCCGCCCCAGCTGCGACCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGACTCgcgc                                                                             >  1:1752713/1‑150 (MQ=255)
                                                                        tactacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGaa                                                                          >  1:348728/1‑150 (MQ=255)
                                                                        tactacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGaa                                                                          <  2:348728/150‑1 (MQ=255)
                                                                          ctacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACt                                                                        >  1:1968623/1‑150 (MQ=255)
                                                                           tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGCGAAAGCTCCCACAGTGATCTCGCATCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                       <  1:2677585/150‑1 (MQ=255)
                                                                           tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGACTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                       <  1:1097007/150‑1 (MQ=255)
                                                                              gagaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGgag                                                                    >  1:101077/1‑150 (MQ=255)
                                                                               agaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAg                                                                                                                   <  1:2148125/102‑1 (MQ=255)
                                                                               agaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAg                                                                                                                   >  2:2148125/1‑102 (MQ=255)
                                                                               agaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGgaga                                                                   <  1:770093/150‑1 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATATGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              >  1:2919467/1‑150 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCGGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGATTTGAGCAAGAGAGTTGCTGCGACAGGGCCTCTCGCGGACTGGAGGTTTAc                                                              >  2:2224277/1‑150 (MQ=39)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              <  1:1003959/150‑1 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              >  1:480719/1‑150 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGATCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGTGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              >  2:1118158/1‑150 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCAGCTGCGACAGTGATCGCGCAACCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              >  1:1752989/1‑150 (MQ=255)
                                                                                     ttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACa                                                             <  2:2757029/150‑1 (MQ=255)
                                                                                     ttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACa                                                             <  1:1535703/150‑1 (MQ=255)
                                                                                     tggcggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCGCCCACAGTGATCCCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACa                                                             <  2:1448291/148‑1 (MQ=255)
                                                                                        cggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCATGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACatat                                                          <  2:1266385/150‑1 (MQ=255)
                                                                                          gctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGa                                                        >  1:966154/1‑150 (MQ=255)
                                                                                           ctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATGTGac                                                       >  2:1446923/1‑149 (MQ=255)
                                                                                           ctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACACATGaa                                                       <  2:2138298/150‑1 (MQ=255)
                                                                                           ctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTGTCATATGaa                                                       <  1:1189459/150‑1 (MQ=255)
                                                                                            tgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACATTGATCGCGCAGCCTTTCAGGCTTTGAGCAATAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaac                                                      <  2:2520241/150‑1 (MQ=255)
                                                                                                ggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaacaac                                                   <  2:1968623/150‑1 (MQ=255)
                                                                                                      aCGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTc                                             <  1:3002481/150‑1 (MQ=255)
                                                                                                        gCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCca                                           >  1:2723614/1‑150 (MQ=255)
                                                                                                             cTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAag                                                                                                       <  2:2310822/85‑1 (MQ=255)
                                                                                                             cTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAag                                                                                                       >  1:2310822/1‑85 (MQ=255)
                                                                                                              tGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggcaa                                     >  2:638826/1‑150 (MQ=255)
                                                                                                                aTTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggcaagg                                   <  1:2607873/150‑1 (MQ=255)
                                                                                                                 ttCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGt                                  >  1:2808475/1‑150 (MQ=255)
                                                                                                                  tCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGtt                                 <  2:1752713/150‑1 (MQ=255)
                                                                                                                    aaaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGc                                                                                                          <  2:1283338/75‑1 (MQ=255)
                                                                                                                    aaaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGc                               >  2:125291/1‑150 (MQ=255)
                                                                                                                     aaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGcc                              >  1:1432993/1‑150 (MQ=255)
                                                                                                                     aaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGcc                              <  1:1654228/150‑1 (MQ=255)
                                                                                                                      aGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaa                                                       <  1:552240/124‑1 (MQ=255)
                                                                                                                      aGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaa                                                       >  2:552240/1‑124 (MQ=255)
                                                                                                                       gCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCt                            <  2:101077/150‑1 (MQ=255)
                                                                                                                             cGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCAc                      <  2:2170832/150‑1 (MQ=255)
                                                                                                                                  aCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTg                 <  2:2723614/150‑1 (MQ=255)
                                                                                                                                   cccTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgc                <  2:2992966/150‑1 (MQ=255)
                                                                                                                                    ccTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTg                           >  2:3062087/1‑138 (MQ=255)
                                                                                                                                      tgtgGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgc             >  2:1287002/1‑150 (MQ=255)
                                                                                                                                      tgtgGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgc             <  1:1970727/150‑1 (MQ=255)
                                                                                                                                      tgtgGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGCCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGACCCAAGGAGTTACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgcagc             <  1:2383510/150‑1 (MQ=255)
                                                                                                                                          ggAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGaa                                                                          >  1:441608/1‑85 (MQ=255)
                                                                                                                                          ggAGCGGGCGAGCCCGCTCCCACAGTGATCGCGTAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGCCACCGGTGCCGCCCaa         <  2:2284579/150‑1 (MQ=255)
                                                                                                                                          ggAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCaa         >  1:2140891/1‑150 (MQ=255)
                                                                                                                                           gAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAg        <  2:1432993/150‑1 (MQ=255)
                                                                                                                                            aGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGt                  <  1:2516014/139‑1 (MQ=255)
                                                                                                                                                ggCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATc   >  1:1671390/1‑150 (MQ=255)
                                                                                                                                                 gCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             >  1:568231/1‑75 (MQ=255)
                                                                                                                                                 gCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCg  <  2:966154/150‑1 (MQ=255)
                                                                                                                                                   gAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      >  1:1461590/1‑64 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCG  >  NC_002947/3584795‑3585088

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