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% 388.2 / 14.7 126intergenic (+25/+81)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (44/72);  new base A (2/7);  total (46/80)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.36e-01

Only 100 of 130 total aligned reads displayed.

CCTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAG  >  NC_002947/1784909‑1785206
                                                                                                                                                    |                                                                                                                                                      
ccTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGgcgc                                                                                                                                                       <  2:6936669/150‑1 (MQ=255)
  tGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGgcgcgc                                                                                                                                                     >  2:443870/1‑150 (MQ=255)
   gAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGcc                                                                                                                                                    <  1:3007670/150‑1 (MQ=255)
   gAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGcc                                                                                                                                                    <  1:5638200/150‑1 (MQ=255)
     gggCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCg                                                                                                                                                  <  1:3521138/150‑1 (MQ=255)
      ggCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGTGCCCGc                                                                                                                                                 <  1:1223293/150‑1 (MQ=255)
        cgcAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTc                                                                                                                                               <  2:6161463/150‑1 (MQ=255)
        cgcAGCGGTTGGAGCTGAACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGTGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTc                                                                                                                                               <  2:4560671/150‑1 (MQ=255)
         gcAGCGGTTGGAGCTGTTCAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  1:2256996/150‑1 (MQ=255)
         gcAGCGGTTGGAGCTGTACAACCTGCATCGGCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCACGGCCCTTCGCGGGCGCGCCCGCTcc                                                                                                                                              <  1:6896693/150‑1 (MQ=255)
           ggCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGCACCTGGGTTCAGTCCAGTAGACCCTATCAAGGCCCGTCGCGGGCGCGCCCGCTCCca                                                                                                                                            <  2:367765/149‑1 (MQ=255)
           aGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTcccc                                                                                                                                            >  1:3720399/1‑149 (MQ=255)
           aGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCca                                                                                                                                            <  2:1804171/150‑1 (MQ=255)
           aGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCca                                                                                                                                            <  2:1804158/150‑1 (MQ=255)
            gCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCCAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCcac                                                                                                                                           >  2:4427714/1‑150 (MQ=255)
            gCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCcac                                                                                                                                           <  2:5762205/150‑1 (MQ=255)
            gCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCcac                                                                                                                                           <  1:7406455/150‑1 (MQ=255)
             cGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCcaca                                                                                                                                          <  1:5277675/150‑1 (MQ=255)
               gTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAgg                                                                                                                                        <  1:443870/150‑1 (MQ=255)
                 tGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTa                                                                                                                                      <  1:1580216/150‑1 (MQ=255)
                 tGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTa                                                                                                                                      <  1:1185593/150‑1 (MQ=255)
                  ggAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAc                                                                                                                                     <  2:3804250/150‑1 (MQ=255)
                   gAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAca                                                                                                                                    <  2:6273917/150‑2 (MQ=255)
                    aGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActg                                                                                                                                   <  1:449673/150‑1 (MQ=255)
                     gCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACCGGTActgc                                                                                                                                  >  1:4808374/1‑150 (MQ=255)
                     gCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgc                                                                                                                                  >  2:2279213/1‑150 (MQ=255)
                     gCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTActtc                                                                                                                                  >  1:6290496/1‑148 (MQ=255)
                      cTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcg                                                                                                                                 <  1:802943/150‑1 (MQ=255)
                        gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGTGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgct                                                                                                                               <  1:3216908/150‑1 (MQ=255)
                        gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgct                                                                                                                               <  1:2751868/150‑1 (MQ=255)
                        gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgct                                                                                                                               <  1:3157098/150‑1 (MQ=255)
                         tACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgctg                                                                                                                              >  2:50549/1‑150 (MQ=255)
                         tACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTActgcgctg                                                                                                                              >  1:4444385/1‑150 (MQ=255)
                          aCAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgctgc                                                                                                                             <  2:3081627/150‑1 (MQ=255)
                            aaCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACtgcgctgcgc                                                                                                                           <  2:1288975/150‑1 (MQ=255)
                            aaCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACtgcgctgcgc                                                                                                                           >  1:6724728/1‑150 (MQ=255)
                              ccTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCtt                                                                                                                         <  2:3831130/150‑1 (MQ=255)
                                tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGGTatgcgctgcgggtta                                                                                                                       >  1:6538663/1‑145 (MQ=255)
                                     cGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGTTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCa                                                                                                                  <  1:14323/150‑1 (MQ=255)
                                      gTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAg                                                                                                                 <  2:49256/150‑1 (MQ=255)
                                      gTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAg                                                                                                                 <  2:5811873/150‑1 (MQ=255)
                                      gTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAg                                                                                                                 <  2:1495737/150‑1 (MQ=255)
                                      gTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACTGCGCTGCGCATGAAATCAg                                                                                                                 >  2:3252185/1‑150 (MQ=255)
                                       tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGCTTGAAATCAgc                                                                                                                >  2:464671/1‑150 (MQ=255)
                                       tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGTACAGCGCAGCGCGTGaaaaaagc                                                                                                                >  1:6991963/1‑145 (MQ=255)
                                        cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCACGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                               <  1:3745807/150‑1 (MQ=255)
                                        cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTAAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                               <  1:1407433/150‑1 (MQ=255)
                                         aGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcgc                                                                                                              <  1:245996/150‑1 (MQ=255)
                                          gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                             <  2:4919880/150‑1 (MQ=255)
                                          gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCt                                                                                                             >  2:3600981/1‑150 (MQ=255)
                                           cACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTg                                                                                                            <  2:4967738/150‑1 (MQ=255)
                                             ccTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                          <  1:5325018/150‑1 (MQ=255)
                                             ccTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCCCTTCGCTTGAAATCAGCGCTGTa                                                                                                          >  1:2378193/1‑150 (MQ=255)
                                              cTGGTGGGGGATGAGGTGGTGCGGCAAGTAGTGGGTGATATTGATATGAGTGAGGCGACCCTGGGTTCAGTCAAGTAGTCGCTATGTAGTTCCTTCTTGGGCGCGCCCGCTCCCACAGGTTCTGCTCTGCGGTTGAAAGCAGCGCTGTAc                                                                                                         <  2:7202934/150‑1 (MQ=255)
                                              cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGTATATGAGTGAGGCGAACCTTGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTACCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                         <  1:2417780/150‑1 (MQ=255)
                                              cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCATGGCCCTTCGCGGGCGCGCCCGCTACCACAGGTACTGCGGTGCGCTTTAAATAAGCGCTGTAc                                                                                                         >  2:2072958/1‑150 (MQ=255)
                                              cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAc                                                                                                         <  1:2144801/150‑1 (MQ=255)
                                               tggtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACCGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAcc                                                                                                        <  1:3039060/150‑1 (MQ=255)
                                                ggtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                                       <  2:3919903/150‑1 (MQ=255)
                                                  tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  2:660850/150‑1 (MQ=255)
                                                  tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:5586965/150‑1 (MQ=255)
                                                  tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:6548355/150‑1 (MQ=255)
                                                  tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:5089347/150‑1 (MQ=255)
                                                   gggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgtg                                                                                                    <  2:6042297/150‑1 (MQ=255)
                                                    ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                   <  2:2315947/150‑1 (MQ=255)
                                                    ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                   <  1:346036/150‑1 (MQ=255)
                                                    ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                   <  1:5328374/150‑1 (MQ=21)
                                                    ggggATGAGGTGGAGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                   <  2:6688368/150‑1 (MQ=255)
                                                       gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAg                                                                                                <  1:2006562/150‑1 (MQ=21)
                                                       gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAg                                                                                                <  1:1107221/150‑1 (MQ=21)
                                                        aggAGGTGGTGCGGCAAGTCCGGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCACGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGc                                                                                               <  1:1896507/148‑1 (MQ=255)
                                                          gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCgg                                                                                             <  1:1143418/150‑1 (MQ=255)
                                                          gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCgg                                                                                             <  1:2826111/150‑1 (MQ=255)
                                                           aGGTGGTGCGGCACGTACTGGGTGAGCTGGATATGCGACAGGCGCACCTGGGTTCAGTCCAGACGACCCTACAACGGCCCTTCGCGCCCCCGCCCGCCCCCACAGGGACTGCGCTGCGCTCGAAATCAGCGCTGTACCTGTCGGACCggg                                                                                            <  2:2927344/150‑1 (MQ=255)
                                                            ggtggtGCGGCAAGTAGTGGGTGTGCTGTATATGAGTGCGGCGACCCTGGGTTCAGTAAAGTAGACCCTATCACGGCCCTTGGCGGGCGCGCCCGCTCCCACTGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGc                                                                                           <  1:4840686/150‑1 (MQ=255)
                                                                                      tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGGACGGCCCTGCGCGTGGAATCACCCCAGTAACAGTGGGAGCGGGCGGGCCCGCGGAGGGGGTGaaaaccgg                                                                >  2:1233644/1‑143 (MQ=255)
                                                                                      tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTAAGGCGCGGCGCGTGCAATCAGCGCGGTAACAGTGGGGGCGGGCGGGCCCGCCGACGGGGTTGTCAGCgg                                                                >  2:1233609/1‑150 (MQ=255)
                                                                                          tatGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACTGCGCTGCGCTTGAAATCAGCCCTGTACCGGGGGGAGCGGGCGTGCCCGCGAAGCGGTTGATCAGCGGtttt                                                            >  2:2626333/1‑150 (MQ=39)
                                                                                                 gagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     >  2:4811360/1‑103 (MQ=255)
                                                                                                 gagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                     <  1:4811360/103‑1 (MQ=255)
                                                                                                   ggCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                          >  2:4699347/1‑96 (MQ=255)
                                                                                                       aaCCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGAAGGcgcg                                               >  2:2291003/1‑150 (MQ=255)
                                                                                                         ccTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACCGGTACGGCCCTGCGCGTGAAATCACCGCAGTACCAGTGGGACCGGGCGTGCCCGCGCCGGGCTTGATCCGCCGGTTTGCAGGATGGcgcgcc                                             >  2:3740593/1‑148 (MQ=25)
                                                                                                            gggTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATgc                                          >  2:5561313/1‑150 (MQ=255)
                                                                                                               ttCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCt                                       >  2:6663333/1‑150 (MQ=255)
                                                                                                               ttCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTTCTGGGCTGCGCTGGGAATAAGCGCGGTAACAGTGGGCGCGGGCGTGCCCCCGGAGGGGTTGTTCAGCGGTTTTTGAGGAAGGCGCGCATgcgca                                       >  1:1116152/1‑149 (MQ=255)
                                                                                                                   gTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGGAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTccgc                                   >  2:1777748/1‑150 (MQ=255)
                                                                                                                   gTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACTTCGCTGCGGGTGAAATCAGCGCAGTACCGGGGGGGGCGGGCGTGCCCGCGGCGGGGGTGATCAGCGGGTTTGCCGGAACGCGCGAAGGCGCTccgc                                   >  2:5669955/1‑150 (MQ=255)
                                                                                                                    tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACCGGTACAGCGCGGCCCGTGACATCCGCGCGGTACCGGTGGGAGCCGGCGTGCCCGCGAAGGGGGTGATCACCGGGTTTGCCGGAAGGCGCGCAGGCGCGccccc                                  >  2:1947314/1‑150 (MQ=255)
                                                                                                                       aGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTACTGCGCTGCGCGTGTAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGTTCAGCGGTTTTTCAGGAAGGCGCGACAGCGCTCCGCCGCt                               >  1:5524781/1‑150 (MQ=37)
                                                                                                                              ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTAACTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGGTTTGCAGGAAGGCGCGAATGCGCACCGCCGCGACGATGc                        >  2:6983717/1‑150 (MQ=38)
                                                                                                                              ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACTGCGCTGCGCTTGAAATCAGCCCGGTACCGGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGCGGGTTTGCAGGAACGCGCGAATGCGCTCCGCCGCGTCGATGc                        >  2:5590999/1‑150 (MQ=37)
                                                                                                                               cTATCAAGGCCCTTCGCGGGCGCGCCCGCTACCACAGGTACTTCGCTGCGCTTGAAATCAGCGCTGTTCCAGTGGGAGCGGGCGGGCCCGCGAACGGGGTGTTCAGCGGTTTTGAAGGAAGGCGCGAATGCGCACCCCCCCGTCGATgga                       >  1:6052523/1‑148 (MQ=255)
                                                                                                                                    aaGGCCCTTCGCGGGCGCGCCCGCTCCCACCGGTACAGCGCTGCGCTTGAACTCAGCGCTGTAACTGTGGGAGCGGGCGTGCCCCCGCAGGGGTTGATCAGCGGTTTTGGACGAATGCGCGAAGGCGCACCGCCGCGACGATGCATTCgg                  >  2:223660/1‑150 (MQ=37)
                                                                                                                                    aaGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGGAAAGCGCTGCGCGTGAAAGCAGCGCAGTAACTGGGGGAGCGGGCGGGCCCGCGAAGGGGTTTATAACCCGTTTTTCAGGGAGGCGCGAATGCGCCCCGCCGCTGCGATGCAATCgg                  >  1:7414048/1‑150 (MQ=255)
                                                                                                                                           ttCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                          >  2:1845452/1‑56 (MQ=255)
                                                                                                                                           ttCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                          <  1:1845452/56‑1 (MQ=255)
                                                                                                                                           ttCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                          <  1:1844949/56‑1 (MQ=255)
                                                                                                                                                gggCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCaa      >  2:5977215/1‑150 (MQ=255)
                                                                                                                                                    gcgcCCGCTCCCACAGGTACTGCGCTGCGCTAGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg  >  2:6414221/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
CCTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAG  >  NC_002947/1784909‑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: