Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,137 C→T 10.6% intergenic (+54/+52) 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,1370CT10.6% 290.1 / 17.2 109intergenic (+54/+52)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (58/36);  new base T (3/7);  total (66/43)
Fisher's exact test for biased strand distribution p-value = 8.80e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.92e-01

Only 100 of 118 total aligned reads displayed.

CCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGG  >  NC_002947/1779990‑1780284
                                                                                                                                                   |                                                                                                                                                    
ccTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCtt                                                                                                                                                    <  1:4261689/150‑1 (MQ=255)
 cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTg                                                                                                                                                   <  1:6297523/150‑1 (MQ=255)
 cTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTg                                                                                                                                                   <  2:2243122/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:5384588/150‑1 (MQ=255)
     aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                                              <  2:2127252/150‑1 (MQ=255)
     aTCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAt                                                                                                                                              <  1:5492267/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)
           aGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcgc                                                                                                                                        <  1:8477398/150‑1 (MQ=255)
           aGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcgc                                                                                                                                        <  1:3872786/150‑1 (MQ=255)
            gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGCTTGAAATCAGCGCt                                                                                                                                       >  2:6283493/1‑150 (MQ=255)
              aCCTGGGGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACGGCCCTGCGCGGGAAAGCAgcgcggg                                                                                                                                     >  2:7108400/1‑147 (MQ=255)
                   gtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtg                                                                                                                                <  1:7567638/150‑1 (MQ=255)
                   gtggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtg                                                                                                                                <  1:7567468/150‑1 (MQ=255)
                    tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                               <  1:1051122/150‑1 (MQ=255)
                    tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgt                                                                                                                               <  1:1051085/150‑1 (MQ=255)
                     gggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgtg                                                                                                                              <  1:2394369/150‑1 (MQ=255)
                     gggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCGATCAACCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                             <  1:3252823/150‑1 (MQ=17)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATc                                                                                                                                             >  2:4467743/1‑134 (MQ=255)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATc                                                                                                                                             <  1:4467743/134‑1 (MQ=255)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:4721781/150‑1 (MQ=255)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:5082139/150‑1 (MQ=255)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:1377198/150‑1 (MQ=255)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                             <  1:8181011/150‑1 (MQ=21)
                      ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTgg                                                                                                                             <  2:1215869/150‑1 (MQ=21)
                       gggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTggg                                                                                                                            <  2:5599227/150‑1 (MQ=21)
                        ggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                                                           <  1:7881299/150‑1 (MQ=255)
                         gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             >  2:5348779/1‑147 (MQ=255)
                         gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTgg                                                                                                                             <  1:5348779/147‑1 (MQ=255)
                         gATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCACCAAGGTAAAGCGCTGCGCGTGAAAACAGCCCGGTAACAGTGGGAg                                                                                                                          >  1:4646281/1‑150 (MQ=255)
                            gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAAACAGCGCTGTACCTGTGGGAGCgg                                                                                                                       >  2:5020077/1‑150 (MQ=255)
                            gAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCgg                                                                                                                       <  1:4374424/150‑1 (MQ=21)
                              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)
                                   tgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAAAGCGCTGCGCTTGAAATCAGCGCAGTAACAGTGGGAGAGGGCGTGcc                                                                                                                >  2:306837/1‑150 (MQ=255)
                                         aaGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                                                          >  2:6817188/1‑150 (MQ=255)
                                            tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCCCGAAggg                                                                                                       >  1:6015958/1‑150 (MQ=255)
                                            tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCACTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                       >  2:5529290/1‑150 (MQ=255)
                                            tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACAGCGCTGCCCATGAAATCAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                       >  2:6177522/1‑150 (MQ=38)
                                               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)
                                                                      gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTa                                                                                                                                    <  2:3623385/95‑1 (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                                                                                                                              <  1:1555957/85‑1 (MQ=255)
                                                                                      tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCtgtg                                                                                                                              >  2:1555957/1‑85 (MQ=255)
                                                                                        aaGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAATGAGCGCGGGACCTGTGGGATCCGGCGTGCCCCCCAACGGGGTGTTCAGCGGTTTTGCAGGCAGGCCCGAAGGCGCTCcgccgc                                                           >  2:4632378/1‑150 (MQ=37)
                                                                                             gACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACGGCGCGGCGCGTGGAATAAGCCCTGTACCATGGGGAGCGGGCGTGCCCGCGCAGGGGGTGGTCAGCGGGGTTGCAGGAAGGCGCGTATGCGCGCCGCCGCGTCGa                                                      >  2:1594695/1‑150 (MQ=255)
                                                                                                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:5744320/1‑150 (MQ=21)
                                                                                                                                ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                   >  2:721602/1‑150 (MQ=255)
                                                                                                                                 cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                  >  2:4271672/1‑150 (MQ=255)
                                                                                                                                 cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                  >  2:700117/1‑150 (MQ=255)
                                                                                                                                  acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                 >  1:8587119/1‑150 (MQ=255)
                                                                                                                                  acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                 >  1:1658038/1‑150 (MQ=255)
                                                                                                                                  acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                 >  2:58036/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)
                                                                                                                                        aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                <  2:2766273/129‑1 (MQ=255)
                                                                                                                                        aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                >  1:2766112/1‑129 (MQ=255)
                                                                                                                                        aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                >  1:2766273/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:2127470/1‑150 (MQ=255)
                                                                                                                                           gcgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGc        >  2:2127452/1‑150 (MQ=255)
                                                                                                                                           gcgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGc        >  1:3753354/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)
                                                                                                                                                   |                                                                                                                                                    
CCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGG  >  NC_002947/1779990‑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: