Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,584,942 A→C 5.5% 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.5% 313.8 / 4.3 110intergenic (+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 (53/50);  new base C (4/3);  total (57/53)
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.99e-01

Only 100 of 110 total aligned reads displayed.

TGCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCG  >  NC_002947/3584795‑3585088
                                                                                                                                                   |                                                                                                                                                  
tgCGGCCCGCCACCCATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTCATTCTGAGCAATACGCCGACCCTGTGGGAGCGGGCGAGc                                                                                                                                                  <  1:552606/150‑1 (MQ=255)
              cATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                                    >  2:745325/1‑150 (MQ=255)
              cATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                                    >  2:1950156/1‑150 (MQ=255)
              cATGGTGGAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTg                                                                                                                                    >  1:716353/1‑150 (MQ=255)
                 ggtggAAGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATc                                                                                                                                 >  2:1231699/1‑150 (MQ=255)
                      aaGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCa                                                                                                                            >  1:1344700/1‑150 (MQ=255)
                      aaGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCa                                                                                                                            <  1:214018/150‑1 (MQ=255)
                       aGCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGACAGCTCCCACAGTGATCgcgcgg                                                                                                                           >  2:134131/1‑148 (MQ=255)
                        gCGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGc                                                                                                                          <  1:1231699/150‑1 (MQ=255)
                         cGGTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGcc                                                                                                                         <  1:1839769/150‑1 (MQ=255)
                          ggTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCTCAGCCt                                                                                                                        >  2:2283548/1‑150 (MQ=255)
                           gTCAGCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGTATGGTCGGTGAAAGCAATACGCCGACGCTGTGGCACACGCCGAGCCCGCTCCCACAGTGATCGCGCAGCCtt                                                                                                                       >  1:2156148/1‑150 (MQ=255)
                               gCCAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGCCCCTGTGGGAGAGGGCGCGACCGCTCCCACAGTGATCGCGCAGCCTTTCAg                                                                                                                   <  1:1806151/150‑1 (MQ=255)
                                ccAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAgg                                                                                                                  >  1:771083/1‑150 (MQ=255)
                                ccAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAgg                                                                                                                  >  2:2275680/1‑150 (MQ=255)
                                ccAGTACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGGGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAgg                                                                                                                  >  1:1093560/1‑150 (MQ=255)
                                    tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                              <  1:511501/150‑1 (MQ=255)
                                    tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCttt                                                                                                              <  2:2222659/150‑1 (MQ=255)
                                    tACGTGCGTGAGCAGGGCATTACCCCGGCCAACTTCTACGACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCTAAGCAATCCGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTGTCAGGAttt                                                                                                              <  2:804560/150‑1 (MQ=255)
                                       gtgcgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAg                                                                                                           >  1:1011488/1‑150 (MQ=255)
                                          cgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGcc                                                                                                                                                 >  1:1340802/1‑109 (MQ=255)
                                          cgtgAGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGcc                                                                                                                                                 <  2:1340802/109‑1 (MQ=255)
                                              aGCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCGGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAagaga                                                                                                    <  1:134131/150‑1 (MQ=255)
                                               gCAGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAagagag                                                                                                   >  1:56036/1‑150 (MQ=255)
                                                 aGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                                    <  1:309784/131‑1 (MQ=255)
                                                 aGGGCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCa                                                                                                                    >  2:309784/1‑131 (MQ=255)
                                                   ggCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCTTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCcaca                                                                                                                                       <  1:333356/110‑1 (MQ=255)
                                                   ggCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCTTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCcaca                                                                                                                                       >  2:333356/1‑110 (MQ=255)
                                                    gCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgct                                                                                              >  2:656526/1‑150 (MQ=255)
                                                    gCATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgct                                                                                              >  1:1002135/1‑150 (MQ=255)
                                                     cATTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgctg                                                                                             >  1:224074/1‑150 (MQ=255)
                                                      aTTACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgctgc                                                                                            <  2:716353/150‑1 (MQ=255)
                                                      aTTACCCCGGCCAACTTCTACTACGAGAAGTCTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGCCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgctgc                                                                                            <  1:1329544/150‑1 (MQ=255)
                                                       ttACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCca                                                                                                                                         <  2:1018428/104‑1 (MQ=255)
                                                       ttACCCCGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCca                                                                                                                                         >  1:1018428/1‑104 (MQ=255)
                                                          ccccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAg                                                                                                                   >  1:1352660/1‑123 (MQ=255)
                                                            ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  1:745325/150‑1 (MQ=255)
                                                            ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  2:1011488/150‑1 (MQ=255)
                                                            ccGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTAAAAGCAATACGCCGACCCTGTGGGAGCGTGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  2:260727/150‑1 (MQ=255)
                                                             cGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGa                                                                                                            <  1:1223536/127‑1 (MQ=255)
                                                             cGGCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGa                                                                                                            >  2:1223536/1‑127 (MQ=255)
                                                              ggCCAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGt                                                                                    <  2:1389920/150‑1 (MQ=255)
                                                                 cAACTTCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCCAGCCCGCTCCCACAGTGATCGCGCCGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCGCAGGTtctc                                                                                 >  2:340665/1‑150 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGa                                                                                                                                   <  2:1399218/96‑1 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGa                                                                                                                                   >  1:1399218/1‑96 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             >  1:106505/1‑150 (MQ=255)
                                                                     ttCTACTACGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             <  1:2299057/150‑1 (MQ=255)
                                                                       ctactaCGAGAAGTTTGAGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTGCGCACAAGAAATTTGCTGCCACAGGGTCTCTCGCGa                                                                           >  2:1191546/1‑150 (MQ=255)
                                                                          ctacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACt                                                                        <  1:2295032/150‑1 (MQ=255)
                                                                           tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTtgctgc                                                                                            <  2:122606/129‑1 (MQ=255)
                                                                           tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                       <  1:1950156/150‑1 (MQ=255)
                                                                           tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTg                                                                       <  2:649736/150‑1 (MQ=255)
                                                                           tacGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTAAGCAAGAGAGTtgctgc                                                                                            >  1:122606/1‑129 (MQ=255)
                                                                            acGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTgg                                                                      <  1:781503/150‑1 (MQ=255)
                                                                             cGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctc                                                                                 <  1:1295695/138‑1 (MQ=255)
                                                                             cGAGAAGTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGtctc                                                                                 >  2:1295695/1‑138 (MQ=255)
                                                                              gagaAGTTTGCGGCTGCGGCCTGACCCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGgag                                                                    >  1:551396/1‑150 (MQ=255)
                                                                                   gTTTGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTa                                                               >  1:251337/1‑150 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              <  2:224074/150‑1 (MQ=255)
                                                                                    tttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGAGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTAc                                                              >  2:764242/1‑150 (MQ=255)
                                                                                     ttGCGGCTGCGGCCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACa                                                             <  1:656526/150‑1 (MQ=255)
                                                                                     ttGCGGCTGCGGCCTGACGCTCCCTGATGCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACa                                                             <  1:441529/150‑1 (MQ=255)
                                                                                         ggctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATg                                                         <  2:56036/150‑1 (MQ=255)
                                                                                         ggccgcggcCTGACGAACTCCGATTCAAAGCAATAAACCGCCACTGCCGGGGGCGGCGCGCCCGCGACCACAGAGATCCCGCAGCCTTTGAAGCATTGAAGAAGAGAGTTGCTGCCACAGGGTCCCTCGCATACTGAAAAAATACATATg                                                         >  2:1479852/1‑150 (MQ=255)
                                                                                          gctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGa                                                        <  2:771083/150‑1 (MQ=255)
                                                                                          gctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGa                                                        >  2:1306782/1‑150 (MQ=255)
                                                                                           ctgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCACCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaa                                                       <  1:1191546/150‑1 (MQ=255)
                                                                                            tgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc                                                                                                                                           <  2:2038648/65‑1 (MQ=37)
                                                                                            tgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCTCGCCCGCTcc                                                                                                                                           >  1:2038648/1‑65 (MQ=255)
                                                                                            tgcggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaac                                                      <  2:27502/150‑1 (MQ=255)
                                                                                               ggcCTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGaacaac                                                   <  2:147223/150‑1 (MQ=255)
                                                                                                 ccTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCt                                                                                                                                             >  2:2292209/1‑58 (MQ=255)
                                                                                                 ccTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCt                                                                                                                                             <  1:2292209/58‑1 (MQ=255)
                                                                                                 ccTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAg                                                 <  2:922700/150‑1 (MQ=255)
                                                                                                 ccTGACGCTCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAg                                                 >  1:922700/1‑150 (MQ=255)
                                                                                                         tCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCTCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaag                                         <  1:2283548/150‑1 (MQ=255)
                                                                                                         tCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaag                                         <  2:2132957/150‑1 (MQ=255)
                                                                                                         tCCCTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaag                                         <  1:1740011/150‑1 (MQ=255)
                                                                                                           ccTGATTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggc                                       >  2:1557824/1‑150 (MQ=255)
                                                                                                               aTTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAg                                                                                       >  2:664887/1‑98 (MQ=255)
                                                                                                               aTTCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCcaaggcaagg                                   <  2:1002135/150‑1 (MQ=255)
                                                                                                                 tCAAAGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAgg                                                                                      <  1:955592/97‑1 (MQ=255)
                                                                                                                   aaaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             >  2:146582/1‑104 (MQ=255)
                                                                                                                   aaaGCAATACGCCGACCCTGTGGCAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTcgc                                                                             <  1:146582/104‑1 (MQ=255)
                                                                                                                    aaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGcc                              >  2:1659715/1‑150 (MQ=255)
                                                                                                                    aaGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGcc                              >  1:1828824/1‑150 (MQ=255)
                                                                                                                     aGCAATACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGccc                             >  1:1165045/1‑150 (MQ=255)
                                                                                                                        aaTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTgg                          >  1:982863/1‑150 (MQ=255)
                                                                                                                         aTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcac                                                                                         <  2:1282316/86‑1 (MQ=255)
                                                                                                                         aTACGCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCcac                                                                                         >  1:1282316/1‑86 (MQ=255)
                                                                                                                             gCCGACCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGGCAGAGGGCCTTTGACCAGTACGAGATAAACAGATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCAcc                     >  1:274546/1‑150 (MQ=255)
                                                                                                                                 aCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAg                                                                                       >  1:1130982/1‑80 (MQ=255)
                                                                                                                                 aCCCTGTGGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAg                                                                                       <  2:1130982/80‑1 (MQ=255)
                                                                                                                                     tgtgGGAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTgccgc             >  1:258110/1‑150 (MQ=255)
                                                                                                                                       tgGGAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGccc           >  1:1253768/1‑150 (MQ=255)
                                                                                                                                         gtaGCGTGAGAGCCCGCTCACACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACATATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCaa         <  1:340665/148‑1 (MQ=255)
                                                                                                                                          gAGCGGGCTCGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGt                                  <  1:408196/124‑1 (MQ=255)
                                                                                                                                          gAGCGGGCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAg        >  1:1936088/1‑150 (MQ=255)
                                                                                                                                               ggCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATc   >  1:1969968/1‑150 (MQ=255)
                                                                                                                                                gCGAGCCCGCTCCCACAGTGATCGCGCAGCCTTTCAGGCTTTGAGCAAGAGAGTTGCTGCCACAGGGTCTCTCGCGAACTGGAGATTTACATATGAACAACAGATTCCAAGGCAAGGTTGCCCTGGTCACCGGTGCCGCCCAAGGCATCg  <  1:1306782/150‑1 (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 < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: