Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,785,120 T→C 6.2% intergenic (+88/+18) 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,1200TC 6.2% 340.8 / 8.7 116intergenic (+88/+18)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (82/21);  new base C (4/3);  total (90/26)
Fisher's exact test for biased strand distribution p-value = 1.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.80e-01

Only 100 of 116 total aligned reads displayed.

CGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGA  >  NC_002947/1784976‑1785266
                                                                                                                                                |                                                                                                                                                  
cGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTGAAATAAGCGCTGTACCTGTGGGAGCGGGCGTGcccc                                                                                                                                               >  1:3188627/1‑149 (MQ=255)
    gcagcggttggagctgtacaacctgcatcgtcagcacctggtgggggctgagggggtgccgcaagtactgggtgagctggagatgagtgaggcgaacctgggtccagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGaa                                                                                                                                           <  2:3921725‑M2/10‑1 (MQ=255)
     cagcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgCGCCCGCGAAg                                                                                                                                          <  1:3320976‑M2/11‑1 (MQ=255)
       tACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAggg                                                                                                                                        >  2:3641592/1‑150 (MQ=255)
           gggTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTg                                                                                                                                    >  2:3757429/1‑150 (MQ=255)
            ttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaccccttcgcgggcaCGCCCgcgaagggccttg                                                                                                                                   <  2:2730231‑M2/18‑6 (MQ=255)
            ttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaccccttcgcgggcaCGCCCgcgaagggccttg                                                                                                                                   <  2:2730293‑M2/18‑6 (MQ=255)
              gagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgtggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATc                                                                                                                                 >  2:4298740‑M2/131‑150 (MQ=255)
                 gCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGc                                                                                                                              >  2:4664975/1‑150 (MQ=255)
                  cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTTCTGCGGTGGGCTTGATATGAGCGCGGTAGCTGTGGGGGCGGGGGTGGCCGCGTAGGGGTTGGTGAGCg                                                                                                                             >  1:2332834/1‑150 (MQ=255)
                  cTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCGTGAAATAAGCGCGGTTCCGGTGGGGGCGGGCGGGCCCGCGAAGGGGGTGTTCAGCg                                                                                                                             >  2:2607961/1‑150 (MQ=37)
                     gATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCATTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGCGGtt                                                                                                                          >  2:2626138/1‑150 (MQ=255)
                          gagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGCTTGAAATCAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCACCGGTGTTGcc                                                                                                                     >  2:323602/1‑149 (MQ=255)
                             tgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTCCTGCGCTGCGCTTGGAATCAGCGCGGTCCCGGGGGGCGCGGGCGGGCCCGCGACGGGGGTGATCAGCGGGTTTGCAGGa                                                                                                                  >  2:5494188/1‑150 (MQ=255)
                               cggCGACCTTGGTTTCCGTCATGTCCACCCTCTCACGGCCCTTCGGGGGCGCGCCCGCCCCCCCCGGTCCGGCGCGGCTATGCAAACCCGCGCTGTACTTGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAg                                                                                                                <  2:3260963/149‑1 (MQ=255)
                                 gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCGACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGgc                                                                                                              >  2:2291344/1‑150 (MQ=255)
                                 gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACGGCGCTGCGCGTGAAATCAGCCCTGTACCTGTGGGAGCGGGCGGGCCCGCGGAGGGGGTGATCAGCGGGTTTGCCGGAAGgc                                                                                                              >  2:675460/1‑150 (MQ=38)
                                  tcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGTCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGgcg                                                                                                             <  2:1544691‑M2/40‑1 (MQ=255)
                                          tggtgggggatgtggtggtgcggcaagtactgtgtgagctggatatgagtgagtcgaagctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGTAAGGCGCGAATgcg                                                                                                     <  1:4524373‑M2/48‑1 (MQ=255)
                                           gTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTTAAATCAGCGCTGTAACTGTGGGCGCGGGCGTGCCCGCGAAGGGGTTGTTCAGCGGGTTTGCAGGAAGGCGCGAATgcgc                                                                                                    >  2:5263736/1‑150 (MQ=38)
                                                 tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCGGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTccgcc                                                                                              >  2:1772601/1‑150 (MQ=255)
                                                 tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACACGTACAGCGCATCGCTTGACAGCCGCGCGGTACCGGTGGGAGCGGGCGTGCCCGCGCAGGGGGTGTTCAGCGGTTTTTCAGGCAGGCGCGAAAGCGCTccccc                                                                                              >  2:1126786/1‑150 (MQ=25)
                                                 tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACAGCGCTGCGCGTTAAATCAGCGCGGTACCAGTGGGGGCGGGCGTGCCCCCGAAGGGGTTGTTCAGCCGTTTTGCAGGAAGGCGCCCAAGCGCAccgcc                                                                                              >  2:1126766/1‑150 (MQ=255)
                                                    aGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACAGCGCTGCGCGTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCGCGTAGGGGGTGTTCAGCGGTTTTGAAGGAAGGCGCGAAAGCGCTCCGCCGCt                                                                                           >  1:4602420/1‑150 (MQ=37)
                                                     gTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCcgccgcgt                                                                                          >  2:1243190/1‑148 (MQ=255)
                                                     gTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACTGCGCTGCGCTTGAAATCAGCGCGGTACCTGTGGGCGCGGGCGTGCCCGCGAAGGGGGTGGTCAGCGGTTTTGCCGGAAGGCGCGCATGCGCGCCGCCGCtg                                                                                          >  2:3541883/1‑149 (MQ=37)
                                                        gACCCTACCACGCCCCTTGGCGGGCGCGCCCGGCCCCACTGGTACGGCGCTGCGCTTTAATTCAGCGCGGTACCTGGGGGCGCGGGGGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGa                                                                                       <  2:556606/150‑1 (MQ=255)
                                                           ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTACTGCGCTGCGCGTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGCAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGc                                                                                    >  2:5296027/1‑150 (MQ=255)
                                                           ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGTAAGGCGCGAATGCGCTCCGCCGCTTCGATGc                                                                                    >  2:4076938/1‑150 (MQ=255)
                                                                 aaGGCCCTTCGCGGGCGCGCCCGCTCCCACCGGGACGGCGCGGCGGCTGCAAACCCCCCGGGTCCCGTGGGTGCGGGCGGGCCCGCGGAGGGGGTGTGGAGCCGTTTTGCAGGCACGCGCGGGGGCGCTCCGCCGGGTAGACGCAGTCgg                                                                              >  2:162110/1‑150 (MQ=255)
                                                                                 gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTAGGCCAGCGGTGCAACCAg                                                              >  1:4062892/1‑150 (MQ=255)
                                                                                 gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTTCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                                              >  1:2156064/1‑150 (MQ=255)
                                                                                  cgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGt                                                             >  2:5587583/1‑150 (MQ=255)
                                                                                   gcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTg                                                            >  2:5162070/1‑150 (MQ=255)
                                                                                     ccGCTCCCACAGGTACTGCGCTGCGCTTGAACTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                                          >  1:1431277/1‑150 (MQ=255)
                                                                                      cGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGACGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCa                                                         >  1:2117344/1‑150 (MQ=255)
                                                                                       tgagtgaggcgaacctgggttcagtcaagtagaccctatccaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCCGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCAt                                                        <  2:5071220‑M2/93‑1 (MQ=255)
                                                                                       tgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCAt                                                        <  2:1817671‑M2/93‑1 (MQ=255)
                                                                                        gagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATg                                                       >  1:289605‑M2/57‑150 (MQ=255)
                                                                                        cTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATg                                                       >  1:4907648/1‑150 (MQ=255)
                                                                                        cTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATg                                                       >  1:4907546/1‑150 (MQ=255)
                                                                                        cTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATg                                                       >  1:1217338/1‑150 (MQ=255)
                                                                                           ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                                                    >  1:5000754/1‑150 (MQ=255)
                                                                                             acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                                                  >  1:543242/1‑150 (MQ=255)
                                                                                             acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                                                  >  1:1692826/1‑150 (MQ=255)
                                                                                             acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                                                  >  1:3273248/1‑150 (MQ=255)
                                                                                                 gTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACgcgc                                              >  2:4949182/1‑150 (MQ=255)
                                                                                                   aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGccc                                            >  2:4623566/1‑150 (MQ=255)
                                                                                                   aCTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGccc                                            >  2:4623962/1‑150 (MQ=255)
                                                                                                    cctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCg                                           <  1:2688765‑M2/106‑1 (MQ=255)
                                                                                                     tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg                                          >  2:4218753/1‑150 (MQ=255)
                                                                                                     tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg                                          >  2:702282/1‑150 (MQ=255)
                                                                                                       cgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCa                                        >  1:2515711/1‑150 (MQ=255)
                                                                                                       cgctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCa                                        >  1:4414728/1‑150 (MQ=255)
                                                                                                        gctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAc                                       >  2:4046235/1‑150 (MQ=255)
                                                                                                        gctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAc                                       >  2:366275/1‑150 (MQ=255)
                                                                                                         ctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAcc                                      >  2:1578817/1‑150 (MQ=255)
                                                                                                          tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc                                     >  2:5542133/1‑150 (MQ=255)
                                                                                                          tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc                                     >  1:3193010/1‑150 (MQ=255)
                                                                                                          tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc                                     >  2:3351373/1‑150 (MQ=255)
                                                                                                              ctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGt                                 >  1:5325709/1‑150 (MQ=255)
                                                                                                               ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt                                >  2:4171952/1‑150 (MQ=255)
                                                                                                               ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt                                >  2:3585154/1‑150 (MQ=255)
                                                                                                               ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATACGCACGCGCCCGGCACCCGGGtt                                >  1:1561839/1‑150 (MQ=255)
                                                                                                                tGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTc                               >  2:184651/1‑150 (MQ=255)
                                                                                                                 gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCa                              >  2:2129806/1‑150 (MQ=255)
                                                                                                                 gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCa                              >  2:2129821/1‑150 (MQ=255)
                                                                                                                 gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCa                              >  2:3410669/1‑150 (MQ=255)
                                                                                                                 gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACACGGGTTCa                              >  2:564944/1‑150 (MQ=255)
                                                                                                                   ttcaaGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                            <  2:4889625/146‑1 (MQ=255)
                                                                                                                   ttcaaGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                            <  2:4889438/146‑1 (MQ=255)
                                                                                                                   aatgcgcggTGTCCCTGTGGGCGGGGGCGGGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                            <  1:1637133/145‑1 (MQ=255)
                                                                                                                   aaTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                            >  2:3723037/1‑150 (MQ=255)
                                                                                                                   aaTCAGCGCGGTCCCTGGGGGCGCGGGGGGGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                            <  1:1637440/150‑1 (MQ=255)
                                                                                                                     tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc                          >  1:264326/1‑150 (MQ=255)
                                                                                                                      cAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCa                         <  1:636081/150‑1 (MQ=255)
                                                                                                                      cAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCa                         <  1:3447987/150‑1 (MQ=255)
                                                                                                                       aGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCag                        >  2:4235289/1‑149 (MQ=255)
                                                                                                                         cctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATcc                      >  1:1894142‑M2/24‑150 (MQ=255)
                                                                                                                          gcTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCa                     >  2:1152864/1‑150 (MQ=255)
                                                                                                                          gcTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCa                     >  1:276213/1‑150 (MQ=255)
                                                                                                                          gcTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCa                     <  1:5162070/150‑1 (MQ=255)
                                                                                                                             gTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACtt                  >  2:4952544/1‑150 (MQ=255)
                                                                                                                                ccTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAc               >  2:3098712/1‑150 (MQ=255)
                                                                                                                                 cggggggcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg              <  2:1332985‑M2/135‑1 (MQ=255)
                                                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg              >  2:3634659/1‑150 (MQ=255)
                                                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg              >  2:4572193/1‑150 (MQ=255)
                                                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg              >  2:3058793/1‑150 (MQ=255)
                                                                                                                                 cTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg              >  1:1972962/1‑150 (MQ=255)
                                                                                                                                  tgtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCAGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcgc             >  2:11774/1‑150 (MQ=255)
                                                                                                                                   cccttcgcgggcaCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcgcg            >  1:1575390‑M2/14‑150 (MQ=255)
                                                                                                                                    ccttcgcgggcaCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGa           >  2:4975552‑M2/13‑150 (MQ=255)
                                                                                                                                     cttcgcgggcaCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGACTGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGAc          >  2:4995664‑M2/12‑150 (MQ=255)
                                                                                                                                     cttcgcgggcaCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGAc          >  1:4872393‑M2/12‑150 (MQ=255)
                                                                                                                                      ggAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACa         >  1:3783116/1‑150 (MQ=255)
                                                                                                                                      ggAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACa         >  1:3783063/1‑150 (MQ=255)
                                                                                                                                           gggggTCCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTAc    <  1:731875‑M2/145‑1 (MQ=255)
                                                                                                                                             gcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGa  <  2:3257400‑M2/147‑1 (MQ=255)
                                                                                                                                             gcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGa  <  2:5373592‑M2/147‑1 (MQ=255)
                                                                                                                                             gcgTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGa  >  1:5357953‑M2/4‑150 (MQ=255)
                                                                                                                                                |                                                                                                                                                  
CGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGA  >  NC_002947/1784976‑1785266

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: