Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,285,613 A→T 7.6% intergenic (+98/+25) 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,6130AT 7.6% 218.5 / 5.4 74intergenic (+98/+25)PP_4658/ggtmethyl‑accepting chemotaxis protein/gamma‑glutamyltranspeptidase
Reads supporting (aligned to +/- strand):  ref base A (35/32);  new base T (2/4);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

GGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAGCCC  >  NC_002947/5285473‑5285761
                                                                                                                                            |                                                                                                                                                    
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                                                                                                                                            >  4:52608/1‑149 (MQ=255)
  gtgggCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCgg                                                                                                                                            >  7:349195/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                                                                                                                                   <  6:231092/133‑1 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGAc                                                                                                                                   >  5:231092/1‑133 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtactac                                                                                                                          >  4:337869/1‑142 (MQ=255)
                           ccAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtactac                                                                                                                          <  3:337869/142‑1 (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                                                                                                                             <  6:269362/114‑1 (MQ=255)
                                                    ccTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGtac                                                                                                                             >  5:269362/1‑114 (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                                                          >  3:27945/1‑149 (MQ=255)
                                                                                    aaGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGa                                                          >  5:77310/1‑149 (MQ=255)
                                                                                       gCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCAt                                                             <  7:90451/143‑1 (MQ=255)
                                                                                       gCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCAt                                                             >  8:90451/1‑143 (MQ=255)
                                                                                         ccccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGAcgcg                                                      >  4:299493/1‑148 (MQ=255)
                                                                                         ccccAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGAcgcg                                                      >  6:107251/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                                   <  8:333699/140‑1 (MQ=255)
                                                                                                                    ccTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGccc                                   >  7:333699/1‑140 (MQ=255)
                                                                                                                    ccTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAg                          <  4:286014/149‑1 (MQ=255)
                                                                                                                      ttGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCa                                                                          <  8:226511/99‑1 (MQ=255)
                                                                                                                      ttGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCa                                                                          >  7:226511/1‑99 (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)
                                                                                                                                            aGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAGccc  <  3:87743/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
GGGTGGGCAGTTGCAGACACTGGTCAGCCAGTTCCGCGTCTGACCGTTGGCACCTGGGGCCGCTTTGCGGCCCTTTCGCGACTCAAGGCCCCCAACGCAGGTCCTCGTGGGAGCGGCCTTGTGTCGCGATAGGAGGGCAAAGCCCTCCCGGCGTTCTGACAGGTACTACTTGACGATCATCCCTACCCCACGCCCACGTGGGTCGGAAGCAGTCTCAAGCTTCGCCCCATCGACGCGGATCGCCTGGATGTCGCCCATCTCCCAGCCCTGGTCCTCCAGCACATAGCCC  >  NC_002947/5285473‑5285761

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: