Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,171 T→C 5.3% 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,780,1710TC 5.3% 291.8 / 4.6 95intergenic (+88/+18)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base T (75/11);  new base C (3/2);  total (81/14)
Fisher's exact test for biased strand distribution p-value = 1.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

TGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGACCC  >  NC_002947/1780022‑1780320
                                                                                                                                                     |                                                                                                                                                     
cctgagggcgcagcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgt                                                                                                                                                       <  1:1678621‑M2/1‑1 (MQ=255)
cctgagggcgcagcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggccctacgcgggcgt                                                                                                                                                       >  2:2767740‑M2/150‑150 (MQ=255)
       gCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACGGCGCTGCGCTTGAAAACAGCGCTGTAACAGGGGGAGCGGGCGTGCCcgcg                                                                                                                                                >  1:4445198/1‑150 (MQ=255)
        cAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGa                                                                                                                                               >  2:1052259/1‑150 (MQ=255)
        cAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCGCTGCGCTTGAAATAAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGa                                                                                                                                               >  2:1462889/1‑150 (MQ=255)
           agcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAgg                                                                                                                                            <  1:3515743‑M2/12‑1 (MQ=255)
           agcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgCGCCCGCGAAgg                                                                                                                                            <  1:189611‑M2/12‑1 (MQ=255)
            gcggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgCGCCCGCGAAggg                                                                                                                                           <  1:3461548‑M2/13‑1 (MQ=255)
             aCTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACAGCCCTGCGCTTGAAAACAGCGCTGTAACAGTGGGAGCGGGCGTTCCCGCGCGgggg                                                                                                                                          >  2:1821883/1‑150 (MQ=38)
             aCTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAATGCGCTGCGCTTGCAATCAGCGCTGTAACAGTGGGCGCGGGCGGGCCCGCGCAgggg                                                                                                                                          >  1:2833932/1‑150 (MQ=255)
              ggttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGt                                                                                                                                         <  1:1480257‑M2/15‑1 (MQ=255)
                ttggagctgtacaacctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTg                                                                                                                                       <  1:750586‑M2/17‑1 (MQ=255)
                     aGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAg                                                                                                                                  >  2:1425459/1‑150 (MQ=255)
                       cTGCATAGGAGTGCGGCCACCCGGGTTTCTGGCACTTAGCCCATAACAAGGCCCTGGGGGGGCGCGCCCGCCCCCACCGGTACGGGGGGGCGCGCAAAAACAGCCGTGAACGGGTGGCCGCGGGCGTGCCCGCGAAGGGGTTGATCAGCg                                                                                                                                <  2:3689921/150‑1 (MQ=255)
                        tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTACTGCGCTGCGCTTGAAATCGGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGACGGGGTTGATCAGCgg                                                                                                                               >  2:3793500/1‑150 (MQ=255)
                        tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACATCGCTGCGCGTGAAATCAGCGCGGTAACAGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATAAGCgg                                                                                                                               >  2:3793355/1‑150 (MQ=37)
                        tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCGCCCACAGGTACAGCGCTGCGCGTGGAATCAGCGCGGTTCCTGGGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCgg                                                                                                                               >  2:3793287/1‑150 (MQ=38)
                        tGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCCCTCCCCCAGGTACAGCGCGGCGCGTGGAATCAGCGCTGTACCTGTGGGCGCGGGCGTGCCCCCGAACGGGTTGATCAGCgg                                                                                                                               >  1:2901813/1‑150 (MQ=255)
                               gagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGGACTGCGCTGCGCATGAAAACAGCGCGGTACCTGTGGGCGCGGGCGTGCCCGCGCAGGGGGGGATCAGCGGTGTTGcc                                                                                                                        >  2:2730181/1‑149 (MQ=37)
                               ctgcatcgtcagcacctggtgggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCa                                                                                                                        <  1:384322‑M2/32‑1 (MQ=255)
                                 gtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCACGTACGGCGCTGCGCGTGGAATCCGCGCGGTACCAGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGGTAAGCCGGGTTGCAgg                                                                                                                      >  2:38338/1‑150 (MQ=25)
                                  tgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCGTGAAAACAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGa                                                                                                                     >  2:2546919/1‑150 (MQ=255)
                                      gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCGGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATGAGCGGTTTTTCAGGAAGgc                                                                                                                 >  2:3629445/1‑150 (MQ=255)
                                                   gggggatgaggtggtgcggcaagtactgggtgagctggatatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTcc                                                                                                    >  2:2897750‑M2/99‑150 (MQ=255)
                                                      tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGTGCGCAAgcgcacggaa                                                                                                 >  1:3754045/1‑144 (MQ=38)
                                                      tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAAACAGCGCTGTAACTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGTGCGCAAgcgcacggaa                                                                                                 >  1:3754064/1‑144 (MQ=37)
                                                           tAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCAGTGGGAGCGGGCGTGCCCGCGAAGGGGGTGGTCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCGTc                                                                                            >  2:2132875/1‑150 (MQ=255)
                                                              aCCATACCACGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACGGCGCTGCGATTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGAt                                                                                         <  2:371697/150‑1 (MQ=37)
                                                                ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACCTGTACAGCGCTGGGCGTGAAATCAGCGCTGTACCTGTGGGGGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGGAGGAAGGCGCGCATTCGCTCCGCCGCGTCGATGc                                                                                       >  1:4038446/1‑150 (MQ=37)
                                                                ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGGGTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCGTCGATGc                                                                                       >  2:1844971/1‑150 (MQ=255)
                                                                ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGTAGGCGCGAATGCGCTCCGCCGCTTCGATGc                                                                                       >  2:3013047/1‑150 (MQ=255)
                                                                ccTATCAAGGCCCTTCGCGGGCGCGCCCGCACCCACCGGTACTGCGCTGCGCGTGAAAACAGCGCTGTACCAGGGGGAGCGGGCGGGCCCGCGAAGGGGGTGATCAGCGGGTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCGTCCATGc                                                                                       >  1:4123899/1‑150 (MQ=255)
                                                                  tATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCCCAGGTACTGCGCTGCGCTTGAATTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAATGGGTTGATCAGCGGTTTTGCCGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCAt                                                                                     >  2:1042902/1‑150 (MQ=255)
                                                                        ggCCCTTCGCGGGCGCGCCCGCTCCCACACGTAAATCGCTGCGCTTTAAATCAGCGCTGTAACAGTGGGAGCGGGCGTGCCCGCGAAGGGGTTTATAAGCGGTTTTGCAGGAAGGCGCGAATGCGCACCGCCGCTTAGATGCATTCGGcc                                                                               >  2:3223195/1‑150 (MQ=38)
                                                                        accccTTCGCGGGCGCGCCCCCGCCCACACGTACAGCGCGGCGGGTGCACACAGCGCGGTACCGGTGGGGGCGGGCGGGCCCCCGAAGGGGGTGGTCAGCCGGTTTTCAGGCAGGCGCGAATGCGCTCCGCCCCGTCGAGGCATTAGGcc                                                                               >  2:480853/3‑150 (MQ=255)
                                                                                   ggCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGGGGGAGCGGGCGTGCCCGCGAAGGGGGTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCCGCGGTGCAAc                                                                    >  2:2519231/1‑150 (MQ=255)
                                                                                   agctggatatgagtgaggcgaccctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAAc                                                                    <  1:212031‑M2/84‑1 (MQ=255)
                                                                                      gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGCAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGCATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                                                 >  2:4194570/1‑150 (MQ=255)
                                                                                      gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAg                                                                 >  2:4194591/1‑150 (MQ=255)
                                                                                      gcgcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAAGAAGGCGCGAATGCGCACCCCCGCTTCGATGCATTAGGCCAGCGGTGCAACCAg                                                                 >  1:338552/1‑150 (MQ=255)
                                                                                        gcCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTg                                                               >  2:3566007/1‑150 (MQ=255)
                                                                                          tatgagtgaggcgaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                                             <  1:2105463‑M2/91‑1 (MQ=255)
                                                                                          ccGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGc                                                                                                                         <  1:3738858/90‑1 (MQ=255)
                                                                                          ccGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGc                                                                                                                         >  2:3738858/1‑90 (MQ=255)
                                                                                          ccGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                                             >  1:2827979/1‑150 (MQ=255)
                                                                                          ccGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                                             >  1:3069978/1‑150 (MQ=18)
                                                                                          ccGCTCCCACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGcc                                                             >  1:2333873/1‑150 (MQ=18)
                                                                                           cGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCa                                                            >  1:3953307/1‑150 (MQ=255)
                                                                                           cGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCa                                                            >  2:264502/1‑150 (MQ=255)
                                                                                             cTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATg                                                          >  2:2813817/1‑150 (MQ=255)
                                                                                               cccACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcg                                                        >  1:2448099/1‑150 (MQ=21)
                                                                                               cccACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcg                                                        >  1:2448115/1‑150 (MQ=21)
                                                                                               cccACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcg                                                        >  1:4234389/1‑150 (MQ=21)
                                                                                               cccACAGGTACAGCGCTGATTTCAAGCGCAGCGCAGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcg                                                        >  1:2448087/1‑150 (MQ=21)
                                                                                                ccACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATgcgc                                                       >  1:1623867/1‑150 (MQ=255)
                                                                                                 cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                                                      >  1:4435471/1‑150 (MQ=255)
                                                                                                 cacaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCa                                                      >  1:4378038/1‑150 (MQ=255)
                                                                                                  acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                                                     >  2:2904521/1‑150 (MQ=255)
                                                                                                  acaGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAc                                                     >  1:4318149/1‑150 (MQ=255)
                                                                                                    aGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgc                                                   >  2:168418/1‑150 (MQ=255)
                                                                                                    aGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgc                                                   >  2:2827138/1‑150 (MQ=255)
                                                                                                     ggTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCAcgcg                                                  >  1:1227757/1‑150 (MQ=255)
                                                                                                      gaacctgggttcagtcaagtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACgcgc                                                 >  2:2005091‑M2/48‑150 (MQ=255)
                                                                                                       tACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGcc                                                >  1:1365628/1‑150 (MQ=255)
                                                                                                          tgcgctgcgcTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCgg                                             >  2:105885/1‑150 (MQ=255)
                                                                                                             gctgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAc                                          >  2:4050564/1‑150 (MQ=255)
                                                                                                               tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc                                        >  1:261342/1‑150 (MQ=255)
                                                                                                               tgcgctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCAccc                                        >  1:3110387/1‑150 (MQ=255)
                                                                                                                  gctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCggg                                     >  1:3595862/1‑150 (MQ=255)
                                                                                                                   ctTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGt                                    >  1:3381216/1‑150 (MQ=255)
                                                                                                                    ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt                                   >  2:865504/1‑150 (MQ=255)
                                                                                                                    ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt                                   >  1:3820697/1‑150 (MQ=255)
                                                                                                                    ttGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGtt                                   >  1:2613964/1‑150 (MQ=255)
                                                                                                                      gAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCa                                 >  1:3703866/1‑150 (MQ=255)
                                                                                                                       aaaTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAc                                >  2:2088198/1‑150 (MQ=255)
                                                                                                                        gtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                               <  1:761151‑M2/121‑1 (MQ=255)
                                                                                                                        gtagaccctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                               <  1:4102935‑M2/121‑1 (MQ=255)
                                                                                                                        aaTCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcc                               >  1:1325150/1‑150 (MQ=255)
                                                                                                                          tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGACGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc                             >  2:522044/1‑150 (MQ=255)
                                                                                                                          tCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCAcccc                             >  2:521930/1‑150 (MQ=255)
                                                                                                                            aGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCAt                           >  2:3389302/1‑150 (MQ=255)
                                                                                                                              cctatcaaggcccttcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATcc                         >  1:80683‑M2/24‑150 (MQ=255)
                                                                                                                                cTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCAc                       >  1:1887899/1‑150 (MQ=255)
                                                                                                                                   tACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTc                    >  2:2475729/1‑150 (MQ=255)
                                                                                                                                    aaggcccgtcgcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCa                   <  1:2498388‑M2/133‑1 (MQ=255)
                                                                                                                                      accccttcgcgggcaCGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcg                 >  2:2790604‑M2/16‑150 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcgc                >  1:3760743/1‑150 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcgc                >  1:3756947/1‑150 (MQ=255)
                                                                                                                                       tgtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCAGCACCCGGGTTCACCCCATCCACTTCAcgc                >  2:2905385/1‑150 (MQ=255)
                                                                                                                                        gtgGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCAcgcg               >  1:915280/1‑150 (MQ=255)
                                                                                                                                           ggAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACa            >  1:2648792/1‑150 (MQ=255)
                                                                                                                                           ggAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACa            >  2:1496896/1‑150 (MQ=255)
                                                                                                                                            gAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAg           >  1:1293581/1‑150 (MQ=255)
                                                                                                                                            gAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAg           >  2:3721892/1‑150 (MQ=255)
                                                                                                                                            gAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAg           >  1:2368190/1‑150 (MQ=255)
                                                                                                                                              gcgggcgTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGt         >  2:566567‑M2/8‑150 (MQ=255)
                                                                                                                                                    gTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGAcc   >  2:2324698‑M2/2‑150 (MQ=255)
                                                                                                                                                     tGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGAccc  >  2:918181/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGTGCAACCAGTGCCATGCGCACGCGCCCGGCACCCGGGTTCACCCCATCCACTTCACGCGACAGGTACGACCC  >  NC_002947/1780022‑1780320

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: