Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,246 T→G 10.7% F168C (TTC→TGC) ‡ murD → UDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,521,2460TG10.7% 159.9 / 9.8 56F168C (TTC→TGC) ‡murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base T (26/24);  new base G (3/3);  total (29/27)
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

CAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGT  >  NC_002947/1521100‑1521390
                                                                                                                                                   |                                                                                                                                                
cAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttcc                                                                                                                                                <  6:193816/149‑1 (MQ=255)
    gCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTtccagcttccagct                                                                                                                                            >  6:120618/1‑149 (MQ=255)
     ccAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTg                                                                                                                                           <  1:243970/149‑1 (MQ=255)
      cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCACGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTACCCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCCGCTg                                                                                                                                           >  1:60917/1‑149 (MQ=255)
      cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTgg                                                                                                                                          <  5:303769/149‑1 (MQ=255)
             aCCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaagctc                                                                                                                                   <  3:83781/149‑5 (MQ=255)
              ccGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                  <  2:320180/149‑1 (MQ=255)
               cGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTACCCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCAAGCTGGAAACCACt                                                                                                                                  >  8:146981/1‑148 (MQ=255)
                         cTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCaa                                                                                                                       <  4:19407/149‑1 (MQ=255)
                          tGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAAc                                                                                                                      >  2:375870/1‑149 (MQ=255)
                                gTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGaa                                                                                                                <  6:251016/149‑1 (MQ=255)
                                  gAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGt                                                                                                              >  1:56767/1‑149 (MQ=255)
                                  gAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGt                                                                                                              >  5:105666/1‑149 (MQ=255)
                                   aaaTGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGt                                                                                                              <  6:294607/148‑1 (MQ=255)
                                    aaTGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                            <  1:375870/149‑1 (MQ=255)
                                     aTGGC‑CGCCTAGTCCGTCAAACGCGTGGCCGTGGGCGGTAACCTCTGTACCCCTGCACTCGACCTGCTGGCCGATGACATCTAGCTGTAGTTGCTGTAGCTTTCCAGCTTCCAGCTGGAAACAACTGACCAGCTCAACGCGGATGTTGc                                                                                                           <  4:243243/149‑1 (MQ=255)
                                        gc‑cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGa                                                                                                         <  1:366242/148‑1 (MQ=255)
                                         c‑cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGAcc                                                                                                       >  5:49167/1‑149 (MQ=255)
                                           cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCg                                                                                                      >  4:256026/1‑149 (MQ=255)
                                           cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCg                                                                                                      >  3:358290/1‑149 (MQ=255)
                                            gccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACTCCGAAGTGGCGACCGt                                                                                                     >  5:298604/1‑149 (MQ=255)
                                               aaGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCt                                                                                                  <  3:256026/149‑1 (MQ=255)
                                                 ggCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCAATCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCCGa                                                                                                <  1:73498/149‑1 (MQ=255)
                                                  gCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGaa                                                                                               <  3:330156/149‑1 (MQ=255)
                                                   ccGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                              <  6:105666/149‑1 (MQ=255)
                                                    cGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACa                                                                                             >  1:324351/1‑149 (MQ=255)
                                                       cAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCa                                                                                          <  3:168111/149‑1 (MQ=255)
                                                          aCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCg                                                                                       <  8:159624/149‑1 (MQ=255)
                                                            gcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCg                                                                                       >  7:242172/1‑147 (MQ=255)
                                                               tGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAc                                                                                  <  7:16055/149‑1 (MQ=255)
                                                                   cGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACcaca                                                                              >  8:127037/1‑149 (MQ=255)
                                                                                cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGGTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                              <  1:119427/120‑1 (MQ=255)
                                                                                cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                              >  2:119427/1‑120 (MQ=255)
                                                                                cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAg                                                                 >  5:286834/1‑149 (MQ=255)
                                                                                   ggTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCg                                                               <  1:174802/148‑1 (MQ=255)
                                                                                    gTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                             <  7:146981/149‑1 (MQ=255)
                                                                                          cGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCg                                                       >  5:32007/1‑149 (MQ=255)
                                                                                                        ctgGCCGATGACATCGAGCTGTACGTTCTTGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTg                                           <  4:233260/147‑1 (MQ=255)
                                                                                                           gCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCa                                       >  8:206817/1‑148 (MQ=255)
                                                                                                                tGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACAcc                                 <  8:242172/149‑1 (MQ=255)
                                                                                                                   cATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGa                              <  6:325588/149‑1 (MQ=255)
                                                                                                                   cATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGa                              >  7:277711/1‑149 (MQ=255)
                                                                                                                   cATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGa                              <  2:56767/149‑1 (MQ=255)
                                                                                                                     tCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAtt                            <  2:324351/149‑1 (MQ=255)
                                                                                                                      cGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAcc                                                                                 >  4:283495/1‑95 (MQ=255)
                                                                                                                      cGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAcc                                                                                 <  3:283495/95‑1 (MQ=255)
                                                                                                                      cGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAttt                           >  5:161717/1‑149 (MQ=255)
                                                                                                                            gTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgcg                     >  6:363298/1‑149 (MQ=255)
                                                                                                                             tACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCgg                    >  5:20663/1‑149 (MQ=255)
                                                                                                                                   cTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCg              >  8:298268/1‑149 (MQ=255)
                                                                                                                                      gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc             >  2:158608/1‑147 (MQ=255)
                                                                                                                                      gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc             >  2:303450/1‑147 (MQ=255)
                                                                                                                                       tggtttCCAGCTGGAAGCTGGAACCCACTGACCAGCTCAACGCCGAAGTGTCGACCGTTCTGAACATCACCGAAGACCACATCGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAGATTCCGCGGTGCCCGCCAg          >  5:195562/4‑149 (MQ=255)
                                                                                                                                         gtttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCTCCAggt        >  6:30913/2‑149 (MQ=255)
                                                                                                                                         cTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggt        >  1:133880/1‑149 (MQ=255)
                                                                                                                                         cTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggt        >  2:130552/1‑149 (MQ=255)
                                                                                                                                          tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGTCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg       >  4:101598/1‑149 (MQ=255)
                                                                                                                                          tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCTCCAggtg       <  5:30913/149‑1 (MQ=255)
                                                                                                                                           ttccagcttccagcTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACCCCGGATTTTCCGCGGTGCCCGCCAggtgg      >  8:341097/1‑149 (MQ=255)
                                                                                                                                            tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                             >  2:366467/1‑93 (MQ=255)
                                                                                                                                            tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                             <  1:366467/93‑1 (MQ=255)
                                                                                                                                            tccagctggaagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggt     >  1:164262/1‑149 (MQ=255)
                                                                                                                                              cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg    <  7:206817/148‑1 (MQ=255)
                                                                                                                                               agcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtggt  >  6:237197/1‑149 (MQ=255)
                                                                                                                                                gcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgc                      >  6:350739/1‑128 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGC‑CGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGT  >  NC_002947/1521100‑1521390

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: ‑