Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,518,231 A→C 7.8% E474A (GAG→GCG) ‡ 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,2310AC 7.8% 181.4 / 7.5 64E474A (GAG→GCG) ‡murEUDP‑N‑acetylmuramoyl‑L‑alanyl‑D‑glutamate‑‑2, 6‑diaminopimelate ligase
Reads supporting (aligned to +/- strand):  ref base A (30/29);  new base C (2/3);  total (32/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 = 7.26e-01

ATTCGCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCA  >  NC_002947/1518083‑1518378
                                                                                                                                                    |                                                                                                                                                   
aTTCGCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGa                                                                                                                                                     <  6:94079/149‑1 (MQ=255)
 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)
                                                                                                                                                    |                                                                                                                                                   
ATTCGCCCCGGTTTCACCAGGCCTGACGACGTCGAGTTCGTCGCCGGCCGTGGCGAAGCCATCGCGCACCTGATTGCCACCGCTGCTGCCAACGATGTGATCGTGCTGGCCGGCAAGGGGCACGAGGATTACCAGGAGATCAATGGCGAGCGCCATGATTTCTCCGATTTGACCGAAGCCGAAAAGGCACTTGCAGCCTGGGAGGCTCCACATGCTTAAGCCGCTGTCACTCAGCCAGCTGACCACAGCACTGAAGGGTCAACACATCGGTGGCGACGCCACGTTCACCGGCGTCA  >  NC_002947/1518083‑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: ‑