Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,248 C→A 8.7% Q169K (CAG→AAG)  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,2480CA 8.7% 202.0 / 8.9 69Q169K (CAG→AAG) murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base C (29/33);  new base A (3/3);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.83e-01

GCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACC  >  NC_002947/1521104‑1521395
                                                                                                                                                |                                                                                                                                                      
gCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTtccagcttccagct                                                                                                                                                    <  2:77271/149‑1 (MQ=255)
  cAAGAGCTCCGGCTCGACGCTGGTCGGTGAAATGGCCGCGATGGCCGGCAATCGCGTGGCGGTGGGCGGTAACCTGGGTTCCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGGTTCCAGCTgg                                                                                                                                                  <  2:384103/149‑1 (MQ=255)
  cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTgg                                                                                                                                                  <  2:2499/149‑1 (MQ=255)
  cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTg                                                                                                                                                   >  8:150198/1‑148 (MQ=255)
  cAAGAGCACCGTCACGACGCTGGTCGGTGAAAAGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCATCCAGCTg                                                                                                                                                   >  1:384103/1‑148 (MQ=255)
      agCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAAc                                                                                                                                              <  4:102840/149‑1 (MQ=255)
           cGTCACGACTCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaagctcca                                                                                                                                         <  2:51899/149‑7 (MQ=255)
           cGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                          >  2:130107/1‑148 (MQ=255)
           cGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                          >  2:33562/1‑148 (MQ=255)
             tCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAc                                                                                                                                       <  7:117161/149‑1 (MQ=255)
                 gacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGc                                                                                                                                   <  5:259246/149‑1 (MQ=255)
                          aatgcccaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTATCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCg                                                                                                                          <  1:176016/142‑1 (MQ=255)
                            gTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGa‑‑a                                                                                                                      <  7:150198/149‑1 (MQ=255)
                               aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTCCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTg                                                                                                                  <  2:361291/149‑1 (MQ=255)
                                aaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTgg                                                                                                                 >  1:212730/1‑149 (MQ=255)
                                       gccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGt                                                                                                          >  4:34766/1‑149 (MQ=255)
                                         cAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTc                                                                                                        >  3:1281/1‑149 (MQ=255)
                                          aaGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCt                                                                                                       <  4:1281/149‑1 (MQ=255)
                                             gCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGaa                                                                                                    <  8:287067/149‑1 (MQ=255)
                                                 gCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                          >  8:305798/1‑110 (MQ=255)
                                                 gCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                          <  7:305798/110‑1 (MQ=255)
                                                    aaCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGc                                                                                             <  5:107100/149‑1 (MQ=255)
                                                       gcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCg                                                                                            >  5:102585/1‑147 (MQ=255)
                                                         gTGGCCGTGGGCGGTAACCTCGGAACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGa                                                                                        <  7:325137/149‑1 (MQ=255)
                                                          tGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGAc                                                                                       <  7:7308/149‑1 (MQ=255)
                                                            gCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACca                                                                                     <  3:34766/149‑1 (MQ=255)
                                                              cGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACcaca                                                                                   >  7:180557/1‑149 (MQ=255)
                                                                tGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATg                                                                                 >  5:388964/1‑149 (MQ=255)
                                                                  ggCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATgg                                                                                <  8:152692/148‑1 (MQ=255)
                                                                   gCGTTAACCTCGGGACCCCTGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGAc                                                                              <  4:236752/149‑1 (MQ=255)
                                                                   gCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGTTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCTTTCTGAACATCAGCGAAGACCACATGGAc                                                                              <  8:364519/149‑1 (MQ=255)
                                                                       tAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGc                                                                           >  1:188849/1‑148 (MQ=255)
                                                                         aCCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCtat                                                                        <  7:260116/149‑1 (MQ=255)
                                                                           cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAg                                                                      >  3:4346/1‑149 (MQ=255)
                                                                                  ccccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCtt                                                               <  5:394404/149‑1 (MQ=255)
                                                                                   cccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTc                                                              >  2:51388/1‑149 (MQ=255)
                                                                                   cccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTc                                                              >  7:137584/1‑149 (MQ=255)
                                                                                         aCTAGACCTGATGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑TGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCt                                                        <  4:118891/149‑1 (MQ=255)
                                                                                          cTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCt                                                                                                       >  4:45319/1‑101 (MQ=255)
                                                                                          cTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCt                                                                                                       <  3:45319/101‑1 (MQ=255)
                                                                                          cTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTa                                                       <  6:102585/149‑1 (MQ=255)
                                                                                             gACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTacca                                                    >  1:257791/1‑149 (MQ=255)
                                                                                             gACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTac                                                      <  1:33562/147‑1 (MQ=255)
                                                                                              aCCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTaccac                                                   <  7:135664/149‑1 (MQ=255)
                                                                                              aCCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTaccac                                                   <  1:153153/149‑1 (MQ=255)
                                                                                                       ccGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCa                                                                                               >  4:240860/1‑96 (MQ=255)
                                                                                                       ccGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCa                                                                                               <  3:240860/96‑1 (MQ=255)
                                                                                                        cGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                                  >  1:272369/1‑124 (MQ=255)
                                                                                                        cGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                                  <  2:272369/124‑1 (MQ=255)
                                                                                                        cGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCaaa                                          >  3:251470/1‑148 (MQ=255)
                                                                                                         gATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAaca                                        >  2:276694/1‑149 (MQ=255)
                                                                                                              cATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCgg                                    <  2:188849/148‑1 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTg                                                <  5:179007/135‑1 (MQ=255)
                                                                                                               aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTg                                                >  6:179007/1‑135 (MQ=255)
                                                                                                                        ttcGTGCTGGTGCTTTTGAGCTTCGAGCTGGTTAAAACTGATCAGCTCAACTCCGATTAGGTGCCCACCGTACTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcg                            <  1:130107/147‑1 (MQ=255)
                                                                                                                        tACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCgg                         >  1:270666/1‑149 (MQ=255)
                                                                                                                        tACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCgg                         >  8:206438/1‑149 (MQ=255)
                                                                                                                            tGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGcc                     >  2:282882/1‑149 (MQ=255)
                                                                                                                               cagtggtttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCgg                         >  5:103605/7‑142 (MQ=255)
                                                                                                                               cagtggtttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCgg                         <  6:103605/136‑1 (MQ=255)
                                                                                                                                ggAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc                  >  5:129194/1‑148 (MQ=255)
                                                                                                                                 gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc                  <  8:180557/147‑1 (MQ=255)
                                                                                                                                 gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc                  >  4:273171/1‑147 (MQ=255)
                                                                                                                                 gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCa                <  2:270666/149‑1 (MQ=255)
                                                                                                                                  aGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGTGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                                  >  1:9075/1‑98 (MQ=255)
                                                                                                                                  aGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGTGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                                  <  2:9075/98‑1 (MQ=255)
                                                                                                                                   gCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAgg              >  8:99114/1‑149 (MQ=255)
                                                                                                                                    gtttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCTGTGCCCGCAAggg             >  6:283401/2‑148 (MQ=255)
                                                                                                                                     tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAac                                         <  4:204793/120‑1 (MQ=255)
                                                                                                                                     tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAac                                         >  3:204793/1‑120 (MQ=255)
                                                                                                                                      ttccagctggaagcTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg           >  6:136243/1‑149 (MQ=255)
                                                                                                                                       tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggt          >  6:284877/1‑149 (MQ=255)
                                                                                                                                       tccagctggaagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg           <  5:136243/148‑1 (MQ=255)
                                                                                                                                         cagcttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtgg        <  7:99114/149‑1 (MQ=255)
                                                                                                                                         cagcttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg         <  4:4346/148‑1 (MQ=255)
                                                                                                                                          agcttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtggt       <  7:206438/149‑1 (MQ=255)
                                                                                                                                            cttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTc      >  6:54197/1‑148 (MQ=255)
                                                                                                                                             ttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAttt                                >  3:338799/1‑121 (MQ=255)
                                                                                                                                             ttccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAttt                                <  4:338799/121‑1 (MQ=255)
                                                                                                                                               ccagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAAcc  >  4:125856/1‑149 (MQ=255)
                                                                                                                                                cagctGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAAcc  <  2:257791/148‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                      
GCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGA‑‑AGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACC  >  NC_002947/1521104‑1521395

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: