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

ACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGG  >  NC_002947/1521102‑1521389
                                                                                                                                                   |                                                                                                                                             
aCGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttccag                                                                                                                                             >  5:90859/1‑149 (MQ=255)
aCGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttccag                                                                                                                                             <  6:90859/149‑1 (MQ=255)
  gCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTtccagcttccagct                                                                                                                                           <  5:198849/149‑1 (MQ=255)
       gagCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                      <  5:145449/149‑1 (MQ=255)
        agCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAAc                                                                                                                                     <  7:250753/149‑1 (MQ=255)
        agCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaag                                                                                                                                     <  2:237023/149‑2 (MQ=255)
            ccGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGCCATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                 <  1:95520/149‑1 (MQ=255)
            ccGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACt                                                                                                                                 >  6:402686/1‑149 (MQ=255)
                cACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAcc                                                                                                                             <  8:352662/149‑1 (MQ=255)
                 acgacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCCTGTACGTGCTGGAGCTCTCCAGCTTCCAGCTGGAAACCACTGAcc                                                                                                                             <  1:291564/149‑1 (MQ=255)
                   gacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGc                                                                                                                          >  2:266025/1‑149 (MQ=255)
                    acgCTGGTCGGTGAAATGACCGCCAAGGCCGGCAAACGCGTGGCCCTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCt                                                                                                                         >  1:292939/1‑149 (MQ=255)
                          gTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAAc                                                                                                                     <  7:386566/147‑1 (MQ=255)
                          gTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACg                                                                                                                    <  6:414794/148‑1 (MQ=255)
                                 aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAcc                                                                                                                             >  8:198938/1‑132 (MQ=255)
                                 aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAcc                                                                                                                             <  7:198938/132‑1 (MQ=255)
                                  aaTGGCCTCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGTCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTg                                                                                                            <  8:262877/148‑1 (MQ=255)
                                   aTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGc                                                                                                          <  1:266025/149‑1 (MQ=255)
                                    tGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCg                                                                                                         <  5:150879/149‑1 (MQ=255)
                                         gccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGt                                                                                                    >  7:88643/1‑149 (MQ=255)
                                                  ggCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                             <  3:21707/147‑1 (MQ=255)
                                                  ggCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                             >  4:21707/1‑147 (MQ=255)
                                                   gCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTg                                                                                                >  6:230126/1‑143 (MQ=255)
                                                   gCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTg                                                                                                <  5:230126/143‑1 (MQ=255)
                                                    cAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCa                                                                                         <  5:402686/149‑1 (MQ=255)
                                                      aaCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGc                                                                                       <  5:151873/149‑1 (MQ=255)
                                                         gcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttccag                                                                                                                                             >  5:141373/1‑92 (MQ=255)
                                                         gcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttccag                                                                                                                                             <  6:141373/92‑1 (MQ=255)
                                                             ggCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGa                                                                                                        <  3:5175/125‑1 (MQ=255)
                                                             ggCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGa                                                                                                        >  4:5175/1‑125 (MQ=255)
                                                              gCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACca                                                                               <  3:174652/149‑1 (MQ=255)
                                                               ccGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                             <  3:215933/134‑1 (MQ=255)
                                                               ccGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAAc                                                                                             >  4:215933/1‑134 (MQ=255)
                                                                    ggCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGa                                                                         >  3:173151/1‑149 (MQ=255)
                                                                           aCCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCta                                                                   <  2:292939/148‑1 (MQ=255)
                                                                            ccTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCtata                                                                 <  4:189313/149‑1 (MQ=255)
                                                                             cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAg                                                                >  4:56562/1‑149 (MQ=255)
                                                                             cTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAg                                                                >  2:390441/1‑149 (MQ=255)
                                                                              tCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGc                                                               >  6:70777/1‑149 (MQ=255)
                                                                                      ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAc                                                                                 >  1:194734/1‑123 (MQ=255)
                                                                                      ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAc                                                                                 <  2:194734/123‑1 (MQ=255)
                                                                                      ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTcc                                                       >  5:423950/1‑149 (MQ=255)
                                                                                        ggCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATATCGGCCTTCCgg                                                     <  6:91561/149‑1 (MQ=255)
                                                                                         gCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGc                                                    <  4:173151/149‑1 (MQ=255)
                                                                                            cTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTa                                                 <  3:56562/149‑1 (MQ=255)
                                                                                             tCGACCTGCAGGCCGATGACATCGAGCTGAACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGACGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCTGCCTTCCGGCCTa                                                 >  5:402659/1‑148 (MQ=255)
                                                                                               gACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAATATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTacca                                              <  1:298539/149‑1 (MQ=255)
                                                                                               gACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTacca                                              >  3:222742/1‑149 (MQ=255)
                                                                                               gACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTac                                                >  1:18471/1‑147 (MQ=255)
                                                                                               gACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTac                                                <  2:18471/147‑1 (MQ=255)
                                                                                                aCCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTaccac                                             <  6:423950/149‑1 (MQ=255)
                                                                                                   tgttgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGa                                                                         <  3:374917/115‑1 (MQ=255)
                                                                                                   tgttgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGa                                                                         >  4:374917/4‑118 (MQ=255)
                                                                                                    gctgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTgg                                         >  4:98773/1‑149 (MQ=255)
                                                                                                     ctgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTg                                          <  4:222742/147‑1 (MQ=255)
                                                                                                           gATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAaca                                  >  7:164563/1‑149 (MQ=255)
                                                                                                           gATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAaca                                  >  5:173494/1‑149 (MQ=255)
                                                                                                             tGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACAcc                                <  8:88643/149‑1 (MQ=255)
                                                                                                                 aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGTAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAt                            >  7:274019/1‑149 (MQ=255)
                                                                                                                 aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGtt                                                                                                   <  4:156065/78‑1 (MQ=255)
                                                                                                                 aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGtt                                                                                                   >  3:156065/1‑78 (MQ=255)
                                                                                                                 aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                            >  8:113025/1‑117 (MQ=255)
                                                                                                                 aTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                            <  7:113025/117‑1 (MQ=255)
                                                                                                                  tCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAtt                           >  8:55883/1‑149 (MQ=255)
                                                                                                                      gCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTcc                       >  7:189400/1‑149 (MQ=255)
                                                                                                                          gTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgcg                    >  6:6012/1‑149 (MQ=255)
                                                                                                                              gTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGc                >  4:233775/1‑149 (MQ=255)
                                                                                                                                    gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGTGGCCTTCCGGCCTACCACCTgg                                         >  3:306029/1‑118 (MQ=255)
                                                                                                                                    gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGTGGCCTTCCGGCCTACCACCTgg                                         <  4:306029/118‑1 (MQ=255)
                                                                                                                                     tggtttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAg         >  4:256631/4‑149 (MQ=255)
                                                                                                                                     aGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAg         <  5:70777/149‑1 (MQ=255)
                                                                                                                                       gtttCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggt       >  5:353139/2‑149 (MQ=255)
                                                                                                                                        tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg      >  2:9773/1‑149 (MQ=255)
                                                                                                                                         ttccagcttccagcTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGATCATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg     >  6:200846/1‑149 (MQ=255)
                                                                                                                                          tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCTTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCgg                                                             <  5:394616/92‑1 (MQ=255)
                                                                                                                                          tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggt    >  7:70734/1‑149 (MQ=255)
                                                                                                                                          tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCATTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCgg                                                             >  6:394616/1‑92 (MQ=255)
                                                                                                                                           ccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg   >  8:91661/1‑149 (MQ=255)
                                                                                                                                            cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtgg  <  6:353139/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                             
ACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGG  >  NC_002947/1521102‑1521389

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