Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 418,296 C→G 6.5% Q27E (CAG→GAG)  PP_0345 → heptose kinase WapP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947418,2960CG 6.5% 249.0 / 5.8 77Q27E (CAG→GAG) PP_0345heptose kinase WapP
Reads supporting (aligned to +/- strand):  ref base C (34/38);  new base G (2/3);  total (36/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.72e-01

GCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAA  >  NC_002947/418149‑418442
                                                                                                                                                   |                                                                                                                                                  
gctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCa                                                                                                                                                   <  7:66363/149‑1 (MQ=255)
 ctCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTcc                                                                                                                                                    <  6:5078/147‑1 (MQ=255)
 ctCGACTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCAATCTGATGGCGAGTTCGACACTGGCGCATGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCTCCCTCCAg                                                                                                                                                  >  7:272317/1‑149 (MQ=255)
    gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGc                                                                                                                                               >  5:111891/1‑149 (MQ=255)
    gCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGACGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAgg                                                                                                                                                 <  6:354056/147‑1 (MQ=255)
    gCCGGGCGGCAGCGCACGGCCCAGAAACTCTACGAGCGCAACCACCGCTTGGGGCATGCACTCTTTGGGGGTGTTGCACACTGGCGCGGGGTTACCAACACCTGGCGGCGCACTTCGGCAGCCTCGCCGCGGTGTTCGCCCTCCAGGGc                                                                                                                                               <  4:344696/149‑1 (MQ=255)
         gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTGGAGGGCGaa                                                                                                                                            <  8:62233/147‑1 (MQ=255)
         gctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAAc                                                                                                                                           <  6:362459/148‑1 (MQ=255)
          ctggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGCACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACgg                                                                                                                                         <  5:347985/149‑1 (MQ=255)
           tggAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGc                                                                                                                                        >  4:224795/1‑149 (MQ=255)
               gcgcAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgac                                                                                                                                    <  5:99462/149‑1 (MQ=255)
                  cAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccg                                                                                                                                 >  7:283836/1‑149 (MQ=255)
                   aaGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                                                <  6:111891/149‑1 (MQ=255)
                   aaGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgc                                                                                                                                >  3:347540/1‑149 (MQ=255)
                     ggCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcga                                                                                                                              >  2:267270/1‑149 (MQ=255)
                         cAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcga                                                                                                                              <  4:185659/145‑1 (MQ=255)
                         cAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTgacccgcga                                                                                                                              >  3:185659/1‑145 (MQ=255)
                           gAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCt                                                                                                                        >  7:9323/1‑149 (MQ=255)
                                    cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGCCCCGCGACCCGCTCAGTgtcgg                                                                                                               >  6:373775/1‑145 (MQ=255)
                                    cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGa                                                                                                                  <  5:58466/146‑1 (MQ=255)
                                    cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAg                                                                                                                <  7:294945/148‑1 (MQ=255)
                                    cGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGAAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAg                                                                                                                <  3:242666/148‑1 (MQ=255)
                                     gATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTg                                                                                                              <  2:109962/149‑1 (MQ=255)
                                       tCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGt                                                                                                            <  5:183742/149‑1 (MQ=255)
                                       tCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGt                                                                                                            <  6:195484/149‑1 (MQ=255)
                                         gcaagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTAc                                                                                                          <  5:148605/149‑1 (MQ=255)
                                         cccagcaaCGCTATGGGGCTGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTGCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTAc                                                                                                          <  8:327163/146‑1 (MQ=255)
                                           aagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAAc                                                                                                                                           >  8:254301/1‑114 (MQ=255)
                                           aagcaaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAAc                                                                                                                                           <  7:254301/114‑1 (MQ=255)
                                               aaCGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCg                                                                                                    >  6:380722/1‑149 (MQ=255)
                                                  gCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGc                                                                                                 <  4:159797/149‑1 (MQ=255)
                                                      tGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGa                                                                                             <  8:9323/149‑1 (MQ=255)
                                                       ggggATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGAc                                                                                            >  6:284337/1‑149 (MQ=255)
                                                            tGCACTCTGATGGCTAGTTGGACACTTGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGAt                                                                                       <  8:272317/149‑1 (MQ=255)
                                                               aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCAt                                                                                                      <  4:265818/131‑1 (MQ=255)
                                                               aCTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCAt                                                                                                      >  3:265818/1‑131 (MQ=255)
                                                                 tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                   >  3:66723/1‑148 (MQ=255)
                                                                 tctGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                   >  8:238348/1‑148 (MQ=255)
                                                                   tGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAac                                                                                   <  4:66723/146‑1 (MQ=255)
                                                                       ggCGAGTTGGACACTGTCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACtactac                                                                             <  8:337948/148‑1 (MQ=255)
                                                                                gACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                    <  4:88134/148‑1 (MQ=255)
                                                                                gACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                    >  2:96795/1‑148 (MQ=255)
                                                                                gACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                    >  5:292290/1‑148 (MQ=255)
                                                                                   acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTaca                                                                >  4:61199/1‑149 (MQ=255)
                                                                                   acTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTaca                                                                >  4:66170/1‑149 (MQ=255)
                                                                                     tGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACAcc                                                              >  7:192022/1‑149 (MQ=255)
                                                                                     tGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACAcc                                                              >  8:166451/1‑149 (MQ=255)
                                                                                      atcgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGCCCCGCGACCCGCTCAGGGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCg                                                             >  4:16820/3‑149 (MQ=255)
                                                                                       gcgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                      >  6:135662/1‑91 (MQ=255)
                                                                                       gcgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCa                                                                                                                      <  5:135662/91‑1 (MQ=255)
                                                                                        cgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCggg                                                           <  5:284337/149‑1 (MQ=255)
                                                                                        cgcCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCggg                                                           <  1:96795/149‑1 (MQ=255)
                                                                                           cGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGcc                                                        <  4:159750/149‑1 (MQ=255)
                                                                                               gTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGc                                                     >  8:103523/1‑148 (MQ=255)
                                                                                                  cGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGc                                                 <  3:66170/149‑1 (MQ=255)
                                                                                                  cGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGc                                                 <  3:61199/149‑1 (MQ=255)
                                                                                                  cGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGc                                                 <  3:224795/149‑1 (MQ=255)
                                                                                                  cGAACACCTGGCGGCGGACTTCGGCACCCTCGACGCGGTGTTCGCCCTCCAGTGCGAACGGCTCACCCGCGACCCGCTCAGTGAAGTGGTTCGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACCCCGGGGCCGGCAAg                                                  >  5:358842/1‑148 (MQ=255)
                                                                                                    aaCACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                    <  2:2377/128‑1 (MQ=255)
                                                                                                    aaCACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAgcgc                                                                    >  1:2377/1‑128 (MQ=255)
                                                                                                     acacCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGcaca                                              >  1:319049/1‑149 (MQ=255)
                                                                                                         cTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATgcg                                          >  3:292542/1‑149 (MQ=255)
                                                                                                         cTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCACTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGAAAGCACATgcg                                          >  6:175690/1‑149 (MQ=255)
                                                                                                                  gACTTCGGCAGCCTAGACGCGGTGTTAGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGTGATCAACTACTACCTCAAGCGCTACACCGGTGCCGGCAAGCACATGCGCCGGTACCt                                 <  5:22716/149‑1 (MQ=255)
                                                                                                                        ggCAGCCTCGACGCGGTGTTCGCCCTCGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggc                            <  5:380722/148‑1 (MQ=255)
                                                                                                                         gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGCCCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg                          >  2:203088/1‑149 (MQ=255)
                                                                                                                         gCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGCCCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGggcgg                          >  2:153530/1‑149 (MQ=255)
                                                                                                                            gCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCa                       >  1:310766/1‑149 (MQ=255)
                                                                                                                               tCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACg                     >  3:305628/1‑148 (MQ=255)
                                                                                                                               tCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACAACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACg                     >  7:334365/1‑148 (MQ=255)
                                                                                                                                 gACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCAt                  <  4:347540/149‑1 (MQ=255)
                                                                                                                                    gcgGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCaa               <  8:283836/149‑1 (MQ=255)
                                                                                                                                     cgGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGAGTGACGGGATCAACTACTAGGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg              <  7:238348/149‑1 (MQ=255)
                                                                                                                                       gtgtTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGc            <  2:319049/149‑1 (MQ=255)
                                                                                                                                        cgtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg              <  1:336372/145‑1 (MQ=255)
                                                                                                                                        cgtTCGCCCTGGAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAg              >  2:336372/2‑146 (MQ=255)
                                                                                                                                             gCCCTCCACGGCGAACGGCTTACCCGCGACCCGCTCAGTGAAGTGGTACGAATAGAGCGTGACGGGATAAACTACTACGTCAAGCACAAAACCGGGGCCGGCAAGCACATGCGCCGGTGCCTGCGGCGGCCACGCATCAACGCCGAATg      >  4:354974/1‑149 (MQ=255)
                                                                                                                                              cccTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATgg     <  7:103523/149‑1 (MQ=255)
                                                                                                                                                  ccAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                           >  1:351964/1‑59 (MQ=255)
                                                                                                                                                  ccAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACg                                                                                           <  2:351964/59‑1 (MQ=255)
                                                                                                                                                   cAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCaa  <  7:166451/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
GCTCGCCTGGCTGGAGCGCAAGGCCCAGAAACTCTACGATCGCAAGCAACGCTATGGGGATGCACTCTGATGGCGAGTTGGACACTGGCGCCGGGGTACGAACACCTGGCGGCGGACTTCGGCAGCCTCGACGCGGTGTTCGCCCTCCAGGGCGAACGGCTGACCCGCGACCCGCTCAGTGAAGTGGTACGCATCGAGCGTGACGGGATCAACTACTACGTCAAGCGCTACACCGGGGCCGGCAAGCACATGCGCCGGTACCTGGGGCGGCCACGCATCAAGGCCGAATGGCAA  >  NC_002947/418149‑418442

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