Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 122,215 C→G 6.3% intergenic (+25/+36) katE → / ← PP_0116 hydroperoxidase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947122,2150CG 6.3% 234.8 / 3.8 74intergenic (+25/+36)katE/PP_0116hydroperoxidase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base C (36/33);  new base G (2/3);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAG  >  NC_002947/122069‑122363
                                                                                                                                                  |                                                                                                                                                    
tAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCcct                                                                                                                                                    >  3:351824/1‑149 (MQ=255)
    ggaggaCAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCcctcct                                                                                                                                                 >  5:80318/1‑148 (MQ=255)
        gaCAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATg                                                                                                                                             >  5:19326/1‑148 (MQ=255)
        gaCAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATg                                                                                                                                             >  3:294664/1‑148 (MQ=255)
          cAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAg                                                                                                                                          >  7:333018/1‑149 (MQ=255)
             ggggCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAg                                                                                                                                       >  3:312120/1‑149 (MQ=255)
             ggggCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAg                                                                                                                                       >  1:54715/1‑149 (MQ=255)
             ggggCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAg                                                                                                                                       >  2:17549/1‑149 (MQ=255)
               ggCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCg                                                                                                                                     >  4:91337/1‑149 (MQ=255)
                 cTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCggg                                                                                                                                   <  3:91337/149‑1 (MQ=255)
                 cTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCggg                                                                                                                                   <  4:343240/149‑1 (MQ=255)
                    ttGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAATCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTAGGGGGTAAACCCCCTCCTGATGTAGGAGCGGGt                                                                                                                                  <  8:214348/147‑1 (MQ=255)
                    ttGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGt                                                                                                                                  <  3:176438/147‑1 (MQ=255)
                       cTCGGTGATGATCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGTGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTAc                                                                                                                              <  2:266968/148‑1 (MQ=255)
                       cTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTAcc                                                                                                                             >  8:15093/1‑149 (MQ=255)
                          ggggATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTAccc                                                                                                                            <  2:307536/147‑1 (MQ=255)
                           gggATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCcgc                                                                                                                          >  3:249606/1‑148 (MQ=255)
                             gATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCGTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGaa                                                                                                                       <  5:27801/149‑1 (MQ=255)
                                  gCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGtt                                                                                                                  <  7:318439/149‑1 (MQ=255)
                                    aGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTTATGTAGGAGCGGGTTTACCCGCGAAGCGTTTa                                                                                                                <  7:4958/149‑1 (MQ=255)
                                    aGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTa                                                                                                                <  6:80318/149‑1 (MQ=255)
                                       aGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTActg                                                                                                             >  2:113860/1‑149 (MQ=255)
                                           aGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTAc                                                                                                               >  7:273728/1‑143 (MQ=255)
                                           aGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTAc                                                                                                               <  8:273728/143‑1 (MQ=255)
                                               gATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCg                                                                                                     >  6:335653/1‑149 (MQ=255)
                                               gATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCg                                                                                                     >  7:107887/1‑149 (MQ=255)
                                                    cTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                                                                 >  8:152453/1‑148 (MQ=255)
                                                     tttGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGAc                                                                                                 >  4:156977/1‑147 (MQ=255)
                                                        gTCAAAGCGGTTCATGGCAACCGGGTGTGGGAGCGGGAGGCGGCGGGTGAAGCTGTTCCTGCCTGTCAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACGCGCGAAGCGTTTACTGCTGATGCGGGACCAGc                                                                                             <  5:39506/148‑1 (MQ=255)
                                                          cAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGGGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcac                                                                                           >  1:125367/1‑148 (MQ=255)
                                                           aaaGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcacca                                                                                         >  7:341097/1‑149 (MQ=255)
                                                             aGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCaccacc                                                                                       <  6:19326/149‑1 (MQ=255)
                                                              gCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCt                                                                                      <  4:249606/149‑1 (MQ=255)
                                                                  tAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGc                                                                                  <  4:312120/149‑1 (MQ=255)
                                                                         ccACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCagc                                                                           <  3:313584/149‑1 (MQ=255)
                                                                          cACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCagcg                                                                          >  7:202574/1‑149 (MQ=255)
                                                                           aCCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCagcga                                                                         <  4:351824/149‑1 (MQ=255)
                                                                           aCCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCagc                                                                           >  4:219206/1‑147 (MQ=255)
                                                                                gtgtGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCt                                                                    <  2:54715/149‑1 (MQ=255)
                                                                                   tgGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGTGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTgg                                                                  <  4:294664/148‑1 (MQ=255)
                                                                                      gAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCt                                                              <  3:219206/149‑1 (MQ=255)
                                                                                           ggAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcacc                                                                                          <  3:377902/116‑1 (MQ=255)
                                                                                           ggAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcacc                                                                                          >  4:377902/1‑116 (MQ=255)
                                                                                             aGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgctt                                                       <  5:137067/149‑1 (MQ=255)
                                                                                                 ggcggcTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttctgc                                                   <  1:113860/149‑1 (MQ=255)
                                                                                                   cggcTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttctgctt                                                 <  7:152453/149‑1 (MQ=255)
                                                                                                    ggcTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGAtgcttctgcttc                                                <  1:17549/149‑1 (MQ=255)
                                                                                                       tGAAGCTGTTCCTGCCTGACAAGGACTTTTCGCGGGTAAACACCCTCCTGAAGTAGGAGCGGGTTTACCCCCGAATCGTTTAATGCTGATTCCGGACCAGCAACCCCAGTGCCGGCAGCCAGCGCATGAtatgattcatatgattcat                                              >  6:40758/1‑147 (MQ=255)
                                                                                                         aaGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTTCCCGCGAAGCGTTTACTGCTGATGCGGGTCCAGCCCCACCAGTGCCGGCAGCGAGCGCTGGATCTGATTGTTCTGGTTACGCTc                                           >  4:347368/1‑149 (MQ=255)
                                                                                                          aGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTc                                           >  3:67346/1‑148 (MQ=255)
                                                                                                           gCTGTTCGTGCCTGTGATGGTGTTTTGGCGGGTTTACGCGCTCGTGATGTTGGTGCGGGTTTACCCGCGTAGCGTTTAGTGGTGATGGGGGAGGAGCACCAGCTGTGGCGGCTGCGAGCGCTGGATCTGATGCTTCTGCTTCAGGTCg                                          <  4:67346/148‑1 (MQ=255)
                                                                                                              gTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAAcc                                      >  7:247287/1‑149 (MQ=255)
                                                                                                               ttCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAAcc                                      >  4:49685/1‑148 (MQ=255)
                                                                                                                   tgcctgACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTg                                 >  6:150351/1‑149 (MQ=255)
                                                                                                                    gcctgACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGt                                >  1:138072/1‑149 (MQ=255)
                                                                                                                     cctgACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGt                                >  7:358892/1‑148 (MQ=255)
                                                                                                                         aCAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAAt                            <  2:125367/148‑1 (MQ=255)
                                                                                                                         aCAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATc                           >  1:144266/1‑149 (MQ=255)
                                                                                                                         aCAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATc                           >  3:103197/1‑149 (MQ=255)
                                                                                                                           aaGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCtt                         <  7:15093/149‑1 (MQ=255)
                                                                                                                             ggCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTgc                       <  8:341097/149‑1 (MQ=255)
                                                                                                                                cTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCa                    >  7:31751/1‑149 (MQ=255)
                                                                                                                                cTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCa                    >  3:123199/1‑149 (MQ=255)
                                                                                                                                 ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCgg                                                                                                    <  1:175503/68‑1 (MQ=255)
                                                                                                                                 ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAc                   >  5:385651/1‑149 (MQ=255)
                                                                                                                                 ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAc                   >  6:186350/1‑149 (MQ=255)
                                                                                                                                    tCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAcgc                 <  3:156977/148‑1 (MQ=255)
                                                                                                                                     cgcgGGTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAcgc                 <  3:311616/147‑1 (MQ=38)
                                                                                                                                      gcgGGTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTgc                       <  4:310924/140‑1 (MQ=38)
                                                                                                                                      gcgGGTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTgc                       >  3:310924/1‑140 (MQ=38)
                                                                                                                                      gcgGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAAc                                       >  1:272296/1‑124 (MQ=255)
                                                                                                                                      gcgGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAAc                                       <  2:272296/124‑1 (MQ=255)
                                                                                                                                        gggTAAACCCGCTCCTACATCAGGAGGGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCAACTGTCAATCTTGCGCACGCGCttt            >  6:50296/1‑149 (MQ=38)
                                                                                                                                          gTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTc          <  3:347368/149‑1 (MQ=255)
                                                                                                                                             aaCCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAGCCCACTGTCAATCTTGCGCACGCGCTTTTCCAg       >  8:187078/1‑149 (MQ=255)
                                                                                                                                             aaCCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCa        <  4:103197/148‑1 (MQ=255)
                                                                                                                                                ccccTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCaa    <  8:107887/149‑1 (MQ=255)
                                                                                                                                                  cctcctGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAg  >  8:239248/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TAAAGGAGGACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAG  >  NC_002947/122069‑122363

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: