Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,521,290 C→A 5.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,2900CA 5.7% 529.9 / 15.4 158Q169K (CAG→AAG) murDUDP‑N‑acetylmuramoylalanine‑‑D‑glutamate ligase
Reads supporting (aligned to +/- strand):  ref base C (79/70);  new base A (3/6);  total (82/76)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 201 total aligned reads displayed.

CCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACCG  >  NC_002947/1521141‑1521438
                                                                                                                                                     |                                                                                                                                                    
ccAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttcc                                                                                                                                                      >  1:4467825/1‑150 (MQ=255)
ccAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttcc                                                                                                                                                      >  1:5384300/1‑150 (MQ=255)
ccAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttcc                                                                                                                                                      >  1:824064/1‑150 (MQ=255)
   aCGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTttccagcttccagc                                                                                                                                                   <  2:2603220/150‑1 (MQ=255)
       cAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGa                                                                                                                                               <  2:1909455/150‑1 (MQ=255)
         agagCACCGTCACGACTCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                             <  1:2549616/150‑1 (MQ=255)
         agagCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAATCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                             <  1:3933307/150‑1 (MQ=255)
         agagCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTGGAAGCTGGaaa                                                                                                                                             <  2:2065328/150‑1 (MQ=255)
            gCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCa                                                                                                                                          >  2:4861137/1‑150 (MQ=255)
             cACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCAc                                                                                                                                         >  2:3845731/1‑150 (MQ=255)
               ccGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTg                                                                                                                                       >  1:1632137/1‑150 (MQ=255)
                cGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGa                                                                                                                                      <  1:40787/150‑1 (MQ=255)
                 gTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAc                                                                                                                                     <  2:4374696/150‑1 (MQ=255)
                 gTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAc                                                                                                                                     <  1:4683595/150‑1 (MQ=255)
                  tCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGAcc                                                                                                                                    <  1:3123528/150‑1 (MQ=255)
                      gacgCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCt                                                                                                                                >  2:2116839/1‑150 (MQ=255)
                            ggTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGc                                                                                                                          <  1:4754615/150‑1 (MQ=255)
                                   gAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTg                                                                                                                   <  2:3722908/150‑1 (MQ=255)
                                    aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                                  >  1:3531350/1‑150 (MQ=255)
                                    aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                                  >  1:98984/1‑150 (MQ=255)
                                    aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                                  >  2:3924275/1‑150 (MQ=255)
                                    aaaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTgg                                                                                                                  <  2:123069/150‑1 (MQ=255)
                                     aaTGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGc                                                                                                                 <  1:2143658/150‑1 (MQ=255)
                                       tGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGa                                                                                                               >  2:4334947/1‑150 (MQ=255)
                                       tGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGa                                                                                                               >  1:4163594/1‑150 (MQ=255)
                                          ccgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCg                                                                                                            >  1:495692/1‑150 (MQ=255)
                                           cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTc                                                                                                                               >  2:1701319/1‑130 (MQ=255)
                                           cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTc                                                                                                                               <  1:1701319/130‑1 (MQ=255)
                                           cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGt                                                                                                           >  1:730780/1‑150 (MQ=255)
                                           cgccAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTGTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGt                                                                                                           >  2:2449645/1‑150 (MQ=255)
                                              cAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACGTCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCt                                                                                                        <  2:1789312/150‑1 (MQ=255)
                                              cAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTTGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCt                                                                                                        <  2:3588071/150‑1 (MQ=255)
                                               aaGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTg                                                                                                       >  2:3447882/1‑150 (MQ=255)
                                                ttggcGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAAc                                                                                                                                            >  1:466090/5‑112 (MQ=255)
                                                     ggCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATc                                                                                                 <  2:3228259/150‑1 (MQ=255)
                                                       cAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAg                                                                                               <  1:4994801/150‑1 (MQ=255)
                                                        aaaCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGc                                                                                              >  1:2222464/1‑150 (MQ=255)
                                                        aaaCGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGc                                                                                              >  1:1426466/1‑150 (MQ=255)
                                                           cgcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGaa                                                                                           <  1:3727681/150‑1 (MQ=255)
                                                            gcgTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAg                                                                                          >  1:1553179/1‑150 (MQ=255)
                                                              gTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAc                                                                                        <  1:458677/150‑1 (MQ=255)
                                                                  ccGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACcacc                                                                                    >  1:1478641/1‑149 (MQ=255)
                                                                      gggCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCt                                                                                                        <  1:2243766/126‑1 (MQ=255)
                                                                         cGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTGCAGCTTCCAGCTGGAAACCACTGAACAGCTCTACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCg                                                                             <  2:3464672/150‑1 (MQ=255)
                                                                         cGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCg                                                                             <  1:3335651/150‑1 (MQ=255)
                                                                              aCCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCtata                                                                        <  2:2224024/150‑1 (MQ=255)
                                                                              aCCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCtata                                                                        <  2:2946676/150‑1 (MQ=255)
                                                                                 tCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTCTAGCg                                                                     >  2:2746592/1‑150 (MQ=255)
                                                                                         ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCg                                                             >  1:1636434/1‑150 (MQ=255)
                                                                                         ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCg                                                             >  2:2517674/1‑150 (MQ=255)
                                                                                         ccGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCg                                                             >  1:1299007/1‑150 (MQ=255)
                                                                                          cGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCgg                                                            >  1:2723012/1‑150 (MQ=255)
                                                                                              aCTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTa                                                        >  1:4925043/1‑150 (MQ=255)
                                                                                                      tgctgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTgg                                                <  2:1299007/150‑1 (MQ=255)
                                                                                                         tgGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCa                                             >  1:4330664/1‑150 (MQ=255)
                                                                                                             cGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAaca                                         >  2:5224316/1‑150 (MQ=255)
                                                                                                             cGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAaca                                         >  2:1047555/1‑150 (MQ=255)
                                                                                                                tGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCg                                      <  1:3447882/150‑1 (MQ=255)
                                                                                                                  aCATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGa                                    <  2:658042/150‑1 (MQ=255)
                                                                                                                   cATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGAt                                   >  1:446881/1‑150 (MQ=255)
                                                                                                                     tCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGtt                                                                                                          <  2:3978181/77‑1 (MQ=255)
                                                                                                                         gCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcg                             >  1:318588/1‑150 (MQ=255)
                                                                                                                          cTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgc                            <  2:3354179/150‑1 (MQ=255)
                                                                                                                          cTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgc                            <  2:1518886/150‑1 (MQ=255)
                                                                                                                          cTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgc                            <  2:4985692/150‑1 (MQ=255)
                                                                                                                           tGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCcgcg                           <  1:535323/150‑1 (MQ=255)
                                                                                                                            gTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGc                                                                   <  1:1003982/109‑1 (MQ=255)
                                                                                                                             tACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGt                         >  1:3973016/1‑150 (MQ=255)
                                                                                                                               cGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGc                       <  1:546723/150‑1 (MQ=255)
                                                                                                                                   cTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc                   <  2:431242/150‑1 (MQ=255)
                                                                                                                                   cTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGc                   >  1:1996076/1‑150 (MQ=255)
                                                                                                                                    tGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGcc                  >  2:4658570/1‑150 (MQ=255)
                                                                                                                                     ggAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCa                 <  1:2449645/150‑1 (MQ=255)
                                                                                                                                      gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAg                >  1:3234754/1‑150 (MQ=255)
                                                                                                                                      gAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAg                >  1:2192745/1‑150 (MQ=255)
                                                                                                                                       tggtttCCAGCTGGAAGCTGGAAAGCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCaa                                            <  2:1129872/118‑1 (MQ=255)
                                                                                                                                       tggtttCCAGCTGGAAGCTGGAAAGCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCaa                                            >  1:1129872/4‑121 (MQ=255)
                                                                                                                                         gtttCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtg             >  2:1016152/2‑150 (MQ=255)
                                                                                                                                          tttCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGTGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg            <  1:693944/150‑1 (MQ=255)
                                                                                                                                          tttCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg            >  2:1607789/1‑150 (MQ=255)
                                                                                                                                          tttCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg            >  1:1824698/1‑150 (MQ=255)
                                                                                                                                          tttCCAGCTGGAAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtgg            >  2:5276671/1‑150 (MQ=255)
                                                                                                                                           ttgcagcttccagctGGAAACCACTGACCAGCTCATCGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGAGCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCATACACCGGATTTTCCGCGGTGCGCGCGAggtggt           >  2:3278584/4‑150 (MQ=255)
                                                                                                                                           ttccagctggaagcTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggt           >  1:609926/1‑150 (MQ=255)
                                                                                                                                            tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCCGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg          <  2:5384300/150‑1 (MQ=255)
                                                                                                                                            tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg          <  1:5213382/150‑1 (MQ=255)
                                                                                                                                            tccagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg          <  1:4945062/150‑1 (MQ=255)
                                                                                                                                              cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACAcc                                       <  2:2286765/119‑1 (MQ=255)
                                                                                                                                              cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtggt        >  2:2699341/1‑150 (MQ=255)
                                                                                                                                              cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtggt        >  1:2953455/1‑150 (MQ=255)
                                                                                                                                              cagcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtggt        >  1:1902165/1‑150 (MQ=255)
                                                                                                                                               agcttccagctGGAAACCTCTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACAGCAGCGAAGACCTCATGGTCCTCTATAGCGGCCTTCCGGCCTAGCACCTGGCGAAACCCCGGATTTTCCGCGGTTCCCGACAGGTGGTGGTc       <  2:1093874/150‑1 (MQ=255)
                                                                                                                                               agcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCCGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTc       <  2:2050529/150‑1 (MQ=255)
                                                                                                                                                gcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAggtggtg          >  1:12040/1‑146 (MQ=255)
                                                                                                                                                gcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCa      <  1:4948343/150‑1 (MQ=255)
                                                                                                                                                gcttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCa      <  2:4467825/150‑1 (MQ=255)
                                                                                                                                                 cttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCaa     >  1:4223333/1‑150 (MQ=255)
                                                                                                                                                 cttccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCaa     <  2:824064/150‑1 (MQ=255)
                                                                                                                                                   tccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAAcc   <  2:4028050/150‑1 (MQ=255)
                                                                                                                                                    ccagctGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACCg  >  2:1376349/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
CCAACGCCAAGAGCACCGTCACGACGCTGGTCGGTGAAATGGCCGCCAAGGCCGGCAAACGCGTGGCCGTGGGCGGTAACCTCGGTACCCCGGCACTCGACCTGCTGGCCGATGACATCGAGCTGTACGTGCTGGAGCTTTCCAGCTTCCAGCTGGAAACCACTGACCAGCTCAACGCCGAAGTGGCGACCGTTCTGAACATCAGCGAAGACCACATGGACCGCTATAGCGGCCTTCCGGCCTACCACCTGGCCAAACACCGGATTTTCCGCGGTGCCCGCCAGGTGGTGGTCAACCG  >  NC_002947/1521141‑1521438

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