Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,731,220 C→T 5.7% intergenic (‑91/+100) dapE ← / ← PP_1526 succinyl‑diaminopimelate desuccinylase/beta‑(1‑3)‑glucosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,731,2200CT 5.7% 210.0 / 2.2 70intergenic (‑91/+100)dapE/PP_1526succinyl‑diaminopimelate desuccinylase/beta‑(1‑3)‑glucosyl transferase
Reads supporting (aligned to +/- strand):  ref base C (33/33);  new base T (2/2);  total (35/35)
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.40e-01

GACGGATCAGGTCGCAGGCCAGTTGAAGGGTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCGGGC  >  NC_002947/1731072‑1731280
                                                                                                                                                     |                                                            
gtggcgatgttgctggtgcagtcgttgccttatgtggcagcgctggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCAc                                                              >  1:2116700‑M1/135‑150 (MQ=255)
gtggcgatgttgctggtgcagtcgttgccttatgtggcagcgctggtgatggcgttcccgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCAc                                                              >  1:1943983‑M1/135‑150 (MQ=255)
   ctaTCAGGTCGCAGGCCAGTTGAAGGGTAGGCGAGAGCTCGGCCCGGGCACGACGCTCTTCCGCTGTGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGccc                                                          <  2:1951200/148‑1 (MQ=255)
     gatgttgctggtgcagtcgttgccttatgtggcagcgctggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCg                                                         <  1:2205340‑M1/21‑1 (MQ=255)
      tCAGGTCGCAGGCCAGTTGAAGGGTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCT                                                       >  2:458827‑M1/1‑150 (MQ=255)
      tCAGGTCGCAGGCCAGTTGAAGGGCGGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGag                                                       >  1:1951200‑M1/1‑148 (MQ=255)
       tgttgctggtgcagtcgttgccttatgtggcagcgctggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCt                                                       >  1:2699322‑M1/128‑150 (MQ=255)
            ctggtgcagtcgttgccttatgtggcagcgctggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCcac                                                  >  1:1275724‑M1/123‑150 (MQ=255)
            cgcgagCCAGTTGAAGGGTAGTCGAGAGCTCGGCATGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGTTGATATCTGAGTGACGTTATCATCTATCAGATCGTTACGAGCAGCACGGCCCATTCGCTGCAACGCGCTCTCCcaca                                                 <  2:1192948/145‑1 (MQ=255)
                 gCCAGTTGAAGGGTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTAc                                            >  1:2328234/1‑150 (MQ=255)
                   cAGTTGAAGGGTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACt                                           >  2:2966906/1‑149 (MQ=255)
                   cAGTTGAAGGGTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACt                                           <  1:2966906/149‑1 (MQ=255)
                        ttgccttatgtggcagcgctggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGcaca                                      <  2:314025‑M1/40‑1 (MQ=255)
                             gTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGGCTTTGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTg                                <  1:2619433/150‑1 (MQ=255)
                              tAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGcac                                       >  2:2596790/1‑142 (MQ=255)
                              tAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGcac                                       <  1:2596790/142‑1 (MQ=255)
                                 gCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcaga                            >  2:668014‑M1/1‑128 (MQ=255)
                                      agCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCtgt                       <  1:212241/150‑1 (MQ=255)
                                       gCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCTTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCtgtg                      <  2:62469/150‑1 (MQ=255)
                                        cTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGc                     <  2:2328234/150‑1 (MQ=255)
                                           tggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCca                                                   <  2:1199723‑M1/27‑1 (MQ=255)
                                           tggtgatggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCca                                                   >  1:1199723‑M1/92‑118 (MQ=255)
                                           ggCAGGGGCCGTCAGGGGGGACGCCGGGGCAAGGCAAGGAGGGGAAATCTGAGAGGCGTTATCTTATATCAAATGGTTTCTGGCAGCACGGCCCATTCGCGGGCACGCCCGGTCCCACAGGTACTGCACAACTCTGAACGCCTGTGCCaa                  <  2:1522199/150‑1 (MQ=255)
                                                gggCCGTCATGGGGGACGCCGGGGCAAGGCAAGGCCGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTCCTGCACAGCTCTGAACGCCTGTGCCAAGCCGg             >  1:286801‑M2/1‑149 (MQ=255)
                                                 ggCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTAGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgt            <  2:2241590‑M2/150‑3 (MQ=255)
                                                 ggCCGTCAGGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgt            <  2:1243293‑M1/150‑39 (MQ=255)
                                                  ggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCGTGcccgcgaat                                                    <  2:666428‑M1/26‑5 (MQ=255)
                                                  ggcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCGTGcccgcgaat                                                    >  1:666428‑M1/85‑106 (MQ=255)
                                                  gCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtg           >  2:1289973‑M2/1‑147 (MQ=255)
                                                  gCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtg           >  2:241808‑M2/1‑147 (MQ=255)
                                                   gcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCca                                                   >  1:2418362‑M1/84‑110 (MQ=255)
                                                   gcgttcctgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCca                                                   <  2:2418362‑M1/27‑1 (MQ=255)
                                                   ccGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgtgg          <  1:668014‑M1/150‑41 (MQ=255)
                                                      gTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtggga        >  1:458827‑M2/1‑144 (MQ=255)
                                                       tCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgt            >  1:1889343‑M2/1‑143 (MQ=255)
                                                       tCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgt            <  2:1889343‑M2/145‑3 (MQ=255)
                                                       tCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTAGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtgggag       >  1:2241590‑M2/1‑143 (MQ=255)
                                                          tgtcgtcgctgcccaagcccgcagaaaaggctgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCgg    >  1:314025‑M1/77‑150 (MQ=255)
                                                          tGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgtgggagcgg    <  1:2469549‑M1/150‑47 (MQ=255)
                                                            ggggACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGaa                              >  2:141748/1‑122 (MQ=255)
                                                            ggggACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGaa                              <  1:141748/122‑1 (MQ=255)
                                                            ggggACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCGTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtgggagcgggc  >  1:147576‑M2/1‑138 (MQ=255)
                                                                             gcaaggcGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgtgggagcgg    >  1:1804971‑M1/1‑85 (MQ=255)
                                                                             gcaaggcGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgtgggagcgg    <  2:1804971‑M1/131‑47 (MQ=255)
                                                                                aggcGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtggga        >  1:1779035‑M2/1‑118 (MQ=255)
                                                                                aggcGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtggga        <  2:1779035‑M2/124‑7 (MQ=255)
                                                                                 ggcGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAg                                                >  1:371311/1‑83 (MQ=255)
                                                                                 ggcGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAg                                                <  2:371311/83‑1 (MQ=255)
                                                                                         tgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcgGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCg     >  2:2289850‑M1/46‑118 (MQ=255)
                                                                                         tgcccaagcgcagcaggcttgaagttgcctgggccgacctgtgcgGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCg     <  1:2289850‑M1/73‑1 (MQ=255)
                                                                                         aaTCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAAc                             <  2:708900/94‑1 (MQ=255)
                                                                                         aaTCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAAc                             >  1:708900/1‑94 (MQ=255)
                                                                                           tCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtgggagc      >  2:1524861‑M2/1‑107 (MQ=255)
                                                                                           tCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGgtgggagc      <  1:1524861‑M2/115‑9 (MQ=255)
                                                                                            ccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCGGGc  >  2:1968832‑M1/43‑118 (MQ=255)
                                                                                            ccaagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCGGGc  <  1:1968832‑M1/76‑1 (MQ=255)
                                                                                            cTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCGTGCCCGCgaa                                                     >  1:1553802‑M1/1‑64 (MQ=255)
                                                                                            cTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCGTGCCCGCgaa                                                     <  2:1553802‑M1/67‑4 (MQ=255)
                                                                                              aagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGa        >  1:1086623‑M1/41‑110 (MQ=255)
                                                                                              aagcgcagcaggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGa        <  2:1086623‑M1/70‑1 (MQ=255)
                                                                                                   gCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCaa                  <  2:1294941/95‑1 (MQ=255)
                                                                                                   gCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCaa                  >  1:1294941/1‑95 (MQ=255)
                                                                                                      caggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAg       <  1:2863446‑M1/71‑1 (MQ=255)
                                                                                                      caggcttgaagttgcctgggccgacctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAg       >  2:2863446‑M1/33‑103 (MQ=255)
                                                                                                        aTCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGaa                              <  1:2147630/78‑1 (MQ=255)
                                                                                                        aTCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGaa                              >  2:2147630/1‑78 (MQ=255)
                                                                                                                   aaTGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgtgg          >  2:100878‑M1/1‑47 (MQ=255)
                                                                                                                   aaTGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACcggcttggcacaggcgttcagagctgtgcagtacctgtgg          <  1:100878‑M1/87‑41 (MQ=255)
                                                                                                                   aaTGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCg              <  2:1001579/83‑1 (MQ=255)
                                                                                                                   aaTGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCg              >  1:1001579/1‑83 (MQ=255)
                                                                                                                               ctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAg       >  1:2000014‑M1/8‑78 (MQ=255)
                                                                                                                               ctgtgcaGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAg       <  2:2000014‑M1/71‑1 (MQ=255)
                                                                                                                                                     |                                                            
GACGGATCAGGTCGCAGGCCAGTTGAAGGGTAGGCGAGAGCTCGGCAGGGGCCGTCATGGGGGACTCCGGGGCAAGGCAAGGCGGGGAAATCTGAGGGGCGTTATCTTATATCAAATGGTTTCGGGCAGCACGGCCCATTCGCGGGCACGCCCGCTCCCACAGGTACTGCACAGCTCTGAACGCCTGTGCCAAGCCGGTGGGAGCGGGC  >  NC_002947/1731072‑1731280

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: