Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,081,173 C→G 7.3% G255G (GGG→GGC dpkA ← delta 1‑piperideine‑2‑carboxylate reductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,081,1730CG 7.3% 200.3 / 4.4 65G255G (GGG→GGCdpkAdelta 1‑piperideine‑2‑carboxylate reductase
Reads supporting (aligned to +/- strand):  ref base C (27/33);  new base G (2/3);  total (29/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCGGTAGTCGGC  >  NC_002947/4081027‑4081320
                                                                                                                                                  |                                                                                                                                                   
cgcTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTTGGCAAACGGCTGGGCTTCGGCCTTGCCCGGGTCGATGACGATGATGAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGccc                                                                                                                                                   <  3:116018/149‑1 (MQ=255)
 gcTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCg                                                                                                                                                  >  3:156155/1‑149 (MQ=255)
     gCCCTGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACGATCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGAc                                                                                                                                               <  6:233304/148‑1 (MQ=255)
     gCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGTCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCTGAc                                                                                                                                               <  7:229562/148‑1 (MQ=255)
     gCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGa                                                                                                                                                <  8:273706/147‑1 (MQ=255)
     gCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGAc                                                                                                                                               <  6:137201/148‑1 (MQ=255)
         ggattGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCAGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATc                                                                                                                                          >  7:274591/5‑149 (MQ=255)
         ggaaTGCGCGTCAGCCACACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCAGGTCCATGGCGTTTTAGCCCCCGGATGCCCGGACCAATc                                                                                                                                          >  7:274592/4‑149 (MQ=255)
                gcgTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGGCCGGACCAATCGAACTcc                                                                                                                                   >  7:265340/1‑149 (MQ=255)
                gcgTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTGCGGCCTTGCCCGGGTCGATGACGATGATCCACTGCCCGGTCCATGGCGTTTTCGCACCCGGATGGCCGGACCAAGCGACCTca                                                                                                                                   >  7:265342/1‑148 (MQ=255)
                     aGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCCGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTAGCCCCCGGATGCCCGGACCAATCGAACGCCCAg                                                                                                                                >  4:147444/1‑147 (MQ=255)
                      gCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAg                                                                                                                                >  5:161275/1‑146 (MQ=255)
                      gCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGgag                                                                                                                             <  2:7278/149‑1 (MQ=255)
                          cACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAAt                                                                                                                         <  1:374/149‑1 (MQ=255)
                          cACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAAt                                                                                                                         <  2:326609/149‑1 (MQ=255)
                          cACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAAt                                                                                                                         >  1:208931/1‑149 (MQ=255)
                            ccGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGGCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGa                                                                                                                       <  8:301243/149‑1 (MQ=255)
                               ccTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGAccg                                                                                                                    >  4:18288/1‑149 (MQ=255)
                                  gCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCCTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGAccgccg                                                                                                                 >  1:355292/1‑149 (MQ=255)
                                   cATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGAccgcc                                                                                                                  <  6:368864/147‑1 (MQ=255)
                                    aTGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGt                                                                                                               <  3:18288/149‑1 (MQ=255)
                                     tgtgCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTc                                                                                                              <  1:276792/149‑1 (MQ=255)
                                        gCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAAt                                                                                                                                           <  3:85724/117‑1 (MQ=255)
                                        gCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAAt                                                                                                                                           >  4:85724/1‑117 (MQ=255)
                                         cTCCACCAGCTCGCGGCTGCACTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAgcg                                                                                                          >  5:437573/1‑149 (MQ=255)
                                            caccaGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGGCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAgcgc                                                                                                         <  5:86831/147‑1 (MQ=255)
                                            caccaGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCCGGAGCAATGACCGCCGGTCAGcgccg                                                                                                       >  7:390418/1‑149 (MQ=255)
                                            caccaGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGcgccg                                                                                                       >  2:323050/1‑149 (MQ=255)
                                              ccaGCTGGCGGTGGCGCTTGGGAAACCGGTGGGCTTCGGACTTGGCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCTATCGAACTCCCAGGAGAAATGACCGCCAGTCAGcgccgcc                                                                                                     <  4:327452/149‑1 (MQ=255)
                                              ccaGCTCGCGGCTGCACTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGcgccgcc                                                                                                     <  6:437573/149‑1 (MQ=255)
                                               caGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAgg                                                                                                                               >  6:54407/1‑122 (MQ=255)
                                               caGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAgg                                                                                                                               <  5:54407/122‑1 (MQ=255)
                                                aGCTCGCGGCTTCGCTGGGCAAACCTCTGGCCTTCGGCCTTGCACGGGTCGATGACTATTATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCGAGGAGAAATGACCGCCGGTCAGcgcctcggc                                                                                                   <  8:274592/149‑6 (MQ=255)
                                                aGCTCGCGGCTGCTCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTAGATGACGATGATCATCTGCACGGTCCATGGCGTTTTCGCCCCCGGATGGCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGcgccgccgc                                                                                                   <  8:274591/149‑1 (MQ=255)
                                                 gCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCgccgccgcc                                                                                                  <  2:157090/149‑1 (MQ=255)
                                                   tCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCcag                                                                                                >  1:420141/1‑149 (MQ=255)
                                                     gcgGCTGCGCTTGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGAGGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGGCCGGACCAATCGACCTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCcagca                                                                                              <  8:265342/149‑1 (MQ=255)
                                                     gcgGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGGCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCcagca                                                                                              <  8:265340/149‑1 (MQ=255)
                                                     gcgGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCcag                                                                                                <  2:234009/147‑1 (MQ=255)
                                                     gcgGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTAGCCCGGGTCGATGACGATGATCAACTGCCCGGTACATGGCGATTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCCGTCAGCGCAGCCGCcagca                                                                                              >  3:288101/1‑149 (MQ=255)
                                                      gggCTGAGCTGGGCAAACCGCTGGCCTTCGGCATTGCCCGGGTCGATGACGATGATCAACTGCCCGGTACATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCcagcag                                                                                             <  8:168989/148‑1 (MQ=255)
                                                       ggCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCagcagc                                                                                            <  5:130549/149‑1 (MQ=255)
                                                        gCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACGAATCGAACTACGAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCt                                                                                           >  8:180054/1‑149 (MQ=255)
                                                                  gCAAACCGCTGGCCTTCGGCCTTGCCGGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACcatcat                                                                                 <  5:395538/149‑1 (MQ=255)
                                                                   cAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCatcatc                                                                                <  6:270816/149‑1 (MQ=255)
                                                                                    gCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCtt                                                               <  2:420141/149‑1 (MQ=255)
                                                                                           ccGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCg                                                        >  5:108223/1‑149 (MQ=255)
                                                                                                 cgatgacgatgaTCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAAt                                                  <  2:208931/149‑1 (MQ=255)
                                                                                                   atgacgatgaTCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGCCCATCATAGACAACGCCGAGCCCTTGTGCCCGCCAAAtcg                                                >  6:264802/1‑147 (MQ=255)
                                                                                                     gacgatgaTCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCAAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACcatcat                                                                                 <  4:31501/114‑1 (MQ=255)
                                                                                                     gacgatgaTCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCAAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACcatcat                                                                                 >  3:31501/1‑114 (MQ=255)
                                                                                                       cgatgaTCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGgcag                                             >  5:200206/1‑148 (MQ=255)
                                                                                                       cgatgaTCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGgcag                                             >  6:24840/1‑148 (MQ=255)
                                                                                                             tCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGc                                       <  5:264802/148‑1 (MQ=255)
                                                                                                               aaCTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGcccc                                    >  3:193725/1‑149 (MQ=255)
                                                                                                                aCTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCg                                   >  1:356533/1‑149 (MQ=255)
                                                                                                                      cGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTcc                             <  6:161275/149‑1 (MQ=255)
                                                                                                                             tGGCGTTTTCGCCCCCGGATGCCCGGACCAATCAAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCg                      <  4:354879/149‑1 (MQ=255)
                                                                                                                                      cGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGt             >  1:376317/1‑149 (MQ=255)
                                                                                                                                          cccGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCGGt         <  1:323050/149‑1 (MQ=255)
                                                                                                                                           ccGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCgg          >  8:300351/1‑147 (MQ=255)
                                                                                                                                           ccGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCGGTa        <  3:147444/149‑1 (MQ=255)
                                                                                                                                                tGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCGGTAGTCgg   <  5:24840/149‑1 (MQ=255)
                                                                                                                                                 gCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCGGTAGTCGGc  >  2:23647/1‑149 (MQ=255)
                                                                                                                                                 gCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAACAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCg           >  1:333015/1‑140 (MQ=255)
                                                                                                                                                 gCCCGGAACAATCGAACTCGCATGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCg           <  2:333015/140‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCTCGCCCGGCATGCGCGTCAGCCCCACCGCCTGCATGTGCTCCACCAGCTCGCGGCTGCGCTGGGCAAACCGCTGGCCTTCGGCCTTGCCCGGGTCGATGACGATGATCAACTGCCCGGTCCATGGCGTTTTCGCCCCCGGATGCCCGGACCAATCGAACTCCCAGGAGAAATGACCGCCGGTCAGCGCCGCCGCCAGCAGCTCGACCATCATCGACAACGCCGAGCCCTTGTGCCCGCCAAATGGCAGCAAGGCCCCGCCTTCCAGGATCGCCTTCGGGTCGGTAGTCGGC  >  NC_002947/4081027‑4081320

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