Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,439,047 C→G 5.7% intergenic (+63/‑32) PP_3058 → / → PP_3059 hypothetical protein/tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,439,0470CG 5.7% 227.5 / 3.2 70intergenic (+63/‑32)PP_3058/PP_3059hypothetical protein/tail sheath protein
Reads supporting (aligned to +/- strand):  ref base C (32/34);  new base G (2/2);  total (34/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 = 2.83e-01

GCCACTCATGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGG  >  NC_002947/3438901‑3439191
                                                                                                                                                  |                                                                                                                                                 
gcCACTCATGACGAGGCCGAAGCTATTGCCGACGAGGTCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGcccc                                                                                                                                                 >  2:201647/1‑149 (MQ=255)
       aTGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCGGCCCcgcagaaa                                                                                                                                          <  1:10139/149‑6 (MQ=255)
       aTGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCGGCCCcgcagaaa                                                                                                                                          >  2:10139/1‑144 (MQ=255)
       aTGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCgg                                                                                                                                           <  5:147530/148‑1 (MQ=255)
       aTGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGc                                                                                                                                         >  7:11434/1‑149 (MQ=255)
       aTGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGc                                                                                                                                         >  2:27114/1‑149 (MQ=255)
       aTGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCGTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGc                                                                                                                                         >  4:218111/1‑149 (MQ=255)
          aCGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGccc                                                                                                                                       <  8:134426/148‑1 (MQ=255)
           cgtggCGGAAGATATTGCCGACGAGGCCATTGCGGTGATCAACGCAAGGGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTGCATCGTGGGGTTTTCTTTTCTGCGGGGCCGCGCGGGGCGGccgc                                                                                                                                      <  4:184322/145‑3 (MQ=255)
            gAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTg                                                                                                                                    >  6:130006/1‑149 (MQ=255)
            gAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTg                                                                                                                                    >  8:172617/1‑149 (MQ=255)
                ccGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCtt                                                                                                                                <  6:4626/149‑1 (MQ=255)
                     gatATTGCCGAAGAGGCAATTGCCGTGATCAACGCAAGAGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCTCCTTCCTTTAGCGCCCCGTCCATCGTGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTg                                                                                                                           <  5:355686/147‑1 (MQ=255)
                           gCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCcgcgc                                                                                                                                             <  4:41612/126‑1 (MQ=255)
                           gCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCcgcgc                                                                                                                                             >  3:41612/1‑126 (MQ=255)
                           gCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGcc                                                                                                                      <  3:93585/148‑1 (MQ=255)
                             cgacgaGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCca                                                                                                                    >  2:97644/1‑148 (MQ=255)
                                   ggCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACt                                                                                                             <  5:130006/149‑1 (MQ=255)
                                    gCCATTGCCGGGATCAACGCAAGCGGGGTCGGGGGCGGCAACGGCTGACCCGGAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGGGTTTTCGTTTCTGCGGGGCCGCCCCGCGGGGACCCTTGTCTTATCTGGAGCCCACATGACt                                                                                                             <  1:27114/149‑1 (MQ=255)
                                      cATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGa                                                                                                          <  1:97644/149‑1 (MQ=255)
                                           ccGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCtt                                                                                                     >  4:138688/1‑149 (MQ=255)
                                           ccGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCtt                                                                                                     <  3:138688/149‑1 (MQ=255)
                                             gTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTc                                                                                                   <  4:5430/149‑1 (MQ=255)
                                                aTCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACg                                                                                                <  7:172617/149‑1 (MQ=255)
                                                 tCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACgg                                                                                               >  5:7486/1‑149 (MQ=255)
                                                   aaCGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCca                                                                                                                    <  5:255253/126‑1 (MQ=255)
                                                   aaCGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCca                                                                                                                    >  6:255253/1‑126 (MQ=255)
                                                   aaCGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCg                                                                                             >  5:90611/1‑149 (MQ=255)
                                                     cGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTa                                                                                          >  4:212197/1‑149 (MQ=255)
                                                      gCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGc                                                                                              <  5:145206/145‑1 (MQ=255)
                                                      gCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGc                                                                                              >  6:145206/1‑145 (MQ=255)
                                                        aaGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCg                                                                                       <  3:212197/149‑1 (MQ=255)
                                                          gCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGt                                                                                      <  8:11434/149‑1 (MQ=255)
                                                            gTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTaa                                                                                    <  6:7486/149‑1 (MQ=255)
                                                                  ggTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTaacgaacg                                                                              >  4:104064/1‑149 (MQ=255)
                                                                    tgtCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTc                                                                            <  1:201647/149‑1 (MQ=255)
                                                                        ggCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCg                                                                                             <  6:72958/128‑1 (MQ=255)
                                                                        ggCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCg                                                                                             >  5:72958/1‑128 (MQ=255)
                                                                           aaCGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCg                                                                     >  3:268520/1‑149 (MQ=255)
                                                                           aaCGCCTGACCCGCAACCCGCCTTCCTTGACCGCCCCGTCCATCGGGGCGTTTTCGTATCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGCATTCTATGACGGCGTTACCGAAACGAACCTCGACAccc                                                                     >  4:133492/1‑148 (MQ=255)
                                                                            aCGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCgg                                                                    >  4:26886/1‑149 (MQ=255)
                                                                            aCGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCgg                                                                    <  5:167571/149‑1 (MQ=255)
                                                                             cGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCggg                                                                   <  3:104064/149‑1 (MQ=255)
                                                                                cTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGgcg                                                                >  5:185513/1‑149 (MQ=255)
                                                                                cTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGgcg                                                                >  6:42134/1‑149 (MQ=255)
                                                                                  gACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTaacgaa                                                                                >  4:26302/1‑131 (MQ=255)
                                                                                     ccGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgcag                                                           >  8:279975/1‑149 (MQ=255)
                                                                                       gCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgcagca                                                         <  6:185513/149‑1 (MQ=255)
                                                                                            ccGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCAtt                                                       <  4:274285/146‑1 (MQ=255)
                                                                                            ccGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCAtt                                                       >  3:274285/1‑146 (MQ=255)
                                                                                            ccGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTgcg                                                    <  3:133492/149‑1 (MQ=255)
                                                                                            ccGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGCAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTgcg                                                    <  7:65987/149‑1 (MQ=255)
                                                                                                cttccttGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                                                            <  4:61568/137‑1 (MQ=255)
                                                                                                cttccttGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                                                            <  4:56850/137‑1 (MQ=255)
                                                                                                cttccttGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                                                            >  3:56850/1‑137 (MQ=255)
                                                                                                cttccttGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                                                            >  3:61568/1‑137 (MQ=255)
                                                                                                 ttccttGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTc                                                                            <  5:53399/120‑1 (MQ=255)
                                                                                                 ttccttGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTc                                                                            >  6:53399/1‑120 (MQ=255)
                                                                                                            ccccGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGa                                    >  7:329828/1‑149 (MQ=255)
                                                                                                             cccGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGAt                                   >  6:364759/1‑149 (MQ=255)
                                                                                                                        gggCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                                                            >  3:11339/1‑113 (MQ=255)
                                                                                                                        gggCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCgca                                                            <  4:11339/113‑1 (MQ=255)
                                                                                                                         ggCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGt                       <  3:26886/149‑1 (MQ=255)
                                                                                                                              tttCGTTTCTGCGGGGCCGCCCCGCGCGGCGCTTGTCTTATCTGGAGCCCACATGACTGGAGTCTTTAACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCGTCGATCAGTGGCTTGGTCGaca                  >  6:273703/1‑149 (MQ=255)
                                                                                                                                  gTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGAGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCtt              >  8:321473/1‑149 (MQ=255)
                                                                                                                                     tCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCa            <  5:42134/148‑1 (MQ=255)
                                                                                                                                       acaagggCCGCGCGGGGCGGCCCCTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCg         >  1:162631/5‑149 (MQ=255)
                                                                                                                                         cGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAg       <  7:279975/149‑1 (MQ=255)
                                                                                                                                         aagggCCGCGCGGGGCGGCCCCTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGaca                  <  8:277782/136‑1 (MQ=255)
                                                                                                                                         aagggCCGCGCGGGGCGGCCCCTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGaca                  >  7:277782/3‑138 (MQ=255)
                                                                                                                                          ggggCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAgg      <  8:329828/149‑1 (MQ=255)
                                                                                                                                              ccgccCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCgg  <  6:90611/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
GCCACTCATGACGAGGCCGAAGCTATTGCCGACGAGGCCATTGCCGTGATCAACGCAAGCGTGGTCGGTGTCGGCAACGCCTGACCCGCAACCCGCCTTCCTTGAGCGCCCCGTCCATCGGGGCGTTTTCGTTTCTGCGGGGCCGCCCCGCGCGGCCCTTGTCTTATCTGGAGCCCACATGACTGGATTCTTTCACGGCGTTACCGTAACGAACGTCGACACCGGGGCGCGCAGCATTGCGCTGCCGTCCTCCTCGATCATTGGCTTGGTCGACACCTTCACCGAGGGCGG  >  NC_002947/3438901‑3439191

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: