Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,780,108 G→A 18.3% intergenic (+25/+81) PP_1587 → / ← dapC NhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,780,1080GA18.3% 134.8 / 23.1 65intergenic (+25/+81)PP_1587/dapCNhaP‑type Na+(K+)/H+ antiporter/N‑succinyl‑L,L‑diaminopimelate aminotransferase
Reads supporting (aligned to +/- strand):  ref base G (29/24);  new base A (6/6);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

CTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGT  >  NC_002947/1779961‑1780249
                                                                                                                                                   |                                                                                                                                             
cTGAGGGCGCAGCGGCTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCG                                                                                                                                             >  1:35908‑M1/1‑150 (MQ=255)
 tGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGC                                                                                                                                            >  2:280537‑M1/1‑150 (MQ=255)
   aGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCC                                                                                                                                          >  2:472754‑M1/1‑150 (MQ=255)
         cAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCC                                                                                                                                    >  2:1295368‑M1/1‑150 (MQ=255)
                       gTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTActgcgct                                                                                                                      <  1:1295368/150‑1 (MQ=255)
                        tACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAct                                                                                                                           <  2:72302/144‑1 (MQ=255)
                        tACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTAct                                                                                                                           >  1:72302/1‑144 (MQ=255)
                           aaCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACtgcgctgcgc                                                                                                                  >  2:942379/1‑150 (MQ=255)
                           aaCCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatt                                                                                                                  >  2:597744‑M1/1‑140 (MQ=255)
                               tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCAC                                                                                                                                  <  1:1049960‑M1/130‑1 (MQ=255)
                               tGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCAC                                                                                                                                  >  2:1049960‑M1/1‑130 (MQ=255)
                                gCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGGGGTTAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaa                                                                                                             <  1:799729/150‑1 (MQ=255)
                                   tCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGAGAAAGGTGAAAATGAGTGAGGCGAACCGGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATc                                                                                                          >  1:1106445/1‑150 (MQ=255)
                                   cgggtgccgggcgcgtgcgcatggcactggttgcaccgctggccgaatgcatcgaagcggcggagcgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATc                                                                                                          >  2:290300‑M1/101‑150 (MQ=255)
                                       cAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAgcg                                                                                                      >  1:631304/1‑150 (MQ=255)
                                         gCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgca                                                                                                    <  2:35908‑M1/150‑25 (MQ=255)
                                          cACCTGGTGGGGGAGGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGCGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcag                                                                                                   <  1:596990‑M1/150‑26 (MQ=255)
                                           gggcgcgtgcgcatggcactggttgcaccgctggccgaatgcatcgaagcggcggagcgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGt                                                                                                  >  2:33105‑M1/93‑150 (MQ=255)
                                           aCCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgt                                                                                                  >  2:203074‑M2/1‑148 (MQ=255)
                                           aCCTGGTGGGGGATGAGGTGGCGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagt                                                                                                  >  1:1337965‑M1/1‑124 (MQ=255)
                                            ccTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagta                                                                                                 <  2:1337965‑M1/150‑28 (MQ=255)
                                                 tggGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCACGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagtacctgt                                                                                            <  2:152164‑M1/150‑33 (MQ=255)
                                                   ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagtacctgtgg                                                                                          <  1:280537‑M1/150‑35 (MQ=255)
                                                   ggggATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgtacctgtgg                                                                                          >  2:282760‑M2/1‑140 (MQ=255)
                                                      catggcactggttgcaccgctggccgaatgcatcgaagcggcggagcgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                        >  1:876192‑M1/82‑149 (MQ=255)
                                                      catggcactggttgcaccgctggccgaatgcatcgaagcggcggagcgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                        <  2:876192‑M1/68‑1 (MQ=255)
                                                          aGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagtacctgtgggagcggg                                                                                   <  1:597744‑M1/150‑42 (MQ=255)
                                                            gtggtgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgtacctgtgggagcgggcg                                                                                 <  2:369780‑M2/150‑20 (MQ=255)
                                                            gtggtgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgtacctgtgggagcgggcg                                                                                 <  2:1033999‑M2/150‑20 (MQ=255)
                                                               gtgCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGc                                                                              >  1:635171/1‑150 (MQ=255)
                                                                    gCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgtacctgtgg                                                                                          >  1:525005‑M2/1‑123 (MQ=255)
                                                                    gCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgtacctgtgg                                                                                          <  2:525005‑M2/133‑11 (MQ=255)
                                                                      aaGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGaa                                                                       >  1:1059106/1‑150 (MQ=255)
                                                                              ggTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAgggcatga                                                               >  1:82129/1‑145 (MQ=255)
                                                                               gTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaaa                                                                                                            >  1:652255/1‑104 (MQ=255)
                                                                               gTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGaaa                                                                                                            <  2:652255/104‑1 (MQ=255)
                                                                                        gaagcggcggagcgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                        >  1:890418‑M1/48‑115 (MQ=255)
                                                                                        gaagcggcggagcgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGa                                                                                        <  2:890418‑M1/68‑1 (MQ=255)
                                                                                        atatGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGttt                                                     <  1:369780/150‑1 (MQ=255)
                                                                                            gagtgagGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCa                                                 <  1:968289/150‑1 (MQ=255)
                                                                                                   gCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGgc                                          >  2:483439/1‑150 (MQ=255)
                                                                                                    cgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                              >  1:804419‑M1/36‑97 (MQ=255)
                                                                                                    cgcattcgcgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                              <  2:804419‑M1/62‑1 (MQ=255)
                                                                                                     gAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGgcgc                                        <  1:472754/150‑1 (MQ=255)
                                                                                                            cgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGt                                                                                                  <  1:501774‑M1/58‑1 (MQ=255)
                                                                                                            cgccttcctgcaaaaccgctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGt                                                                                                  >  2:501774‑M1/28‑85 (MQ=255)
                                                                                                                  gTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagtac                                                                                                >  1:932822‑M1/1‑53 (MQ=255)
                                                                                                                  gTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACagcgctgatttcaagcgcagcgcagtac                                                                                                <  2:932822‑M1/81‑29 (MQ=255)
                                                                                                                   tCAAGTAGACCCTATCAAGTCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTccgcc                          >  2:1090777/1‑150 (MQ=255)
                                                                                                                   tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGt                                                                  >  1:20046/1‑110 (MQ=255)
                                                                                                                   tCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGt                                                                  <  2:20046/110‑1 (MQ=255)
                                                                                                                          gACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGa                                                               <  2:240986/106‑1 (MQ=255)
                                                                                                                          gACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGa                                                               >  1:240986/1‑106 (MQ=255)
                                                                                                                           aCCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGCGCGGGCGTGCCCGCGACGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGAt                  >  2:117637/1‑150 (MQ=255)
                                                                                                                             gctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                              <  1:552993‑M1/62‑1 (MQ=255)
                                                                                                                             gctgatcaacCCCTTCGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCt                                                                                              >  2:552993‑M1/11‑72 (MQ=255)
                                                                                                                             ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                      >  1:1011235/1‑96 (MQ=255)
                                                                                                                             ccTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAg                                                                      <  2:1011235/96‑1 (MQ=255)
                                                                                                                                  cAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCg           <  1:1090777/150‑1 (MQ=255)
                                                                                                                                   aaGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCgg          >  2:694602/1‑150 (MQ=255)
                                                                                                                                    aaccccttcGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgt                                                                                                  >  2:395502‑M2/4‑59 (MQ=255)
                                                                                                                                    aaccccttcGCGGGCACGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTgt                                                                                                  <  1:395502‑M2/58‑3 (MQ=255)
                                                                                                                                         CTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCg    >  1:483439‑M1/1‑150 (MQ=255)
                                                                                                                                          TTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGGCGGGGAAGGGGTTGATACGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCgg   <  2:338410‑M1/150‑1 (MQ=255)
                                                                                                                                           TCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGcggg  <  2:856737‑M1/150‑2 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
CTGAGGGCGCAGCGGTTGGAGCTGTACAACCTGCATCGTCAGCACCTGGTGGGGGATGAGGTGGTGCGGCAAGTACTGGGTGAGCTGGATATGAGTGAGGCGAACCTGGGTTCAGTCAAGTAGACCCTATCAAGGCCCTTCGCGGGCGCGCCCGCTCCCACAGGTACTGCGCTGCGCTTGAAATCAGCGCTGTACCTGTGGGAGCGGGCGTGCCCGCGAAGGGGTTGATCAGCGGTTTTGCAGGAAGGCGCGAATGCGCTCCGCCGCTTCGATGCATTCGGCCAGCGGT  >  NC_002947/1779961‑1780249

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