Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,135 C→A 8.6% intergenic (+52/+54) 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,780,1350CA 8.6% 309.8 / 17.5 108intergenic (+52/+54)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (60/36);  new base A (2/7);  total (65/43)
Fisher's exact test for biased strand distribution p-value = 3.02e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

Only 100 of 120 total aligned reads displayed.

ACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGG  >  NC_002947/1779986‑1780284
                                                                                                                                                     |                                                                                                                                                      
aCAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgctgc                                                                                                                                                        <  2:1475401/150‑1 (MQ=255)
 cAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTActgggcggcg                                                                                                                                                       >  1:1216579/1‑150 (MQ=255)
    ccTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCtt                                                                                                                                                    <  1:4261689/150‑1 (MQ=255)
     cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTg                                                                                                                                                   <  1:6297523/150‑1 (MQ=255)
      tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTTGATATGAGTGAGGCTAACCTGGGTTCAGTCAAGTAGACCATATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCTCTTCGCTTGa                                                                                                                                                  <  2:7882520/150‑1 (MQ=255)
      tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGa                                                                                                                                                  >  2:7779558/1‑150 (MQ=255)
       gCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCTGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaa                                                                                                                                                 <  1:1580818/150‑1 (MQ=255)
       gCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCACGTAAAGCGCTGCGCGTGaa                                                                                                                                                 >  1:2590039/1‑150 (MQ=255)
         aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                                              <  2:2127252/150‑1 (MQ=255)
         aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                                              <  2:2126687/150‑1 (MQ=255)
         aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                                              <  1:6612096/150‑1 (MQ=255)
           cGTCAGCACGTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCa                                                                                                                                            <  1:4152292/150‑1 (MQ=255)
             tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTTCGCTGCGCGTGTAATAAgc                                                                                                                                          >  1:8798294/1‑150 (MQ=255)
             tCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCCCTGCGGGGGCAATCAgc                                                                                                                                          >  1:5901363/1‑150 (MQ=255)
              cAGCACGTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                                                         <  1:3520604/150‑1 (MQ=255)
              cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                                                         <  2:1080250/150‑1 (MQ=255)
              cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                                                         <  1:2704526/150‑1 (MQ=255)
               aGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcgc                                                                                                                                        <  1:8477398/150‑1 (MQ=255)
                gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGCTTGAAATCAGCGCt                                                                                                                                       >  2:6283493/1‑150 (MQ=255)
                  aCCTGGGGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACGGCCCTGCGCGGGAAAGCAgcgcggg                                                                                                                                     >  2:7108400/1‑147 (MQ=255)
                    cTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCGCCCCCAAGTAATGCGCGGCGCGTGTAAACAGCGCGGTAc                                                                                                                                   >  1:1596161/1‑150 (MQ=255)
                       gtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtg                                                                                                                                <  1:7567638/150‑1 (MQ=255)
                        tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                               <  1:1051122/150‑1 (MQ=255)
                        tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                               <  1:1051085/150‑1 (MQ=255)
                          ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATc                                                                                                                                             <  1:4467743/134‑1 (MQ=255)
                          ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:1377198/150‑1 (MQ=255)
                          ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:5082139/150‑1 (MQ=255)
                          ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:3683955/150‑1 (MQ=255)
                          ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                             <  2:1215869/150‑1 (MQ=21)
                          ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                             <  1:8181011/150‑1 (MQ=21)
                           gggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTggg                                                                                                                            <  2:5599227/150‑1 (MQ=21)
                            ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                                                           <  1:7881299/150‑1 (MQ=255)
                            ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGGACTGCGCTGCGCGGGAAAACAGCCCGGGACCGGTgggg                                                                                                                           >  1:5463551/1‑149 (MQ=255)
                             gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:5348779/147‑1 (MQ=255)
                             gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             >  2:5348779/1‑147 (MQ=255)
                             gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCACCAAGGTAAAGCGCTGCGCGTGAAAACAGCCCGGTAACAGTGGGAg                                                                                                                          >  1:4646281/1‑150 (MQ=255)
                                gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAAACAGCGCTGTACCTGTGGGAGCgg                                                                                                                       >  2:5020077/1‑150 (MQ=255)
                                  ggtggtGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCGCCCACCCGGAATGCGCGGCGGGTGTAATAAGCCCGGGAACAGTGGGGGCgggg                                                                                                                     >  1:3800269/1‑149 (MQ=255)
                                   gtggtgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTGAAATCAGCGCGGTACCAGTGGGAGCGGGCg                                                                                                                    >  1:2202605/1‑150 (MQ=255)
                                    tggtgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGcac                                                                                                                   <  1:1967901/150‑3 (MQ=255)
                                    tggtgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCGGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGcac                                                                                                                   <  1:6590534/150‑3 (MQ=255)
                                                tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCCCGAAggg                                                                                                       >  1:6015958/1‑150 (MQ=255)
                                                tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCACTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                       >  2:5529290/1‑150 (MQ=255)
                                                tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACAGCGCTGCCCATGAAATCAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                       >  2:6177522/1‑150 (MQ=38)
                                                   tGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCa                                                                                                                                            >  2:5137558/1‑110 (MQ=255)
                                                   tGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCa                                                                                                                                            <  1:5137558/110‑1 (MQ=255)
                                                          gCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGCGCGGGCGTGCCCGCGAAGGGGGTGATCAGc                                                                                             >  1:7592934/1‑150 (MQ=255)
                                                          gCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGc                                                                                             >  1:7592952/1‑150 (MQ=255)
                                                            tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGCAGGGGGTGATCCGCgg                                                                                           >  2:6524205/1‑150 (MQ=255)
                                                                tatGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCGGTTCCGGTGGGAGCGGGCGGGCCCGCGCAGGGGTTGATCAGCGGgttt                                                                                       >  2:4609652/1‑150 (MQ=38)
                                                                          gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                                                    >  1:3623385/1‑95 (MQ=255)
                                                                          gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTACCACAGGTACAGCGCTGCGCTTGAAATAAGCGCAGAACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGAAGGAAGgc                                                                             >  2:4237873/1‑150 (MQ=38)
                                                                          gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTACTGGGCGGCGGGTGGAAACAGCGCGGTGCCGGGGGGCGCGGGCGTGCCCCCGGAGGGGGTGGTCAGCGGGTTTTGAGGTAGgc                                                                             >  2:25797/1‑150 (MQ=255)
                                                                          gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCGCCCACCGGTACACCGCGGCGCGTGAAAGCAGCGCGGTACCGGGGGGCGCGGGCGGGCCCGCCCCCGGGGTGTTCAGCGGTTTTGCCGGAACgc                                                                             >  2:25788/1‑150 (MQ=255)
                                                                           cGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTAATGCGCTGCGCTGGTAATCAGCGCGGTACCGGGGGGAGCGGGCGTGCCCCCGGAGGGGTTGATAAGCGGTTTTGCAGGAAGgcg                                                                            >  1:4405917/1‑150 (MQ=38)
                                                                           cGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGGTCTGGGCTGGGCTTGGAATGAGCGCGGTACCAGTGGGAGCGGGCGTGCCCCCGCAGGGGGTGATGAGCCGTTTTGGGGGGAGgcg                                                                            >  1:4406195/1‑150 (MQ=255)
                                                                            gAAGCTGGGTTCAGTCAAGTAGACGCTATCAAGGCCCTTCGCGGGGGCGCCCGCTGGCACAGGTACTGCGGTGCGGTTGAAATGAGCGCTGTACCTGTGGGAGCGGGCGTGGCCGCGAAGGGGTTGTTCAGCGGTTTTGCAGGGAGgcgg                                                                           >  2:325877/1‑149 (MQ=255)
                                                                            gAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGGCCGCGGAGGGGTTGATCAGCGGTTTTGCAGGGAGgcgc                                                                           >  2:3964349/1‑150 (MQ=38)
                                                                            gAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACGGCGCTGCGCGTGAAATCAGCGCGGTACCGGGGGGCGCGGGCGGGCCCGCGCACGGGGGGATCCGCGGTTTTGCACGAAGgcgc                                                                           >  2:3964798/1‑150 (MQ=37)
                                                                                cTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAAt                                                                       >  2:2883991/1‑150 (MQ=255)
                                                                                     ttCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGTACAGCGCTGCGCGGGAAAACAGCGCGGTACCTGTGGGGGCGGGCGGGCCCGCGAAGGGGTTGTTCAGCCGTTTTGCAGGAAGGCGCCAATgcgcg                                                                  >  2:1745185/1‑149 (MQ=255)
                                                                                      tCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTACTGCGCTGCGCGTGCCATCCGCGCGGTACCTGTGGGAGCGGGCGTGCCCGCGACGGGGGTGATCACCGGCGTTGCAGGAAGGCGCGAAGGCGCTc                                                                 >  2:7074332/1‑150 (MQ=25)
                                                                                      tCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCGGGGGGAGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGACGGCGCGCATGCGCTc                                                                 >  2:7073964/1‑150 (MQ=38)
                                                                                          tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgtg                                                                                                                              >  2:1555957/1‑85 (MQ=255)
                                                                                            aaGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAATGAGCGCGGGACCTGTGGGATCCGGCGTGCCCCCCAACGGGGTGTTCAGCGGTTTTGCAGGCAGGCCCGAAGGCGCTCcgccgc                                                           >  2:4632378/1‑150 (MQ=37)
                                                                                              gTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGTAATAAGCGCTGTACCGGTGGGAGCCGGCGTGCCCGCGCAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGCATGCGCTCCGCCGCtt                                                         >  2:3265443/1‑150 (MQ=39)
                                                                                                    ccTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCCCAGGTACGGCGCGGCGCGTGTAAACACCGCGGTACCGGTGGGAGCGGGCGGGCCCCCGCAGGGGTTTGTAAGCCGGTTTTCAGGCAGGCGCGAAAGCGCACCCCCGCGTCGCTGc                                                   >  2:3077475/1‑150 (MQ=255)
                                                                                                    ccTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACACGGACTGCGCTGGGGGGGCAATCAGAGCAGTACCGGGGGGAGCGGGCGTGCCCGCGCAGGGGTTGGTCAGCGGTGTTGCAGGAACGCGCGAAGGCGCTCCGCCGCGGCGATGc                                                   >  2:3367330/1‑150 (MQ=255)
                                                                                                                    gcgGGCACGCCCGCTCCCACAGGTAAAGCGCTGGTTTAAAGCGCAGCGCCCTACCGGGGGGCCCGGGCCCGCCCGCCCAGGGGGTGATCAGCGGGTTTGCAGGCAGGCGCGAATGCGCTCCGCCGCGTCGATGCAGTCGGCCCGccgggc                                   >  2:8264164/1‑145 (MQ=9)
                                                                                                                     cgGGCGCGCCCCCTCCCACAGGTACAGCGCTGCGCTTGAAATCAGCGCTGTACCGGTGGGATCCGGCGTGCCCCCGAACGGGTTGATCATCGGTTTTGCAGGAAGGCGCGAATGCGCACCGCCCCGTCCATGCATTCGGCCAGCGGTGCa                                  >  2:627128/1‑150 (MQ=255)
                                                                                                                            gcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGTGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTg                           >  2:7993439/1‑150 (MQ=255)
                                                                                                                               cGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCa                        >  1:2690311/1‑150 (MQ=255)
                                                                                                                                  tCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgc                     >  2:5744341/1‑150 (MQ=21)
                                                                                                                                  tCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgc                     >  2:5744320/1‑150 (MQ=21)
                                                                                                                                    ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                   >  2:721602/1‑150 (MQ=255)
                                                                                                                                     cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                  >  2:4271654/1‑150 (MQ=255)
                                                                                                                                     cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                  >  2:4271672/1‑150 (MQ=255)
                                                                                                                                     cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                  >  2:700117/1‑150 (MQ=255)
                                                                                                                                      acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                 >  1:1658038/1‑150 (MQ=255)
                                                                                                                                      acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                 >  1:8587119/1‑150 (MQ=255)
                                                                                                                                       caGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcg                >  1:8779212/1‑150 (MQ=255)
                                                                                                                                        aGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgc               >  1:1373050/1‑150 (MQ=255)
                                                                                                                                        aGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgc               >  2:8343945/1‑150 (MQ=255)
                                                                                                                                          gTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACgcgc             >  1:4493678/1‑150 (MQ=255)
                                                                                                                                           tACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGcc            >  2:6538831/1‑150 (MQ=255)
                                                                                                                                            aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                <  2:2766273/129‑1 (MQ=255)
                                                                                                                                            aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                <  2:2766112/129‑1 (MQ=255)
                                                                                                                                            aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                >  1:2766273/1‑129 (MQ=255)
                                                                                                                                            aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                >  1:2766112/1‑129 (MQ=255)
                                                                                                                                            aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGccc           >  1:4906696/1‑150 (MQ=255)
                                                                                                                                             ctgcgctgcgCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCg          >  2:1293962/1‑150 (MQ=255)
                                                                                                                                              tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg         >  2:354317/1‑150 (MQ=255)
                                                                                                                                               gcgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGc        >  2:2127452/1‑150 (MQ=255)
                                                                                                                                               gcgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGc        >  1:3753354/1‑150 (MQ=255)
                                                                                                                                               gcgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGc        >  2:2127470/1‑150 (MQ=255)
                                                                                                                                                   tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc    >  2:5504698/1‑150 (MQ=255)
                                                                                                                                                   tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc    >  1:4325201/1‑150 (MQ=255)
                                                                                                                                                   tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc    >  1:8927225/1‑150 (MQ=255)
                                                                                                                                                     cgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCgg  <  1:4271654/150‑1 (MQ=255)
                                                                                                                                                     cgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCgg  <  1:4271672/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
ACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGG  >  NC_002947/1779986‑1780284

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