Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,785,057 G→A 7.2% intergenic (+25/+81) PP_1587 → / ← dapC NhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,785,0570GA 7.2% 418.1 / 15.8 138intergenic (+25/+81)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (57/71);  new base A (8/2);  total (65/73)
Fisher's exact test for biased strand distribution p-value = 4.59e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

Only 100 of 150 total aligned reads displayed.

CCCTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAG  >  NC_002947/1784908‑1785206
                                                                                                                                                     |                                                                                                                                                      
cccTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCG                                                                                                                                                        >  1:4008825‑M1/1‑150 (MQ=255)
cccTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCG                                                                                                                                                        >  2:3898634‑M1/1‑150 (MQ=255)
 ccTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGgcgc                                                                                                                                                       <  1:5203824/150‑1 (MQ=255)
 ccTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGgcgc                                                                                                                                                       >  1:6138530/1‑150 (MQ=255)
 ccTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGgcgc                                                                                                                                                       <  1:5128459/150‑1 (MQ=255)
   tcaggGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgc                                                                                                                                                     <  2:6063361‑M1/148‑4 (MQ=255)
    gAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGcc                                                                                                                                                    <  2:4393774/150‑1 (MQ=255)
    gAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTCGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGcc                                                                                                                                                    <  1:2136185/150‑1 (MQ=255)
       ggCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGc                                                                                                                                                 <  1:822311/150‑1 (MQ=255)
        gcgcCGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCt                                                                                                                                                <  1:6509654/150‑1 (MQ=255)
         cgcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTc                                                                                                                                               <  1:685293/150‑1 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCCTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  2:1977317/150‑1 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGGGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGcgaa                                                                                                                                              <  2:5671384/150‑4 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTCGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              >  2:4965618/1‑150 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  2:1231100/150‑1 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              >  2:4965458/1‑150 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  2:1231136/150‑1 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  2:3784973/150‑1 (MQ=255)
          gcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  2:3393974/150‑1 (MQ=255)
            aGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGcgcccgcgaagg                                                                                                                                            <  1:363956‑M1/150‑13 (MQ=255)
               ggTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAg                                                                                                                                         <  2:6001732/150‑1 (MQ=255)
               ggTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAg                                                                                                                                         <  2:2927996/150‑1 (MQ=255)
               ggTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAg                                                                                                                                         <  2:3271965/150‑1 (MQ=255)
                gTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCgcgcccgcgcccccccg                                                                                                                                        >  1:6349386/1‑145 (MQ=255)
                  tGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACGGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTa                                                                                                                                      <  1:3840248/150‑1 (MQ=255)
                     aGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActg                                                                                                                                   <  2:5183585/150‑1 (MQ=255)
                      gCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgc                                                                                                                                  <  2:6405653/150‑1 (MQ=255)
                        tGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgc                                                                                                                                >  2:5339566/1‑150 (MQ=255)
                        tGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgc                                                                                                                                >  2:5342889/1‑150 (MQ=255)
                         gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgct                                                                                                                               <  1:3257613/150‑1 (MQ=255)
                         gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgct                                                                                                                               <  1:3186840/150‑1 (MQ=255)
                         gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCACGTAcagcgcg                                                                                                                               >  1:4679873/1‑149 (MQ=255)
                          tACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTAcagcgctg                                                                                                                              >  2:4787466/1‑150 (MQ=255)
                            cAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgctgcg                                                                                                                            >  2:2718813/1‑150 (MQ=255)
                            cAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgctgcg                                                                                                                            >  2:6722435/1‑150 (MQ=255)
                              aCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTgcgctgcgct                                                                                                                          <  1:2190046/150‑1 (MQ=255)
                               ccTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCtt                                                                                                                         <  1:1154451/150‑1 (MQ=255)
                                 tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                       <  2:4909292/150‑1 (MQ=255)
                                  gCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaa                                                                                                                      <  1:4219266/150‑1 (MQ=255)
                                  gCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaa                                                                                                                      <  1:6871150/150‑1 (MQ=255)
                                    aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                    <  2:1330112/150‑1 (MQ=255)
                                       gTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAg                                                                                                                 <  1:6655914/150‑1 (MQ=255)
                                        tCAGCACCTGGTGGGGGTTGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgc                                                                                                                <  1:187069/150‑1 (MQ=255)
                                           gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                             >  2:3307142/1‑150 (MQ=255)
                                           gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                             >  2:6766636/1‑150 (MQ=255)
                                           gCACCTGGGGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTGAAATCAGCGCt                                                                                                             >  2:5816923/1‑150 (MQ=255)
                                             aCCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGt                                                                                                           <  1:5558728/150‑1 (MQ=255)
                                             aCCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGACCACAGGTACTTCGCTGCGCTTGAAAAAAGCGCTGt                                                                                                           >  1:6001335/1‑150 (MQ=255)
                                             aCCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCACACGGAATGCGCGGCGCTTGAAATAAGCGCTGt                                                                                                           >  1:6001469/1‑150 (MQ=255)
                                               cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                         >  2:6221690/1‑150 (MQ=255)
                                                 ggtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                                       <  1:1730990/150‑1 (MQ=255)
                                                  gtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACGGCGCTGCGCTTGAAATCAGCGCGGTCCCtg                                                                                                      >  1:522333/1‑150 (MQ=255)
                                                   tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:4695756/150‑1 (MQ=255)
                                                   tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:3985933/150‑1 (MQ=255)
                                                    gggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgtg                                                                                                    <  2:3670502/150‑1 (MQ=255)
                                                     ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                   <  1:3554545/150‑1 (MQ=255)
                                                         taaccgtggtgCGGCAAGTACTGGGGGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGGGTGAAATCAGCGCTGTACCAGGGGgcgc                                                                                               >  2:4177999/6‑150 (MQ=255)
                                                                   gCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcat                                                                                     >  2:739639‑M1/1‑83 (MQ=255)
                                                                    cGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCCCTTGAAAACAGCGCAGTACCGGTGGGAGCGGGCGTGCCcg                                                                                    >  2:6064351/1‑150 (MQ=38)
                                                                     ggCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACATGTACTGCGCTGCGCTTGAAATAAGCGCTGTACCTGTGGGAGCGGGCGTGCCcgc                                                                                   >  2:3057973/1‑150 (MQ=255)
                                                                      gCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCcgcg                                                                                  >  2:5487065/1‑150 (MQ=255)
                                                                       cAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcgg                                                                                 <  2:806126‑M1/150‑72 (MQ=255)
                                                                       cAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcgg                                                                                 <  2:165926‑M1/150‑72 (MQ=255)
                                                                           tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATc                                                                                                                   >  2:2124105/1‑112 (MQ=255)
                                                                             cTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGTAAGGCGCGGCCCGTGTAAAAAGCGCCGTACCGGTGGGAGCGGGCGTGCCCGCGAAGGGGt                                                                           >  1:5210700/1‑150 (MQ=255)
                                                                              tGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTg                                                                         <  2:6634288/150‑1 (MQ=255)
                                                                                   gAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaa                                                                     <  1:1590155‑M1/150‑84 (MQ=255)
                                                                                       tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCGGCGCGTGGAATCAGCGCGGTTCCGGTGGGGGCGGGCGTGGCCGCGCAGGGGGTGGTGAGCgg                                                                >  1:4283870/1‑150 (MQ=37)
                                                                                           tatGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGCGCGGGCGTGCCCGCGACGGGGTTGATCAGCGGtttt                                                            >  2:6523695/1‑150 (MQ=255)
                                                                                            atGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                             >  1:2899541/1‑101 (MQ=255)
                                                                                            atGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACTGCGCTGCGCGTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGGAGGGGGTGATCAGCCGGtttt                                                            >  2:287120/1‑150 (MQ=38)
                                                                                              gtgtgcgGCGAACCGGGTTTCAGTCACGTCGCCCCTATCACGGCCCTTGGCGGGCGCGCCCGCCCCCCCCGGTACTGCGGGCCCTTTTAAACCCGCGCTGTCCCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCa                                                         <  2:6291731/148‑1 (MQ=255)
                                                                                                 tgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTAATGCGCTGCGCTTTAAATCACCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGCGGTTTTTCAGGa                                                      >  2:1936959/1‑150 (MQ=38)
                                                                                                          ccTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCGCTGCGCTTGGAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGaa                                             >  2:6555375/1‑150 (MQ=255)
                                                                                                           cTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAAt                                            >  2:2774106/1‑150 (MQ=255)
                                                                                                            tGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCCGGTACTGCGGTGCGCGTGGAAGCAGCGCTGTACCTGGGGGGGCGGGCGTGCCCGCGAAGGGGTTGGTCAGCGGTTTTGCAGGAAGGCGCGCATg                                           >  2:762457/1‑150 (MQ=37)
                                                                                                            tGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGGCCCCAGGTACTGCCCTGCGCGTGAAATCAGCGCGGTTACGGTGGGGGCGGGCGGGCCCCCGCAGGGGGTGTTGAGCGGTTTTGGAGGCAGGCGCGAATg                                           >  1:2740959/1‑150 (MQ=255)
                                                                                                             gggTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGGTCAGCGGTTTTGCAGGAAGGCGCGAATgc                                          >  2:2247382/1‑150 (MQ=255)
                                                                                                                   aGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccggg                                     >  2:3070938‑M1/1‑35 (MQ=255)
                                                                                                                   aGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCGGTACCAGTGGGCGCGGGCGTGCCCGCGAAGGGGTTTATCAGCGGTTTTTCAGGTAGGCGCGGATGCGCTccg                                    >  1:1073017/1‑150 (MQ=37)
                                                                                                                   aGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACTGCGCTGCGCGTGGAAACCGCGCGGTAACTGGGGGCGCGGGCGTGCCCGCGCAGGGGTGGATCAGCGGGTTTGGAGGCAGGCGCGGATGCGCTccg                                    >  2:3050319/1‑150 (MQ=25)
                                                                                                                     tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTACTGCGCTGCGCGTGGAAACAGCGCTGTACCAGTGGGAGCGGGCGGGCCCGCGGAGGGGGTGTTGAGCGGTTTTTCAGGAAGGCGCGAATTCGCgccgcc                                  >  1:5536307/1‑150 (MQ=255)
                                                                                                                              cccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGGAATCCGCGCTGTTCCTGTGGGCGCGGGCGTGCCCGCGCAGGGGGTGATCAGCGGTTTTGCAGGCAGGCGCGAATTCGCTCCGCCGCTTCGagg                         >  2:4836016/1‑148 (MQ=38)
                                                                                                                                  aTCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACTGCCCTGCGCTTGAAATCAGCGCTGGACCGGGGGGCGCGGGCGGGCCCGCGCACGGGTTGATCAGCGGTTTTGCACGAAGGCGCGCATGCGCTCCGCCGCGACCATGCagt                     >  2:6424805/1‑148 (MQ=255)
                                                                                                                                   tCATGGCCCTGCGGGGGGGCGCCCGCTCCCACTGTTACTGCGCTGCGATTTAAATCCGCGCTTTACGTGGGGGCGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTc                    <  2:3879124/150‑1 (MQ=255)
                                                                                                                                        ggggggcgcgggcgtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgc                <  1:4965618‑M1/144‑137 (MQ=255)
                                                                                                                                        cggggccgcgggcgggcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgc                <  1:4965458‑M1/144‑137 (MQ=255)
                                                                                                                                           cTTCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgac             >  2:6866153‑M1/1‑11 (MQ=255)
                                                                                                                                            ttCGCGGGCGCGCCCGCGCCCACAGGTACAGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCgg           >  2:4358519/1‑150 (MQ=255)
                                                                                                                                            ttCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgaca            >  1:4106580‑M1/1‑10 (MQ=255)
                                                                                                                                            ttCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgaca            >  1:6146300‑M1/1‑10 (MQ=255)
                                                                                                                                            ttCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgaca            >  1:6486461‑M1/1‑10 (MQ=255)
                                                                                                                                             tCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:374478/60‑1 (MQ=255)
                                                                                                                                             tCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacag           >  2:5885820‑M1/1‑9 (MQ=255)
                                                                                                                                             tCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacag           >  2:2673862‑M1/1‑9 (MQ=255)
                                                                                                                                             tCGCGGGCAcgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacag           >  2:5885420‑M1/1‑9 (MQ=255)
                                                                                                                                               gggggggtccccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggt         <  1:1936959‑M1/148‑144 (MQ=255)
                                                                                                                                                 gggCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggtac       >  1:3554135‑M1/1‑5 (MQ=255)
                                                                                                                                                   gCGtgcccgcgaaggggttgatcagcggttttgcaggaaggcgcgaatgcgctccgccgcttcgatgcattcggccagcggtgcaaccagtgccatgcgcacgcgcccggcacccgggttcaccccatccacttcacgcgacaggtacga     >  2:6278965‑M1/1‑3 (MQ=255)
                                                                                                                                                     gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg  >  1:6286521/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
CCCTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAG  >  NC_002947/1784908‑1785206

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