Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,518,232 G→T 8.3% E474D (GAG→GAT) ‡ murE → UDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,518,2320GT 8.3% 183.5 / 5.3 60E474D (GAG→GAT) ‡murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base G (26/29);  new base T (2/3);  total (28/32)
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.53e-01

TTCGCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCA  >  NC_002947/1518084‑1518378
                                                                                                                                                    |                                                                                                                                                  
ttCGCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAg                                                                                                                                                    >  3:180553/1‑149 (MQ=255)
   gCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAgcgc                                                                                                                                                 >  7:172376/1‑149 (MQ=255)
    ccccGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAgcgc                                                                                                                                                 >  4:196565/1‑148 (MQ=255)
    ccccGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAgcgc                                                                                                                                                 >  1:162688/1‑148 (MQ=255)
       cGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCAt                                                                                                                                              >  5:206989/1‑148 (MQ=255)
        ggTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGa                                                                                                                                            <  2:175808/149‑1 (MQ=255)
         gTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGAt                                                                                                                                           <  4:266132/149‑1 (MQ=255)
          tttCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGCTCGCCAttgat                                                                                                                                          <  8:141082/149‑5 (MQ=255)
          tttCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGAtt                                                                                                                                          >  1:159857/1‑149 (MQ=255)
          tttCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTACCAACGATGTGATCGTGCTGGCCGCCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGCTCGCCAttgat                                                                                                                                          <  7:132614/149‑5 (MQ=255)
              aCCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTtctc                                                                                                                                      <  5:216905/149‑1 (MQ=255)
               ccAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTcc                                                                                                                                     >  4:167596/1‑149 (MQ=255)
               ccAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTcc                                                                                                                                     <  6:55861/149‑1 (MQ=255)
                   gCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGAtt                                                                                                                                 <  4:180553/149‑1 (MQ=255)
                   gCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGAt                                                                                                                                  <  2:69610/148‑1 (MQ=255)
                   gCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGAt                                                                                                                                  >  1:69610/1‑148 (MQ=255)
                        acgacgTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGAcc                                                                                                                            >  2:174586/1‑149 (MQ=255)
                          gacgTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCg                                                                                                                           >  2:113283/1‑148 (MQ=255)
                            cgTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAg                                                                                                                        <  7:60545/149‑1 (MQ=255)
                             gTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGc                                                                                                                       <  2:159857/149‑1 (MQ=255)
                                   ttCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGaa                                                                                                                   >  4:300445/1‑147 (MQ=255)
                                   ttCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGaa                                                                                                                   >  6:229648/1‑147 (MQ=255)
                                   ttCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGaa                                                                                                                   >  4:199980/1‑147 (MQ=255)
                                         gccggccgTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTACGATTTGACCGAAGCCGAAAAGGCACt                                                                                                           >  8:218908/1‑149 (MQ=255)
                                          ccggccgTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACtt                                                                                                          >  8:97622/1‑149 (MQ=255)
                                           cggccgTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTg                                                                                                         <  6:31734/149‑1 (MQ=255)
                                              ccgTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAg                                                                                                      <  1:174586/149‑1 (MQ=255)
                                               cgTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGc                                                                                                     <  1:113283/149‑1 (MQ=255)
                                                   gCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTgg                                                                                                 <  7:218908/149‑1 (MQ=255)
                                                     gAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTgg                                                                                                 >  1:194871/1‑147 (MQ=255)
                                                                 cACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGc                                                                                   <  2:162688/149‑1 (MQ=255)
                                                                   ccTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCtt                                                                                 <  1:15608/149‑1 (MQ=255)
                                                                       aTTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAg                                                                              <  7:97622/148‑1 (MQ=255)
                                                                              ccGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGt                                                                       >  5:192904/1‑148 (MQ=255)
                                                                               cGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCa                                                                     <  3:167596/149‑1 (MQ=255)
                                                                                gctgctGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCAc                                                                    >  6:162726/1‑149 (MQ=255)
                                                                                         aaCGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTcagccagc                                                           >  5:250867/1‑149 (MQ=255)
                                                                                             aTGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGAc                                                       >  7:111797/1‑149 (MQ=255)
                                                                                                 gATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACca                                                     >  3:186843/1‑147 (MQ=255)
                                                                                                         tGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTg                                             >  3:205339/1‑147 (MQ=255)
                                                                                                         tGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACAACAGCACTg                                             >  7:54580/1‑147 (MQ=255)
                                                                                                           gCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAgg                                         <  5:229648/149‑1 (MQ=255)
                                                                                                             cGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCca                                                                                        <  5:162823/100‑1 (MQ=255)
                                                                                                             cGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCca                                                                                        >  6:162823/1‑100 (MQ=255)
                                                                                                               gCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAgg                                                                                                               >  8:39251/1‑75 (MQ=255)
                                                                                                               gCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAgg                                                                                                               <  7:39251/75‑1 (MQ=255)
                                                                                                               gCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCCGACCACAGCACTGAAGGGTc                                      >  2:136311/1‑148 (MQ=255)
                                                                                                                cAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCt                                                                         >  4:5929/1‑112 (MQ=255)
                                                                                                                cAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCt                                                                         <  3:5929/112‑1 (MQ=255)
                                                                                                                   ggggCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAaca                                  >  5:106486/1‑148 (MQ=255)
                                                                                                                    gggCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACt                                              <  4:1658/135‑1 (MQ=255)
                                                                                                                    gggCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACt                                              >  3:1658/1‑135 (MQ=255)
                                                                                                                    gggCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAAcaca                                <  6:250867/149‑1 (MQ=255)
                                                                                                                               ttACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGAc                     >  3:189208/1‑149 (MQ=255)
                                                                                                                                  ccAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGcc                  <  7:274711/149‑1 (MQ=255)
                                                                                                                                   cAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCa                 <  6:192904/149‑1 (MQ=255)
                                                                                                                                   cAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCa                 <  8:54580/149‑1 (MQ=255)
                                                                                                                                       gaaatcaTGGCGCTCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTgg                         >  4:91753/7‑137 (MQ=255)
                                                                                                                                       gaaatcaTGGCGCTCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTgg                         <  3:91753/131‑1 (MQ=255)
                                                                                                                                       agaTCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGtt             <  4:205339/149‑1 (MQ=255)
                                                                                                                                         aatcaTGGCGCTCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCa           >  6:140467/5‑149 (MQ=255)
                                                                                                                                          tCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCCGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCAc          <  1:136311/149‑1 (MQ=255)
                                                                                                                                                  gAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCa  <  3:199980/149‑1 (MQ=255)
                                                                                                                                                  gAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCa  <  8:172376/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
TTCGCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCA  >  NC_002947/1518084‑1518378

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