Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,584,942 A→C 5.6% 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 5.6% 266.2 / 6.6 89intergenic (+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 (40/44);  new base C (2/3);  total (42/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.05e-01

TGTGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAA  >  NC_002947/3584793‑3585081
                                                                                                                                                      |                                                                                                                                           
tgtgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGa                                                                                                                                             <  2:236112/150‑1 (MQ=255)
tgtgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGa                                                                                                                                             <  1:1405665/150‑1 (MQ=255)
  tgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGTATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCGGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGc                                                                                                                                           <  1:2295431/150‑1 (MQ=255)
       cccGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCt                                                                                                                                      <  2:2361719/150‑1 (MQ=255)
        ccGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTc                                                                                                                                     <  1:986525/150‑1 (MQ=255)
         cGCCACCAATCGTAGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCGGGCACAGGCGCGAGGCGATAGGTTTGCGGCTCCGGCCTGACGCTCCCTCCTTCAAAGTATGCCGCCGACCCTGTGGGAGCGGGCGAGCCCGCTcc                                                                                                                                    <  1:1749465/150‑1 (MQ=255)
             ccccATGGTGGAAGCGGTCAGCCAGTACGGGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCcaca                                                                                                                                <  2:2360550/149‑1 (MQ=255)
               ccATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGt                                                                                                                              >  1:681045/1‑150 (MQ=255)
                cATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                             >  2:2102509/1‑150 (MQ=255)
                     tggAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATcgc                                                                                                                        <  2:1462845/150‑1 (MQ=255)
                        aaGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCa                                                                                                                     <  1:2209909/150‑1 (MQ=255)
                              tCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGATGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCttt                                                                                                               <  1:138088/150‑1 (MQ=255)
                                aGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                             >  1:660785/1‑150 (MQ=255)
                                aGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                             >  2:2380337/1‑150 (MQ=255)
                                 gCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAg                                                                                                            <  2:691026/150‑1 (MQ=255)
                                  ccAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCg                                                                                                                                        >  1:1025303/1‑121 (MQ=255)
                                  ccAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCg                                                                                                                                        <  2:1025303/121‑1 (MQ=255)
                                      tACGTGCGTGAGCAGGGCATTCCCCCGTCCAACTTCTAGTAGGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                       <  1:2061260/150‑1 (MQ=255)
                                      tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                       <  1:2380337/150‑1 (MQ=255)
                                      tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                       >  2:1965766/1‑150 (MQ=255)
                                      tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                       <  1:1545875/150‑1 (MQ=255)
                                        cgtgcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGa                                                                                                     >  2:2001596/1‑150 (MQ=255)
                                          tgcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGc                                                                                                   <  1:401301/150‑1 (MQ=255)
                                              tgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCcaca                                                                                                                                >  2:764785/1‑117 (MQ=255)
                                              tgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCcaca                                                                                                                                <  1:764785/117‑1 (MQ=255)
                                                  cAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGt                                                                                           <  1:2102509/150‑1 (MQ=255)
                                                   agtgCATCACCCCGGCCAACATCTACGACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGATCGGGCGAGACCGCGCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAATAGAGtt                                                                                          <  2:818390/147‑1 (MQ=255)
                                                     ggCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgc                                                                                        >  1:2369550/1‑150 (MQ=255)
                                                     ggCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGACGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgc                                                                                        >  1:1194854/1‑150 (MQ=255)
                                                       cATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGACTtgctg                                                                                      <  2:1901030/150‑1 (MQ=255)
                                                       cATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCAGACCCTGTGGGAGCGGGCGAGCCCGc                                                                                                                                       >  2:626288/1‑101 (MQ=255)
                                                       cATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCAGACCCTGTGGGAGCGGGCGAGCCCGc                                                                                                                                       <  1:626288/101‑1 (MQ=255)
                                                         ttACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGcc                                                                                    >  2:1995982/1‑150 (MQ=255)
                                                           aCCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGTGCAAGAGAGTTGCTGCcac                                                                                  >  2:1161908/1‑150 (MQ=255)
                                                            ccccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                       <  2:1308151/128‑1 (MQ=255)
                                                            ccccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                       >  1:1308151/1‑128 (MQ=255)
                                                            ccccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcaca                                                                                 >  1:335855/1‑150 (MQ=255)
                                                             cccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTGAGCCTGTGAGCAAGAGAGTTGCTGCCACAg                                                                                >  1:2365/1‑150 (MQ=255)
                                                              ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                               <  1:2001596/150‑1 (MQ=255)
                                                              ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                               <  1:2369082/150‑1 (MQ=255)
                                                              ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGCCGCGGCCACAgg                                                                               <  1:329393/150‑1 (MQ=255)
                                                                    cAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctc                                                                          >  2:1266034/1‑150 (MQ=255)
                                                                    aaaCTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGGCACAGGGtctc                                                                          >  2:580598/2‑150 (MQ=255)
                                                                      aCTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctctc                                                                        <  1:1266034/150‑1 (MQ=255)
                                                                           tactacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGaa                                                                   <  1:2391719/150‑1 (MQ=255)
                                                                              tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCaa                                                                                                 >  2:1153183/1‑117 (MQ=255)
                                                                              tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCaa                                                                                                 <  1:1153183/117‑1 (MQ=255)
                                                                              tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                <  1:621307/150‑1 (MQ=255)
                                                                              tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                >  2:621307/1‑150 (MQ=255)
                                                                                cGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGga                                                              >  2:1352101/1‑150 (MQ=255)
                                                                                cGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGga                                                              <  2:974559/150‑1 (MQ=255)
                                                                                cGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGga                                                              <  1:1995982/150‑1 (MQ=255)
                                                                                 gagaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtct                                                                           >  1:1066073/1‑136 (MQ=255)
                                                                                 gagaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtct                                                                           <  2:1066073/136‑1 (MQ=255)
                                                                                 gagaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGgag                                                             >  2:326116/1‑150 (MQ=255)
                                                                                  agaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCt                                                                                                                                      <  2:2397910/76‑1 (MQ=38)
                                                                                  agaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCt                                                                                                                                      >  1:2397910/1‑76 (MQ=38)
                                                                                   gaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                               <  2:408066/130‑1 (MQ=255)
                                                                                   gaAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                               >  1:408066/1‑130 (MQ=255)
                                                                                        ttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTCCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACa                                                      <  2:681045/150‑1 (MQ=255)
                                                                                         tgcggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGACGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACat                                                     <  2:1194854/150‑1 (MQ=255)
                                                                                           cggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACatat                                                   <  2:2369550/150‑1 (MQ=255)
                                                                                               tgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaac                                               >  1:518987/1‑150 (MQ=255)
                                                                                                 cggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCCCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                      >  2:557830/1‑125 (MQ=255)
                                                                                                 cggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCCCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                      <  1:557830/125‑1 (MQ=255)
                                                                                                  ggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaacaac                                            <  1:326116/150‑1 (MQ=255)
                                                                                                   gcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAacaaca                                           >  2:1248692/1‑150 (MQ=255)
                                                                                                          gCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCttt                                                                                                               <  1:1263808/75‑1 (MQ=255)
                                                                                                          gCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCttt                                                                                                               >  2:1263808/1‑75 (MQ=255)
                                                                                                           cTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaa                                   >  2:1859993/1‑150 (MQ=255)
                                                                                                            tCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaag                                  <  1:1352101/150‑1 (MQ=255)
                                                                                                              ccTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggg                                >  1:2317784/1‑149 (MQ=255)
                                                                                                                tGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCACGGGCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggcaa                              >  1:561116/1‑150 (MQ=255)
                                                                                                                   ttCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCa                                                                                                                     <  1:1743509/60‑1 (MQ=255)
                                                                                                                   ttCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCa                                                                                                                     >  2:1743509/1‑60 (MQ=255)
                                                                                                                      aaaGCAGTAGGCCGACCCTGCTGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGc                        <  2:518987/150‑1 (MQ=255)
                                                                                                                       aaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                             <  1:1204943/64‑1 (MQ=255)
                                                                                                                       aaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                             >  2:1204943/1‑64 (MQ=255)
                                                                                                                         gCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATGTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCt                     >  1:1230628/1‑150 (MQ=255)
                                                                                                                           aaTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTgg                   >  1:277155/1‑150 (MQ=255)
                                                                                                                           aaTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTgg                   >  1:1788093/1‑150 (MQ=255)
                                                                                                                                gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCAcc              >  1:576014/1‑150 (MQ=255)
                                                                                                                                  cGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACAAATGAACTACAGATTCCAAGGCAAGGTTGCCCTGGTCACCgg            >  1:739046/1‑150 (MQ=255)
                                                                                                                                   gACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTg                    >  1:773120/1‑141 (MQ=255)
                                                                                                                                   gACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTg                    <  2:773120/141‑1 (MQ=255)
                                                                                                                                    aCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAAc                                                                  >  1:98354/1‑94 (MQ=255)
                                                                                                                                    aCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAAc                                                                  <  2:98354/94‑1 (MQ=255)
                                                                                                                                        tgtgGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgc      <  2:335855/150‑1 (MQ=255)
                                                                                                                                          tgGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaag                                  >  2:132172/1‑120 (MQ=255)
                                                                                                                                          tgGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaag                                  <  1:132172/120‑1 (MQ=255)
                                                                                                                                            gaaGCGGGCGAGCACGCACCCACAGTGATCGCACGGCGTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCaa  <  2:1230628/148‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                           
TGTGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAA  >  NC_002947/3584793‑3585081

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