Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,285,612 A→T 8.1% intergenic (+97/+26) PP_4658 → / ← ggt methyl‑accepting chemotaxis protein/gamma‑glutamyltranspeptidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,285,6120AT 8.1% 210.9 / 6.6 74intergenic (+97/+26)PP_4658/ggtmethyl‑accepting chemotaxis protein/gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (36/32);  new base T (2/4);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

GGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAGC  >  NC_002947/5285473‑5285759
                                                                                                                                           |                                                                                                                                                   
gggtgggCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTccc                                                                                                                                            >  3:354467/1‑149 (MQ=255)
 ggtgggCAGTTGCAGACACTGGTCAGTCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCTTTGCCCTccta                                                                                                                                           <  6:310407/149‑3 (MQ=255)
 ggtgggCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCTTTGCCCTccta                                                                                                                                           <  2:326908/149‑3 (MQ=255)
 ggtgggCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCg                                                                                                                                           >  8:3310/1‑149 (MQ=255)
  gtgggCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCgg                                                                                                                                          >  7:349195/1‑149 (MQ=255)
  gtgggCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCgg                                                                                                                                          >  4:52608/1‑149 (MQ=255)
       cAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCGTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTc                                                                                                                                     >  5:94961/1‑149 (MQ=255)
            gCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGa                                                                                                                                  >  8:269249/1‑147 (MQ=255)
                  acTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCTTTGCCCTcc                                                                                                                                             <  4:367858/130‑1 (MQ=255)
                       tCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACtt                                                                                                                      >  6:49448/1‑148 (MQ=255)
                         aGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGa                                                                                                                    <  8:349195/148‑1 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGAc                                                                                                                                 >  5:231092/1‑133 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGAc                                                                                                                                 <  6:231092/133‑1 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtactac                                                                                                                        <  3:337869/142‑1 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtactac                                                                                                                        >  4:337869/1‑142 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGa                                                                                                                 >  3:286014/1‑149 (MQ=255)
                                  cgcgggtGAGCGTTGGCACGTGGGGCCGGTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAAGGCAGGTCCTCGTGGGAGCGTGCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATcc                                                                                                           <  4:44048/142‑1 (MQ=255)
                                      tCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTAcc                                                                                                      <  6:342190/149‑1 (MQ=255)
                                         gACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACccca                                                                                                   >  3:330781/1‑149 (MQ=255)
                                          aCCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACcccac                                                                                                  <  7:268962/149‑1 (MQ=255)
                                           ccGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGTCAAAGCCCTCCCGTCGTTCTGACAGGTACTACTTGACGATCATCCCTACcccacg                                                                                                 <  8:77539/149‑1 (MQ=255)
                                            cGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACcccacgc                                                                                                >  7:32966/1‑149 (MQ=255)
                                              ttGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACcccacgccc                                                                                              <  5:49448/149‑1 (MQ=255)
                                                 gCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAACGCCCTCCCGGCGTTCTGACAGGTACTACTTGCCGATCATCCCTACcccacgcccacg                                                                                           >  4:218778/1‑149 (MQ=255)
                                                    ccTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtac                                                                                                                           >  5:269362/1‑114 (MQ=255)
                                                    ccTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtac                                                                                                                           <  6:269362/114‑1 (MQ=255)
                                                           ccGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCg                                                                                                                                        >  4:234196/1‑94 (MQ=255)
                                                           ccGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCg                                                                                                                                        <  3:234196/94‑1 (MQ=255)
                                                           ccGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGaa                                                                                 <  4:330781/149‑1 (MQ=255)
                                                           ccGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGaa                                                                                 <  8:32966/149‑1 (MQ=255)
                                                                 tGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGtc                                                                           >  8:66287/1‑149 (MQ=255)
                                                                 tGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGt                                                                            <  3:52608/148‑1 (MQ=255)
                                                                      cccTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCTCGATAGGACGGCGAAGCCCTCCCGGCGTTCTGGCAGGTACTACTTGACGATCATCCATACCCCACGCCCACGTAGGTAGGAAGCAGTCTAAAg                                                                      >  2:105541/1‑149 (MQ=255)
                                                                             gcgACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGcc                                                               <  7:3310/149‑1 (MQ=255)
                                                                              cgACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGccc                                                              >  3:364935/1‑149 (MQ=255)
                                                                                    aaGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGa                                                        >  5:77310/1‑149 (MQ=255)
                                                                                    aaGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGa                                                        >  3:27945/1‑149 (MQ=255)
                                                                                       gCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCAt                                                           >  8:90451/1‑143 (MQ=255)
                                                                                       gCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCAt                                                           <  7:90451/143‑1 (MQ=255)
                                                                                         ccccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGAcgcg                                                    >  6:107251/1‑148 (MQ=255)
                                                                                         ccccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGAcgcg                                                    >  4:299493/1‑148 (MQ=255)
                                                                                         ccccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCgg                                                   <  4:354467/149‑1 (MQ=255)
                                                                                          cccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCCGTCTCAAGCTTCGCCCCATCGACGCGGa                                                  >  3:180214/1‑149 (MQ=255)
                                                                                           ccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGAt                                                 >  3:360692/1‑149 (MQ=255)
                                                                                                gCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCt                                            <  6:77310/149‑1 (MQ=255)
                                                                                                    gTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGtc                                                                           <  2:158328/114‑1 (MQ=255)
                                                                                                        tCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCg                                    >  2:192819/1‑149 (MQ=255)
                                                                                                          gTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGc                                   <  7:269249/148‑1 (MQ=255)
                                                                                                           tGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGccc                                 <  1:227619/149‑1 (MQ=255)
                                                                                                            gggAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCa                                >  6:93166/1‑149 (MQ=255)
                                                                                                                gCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCAtctc                            >  5:8/1‑149 (MQ=255)
                                                                                                                   gCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTc                                     <  3:203053/137‑1 (MQ=255)
                                                                                                                   gCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTc                                     >  4:203053/1‑137 (MQ=255)
                                                                                                                   gCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTcccg                         >  6:314984/1‑148 (MQ=255)
                                                                                                                   gCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCa                         <  6:8/149‑1 (MQ=255)
                                                                                                                    ccTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGccc                                 >  7:333699/1‑140 (MQ=255)
                                                                                                                    ccTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGccc                                 <  8:333699/140‑1 (MQ=255)
                                                                                                                    ccTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAg                        <  4:286014/149‑1 (MQ=255)
                                                                                                                      ttGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCa                                                                        >  7:226511/1‑99 (MQ=255)
                                                                                                                      ttGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCa                                                                        <  8:226511/99‑1 (MQ=255)
                                                                                                                        gtgtCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCGGccc                     <  5:314984/148‑1 (MQ=255)
                                                                                                                         tgtCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGACTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTg                   >  1:5779/1‑149 (MQ=255)
                                                                                                                           tCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTgg                  <  4:364935/148‑1 (MQ=255)
                                                                                                                           tCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGt                 >  4:87743/1‑149 (MQ=255)
                                                                                                                             gcgATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGtc                <  7:66287/148‑1 (MQ=255)
                                                                                                                             gcgATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGtc                <  5:107251/148‑1 (MQ=255)
                                                                                                                                 tAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCa           <  3:240209/149‑1 (MQ=255)
                                                                                                                                  gggaggGCTTTGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAg          >  6:302541/2‑149 (MQ=255)
                                                                                                                                   ggaggGCTTTGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGAt                                        <  1:335978/118‑1 (MQ=255)
                                                                                                                                   ggaggGCTTTGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGAt                                        >  2:335978/1‑118 (MQ=255)
                                                                                                                                    gaggGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCg                                    >  4:368149/1‑121 (MQ=255)
                                                                                                                                    gaggGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCg                                    <  3:368149/121‑1 (MQ=255)
                                                                                                                                    gaggGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGca        >  4:292950/1‑149 (MQ=255)
                                                                                                                                         cAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAg   <  4:27945/149‑1 (MQ=255)
                                                                                                                                          aaaGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAGc  <  1:192819/149‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                                   
GGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAGC  >  NC_002947/5285473‑5285759

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