Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 122,222 A→C 5.1% intergenic (+32/+29) katE → / ← PP_0116 hydroperoxidase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947122,2220AC 5.1% 245.0 / 4.0 78intergenic (+32/+29)katE/PP_0116hydroperoxidase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (38/36);  new base C (2/2);  total (40/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.46e-01

GACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTC  >  NC_002947/122077‑122365
                                                                                                                                                 |                                                                                                                                               
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                                                                                                                                     <  4:343240/149‑1 (MQ=255)
         cTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCggg                                                                                                                                     <  3:91337/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                                                                                                       >  7:107887/1‑149 (MQ=255)
                                       gATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCg                                                                                                       >  6:335653/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                                                                                            >  4:377902/1‑116 (MQ=255)
                                                                                   ggAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGcacc                                                                                            <  3:377902/116‑1 (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                      >  3:123199/1‑149 (MQ=255)
                                                                                                                        cTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCa                      >  7:31751/1‑149 (MQ=255)
                                                                                                                         ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCgg                                                                                                      <  1:175503/68‑1 (MQ=255)
                                                                                                                         ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAc                     >  6:186350/1‑149 (MQ=255)
                                                                                                                         ttttCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCAc                     >  5:385651/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                                         <  2:272296/124‑1 (MQ=255)
                                                                                                                              gcgGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAAc                                         >  1:272296/1‑124 (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)
                                                                                                                                           ctcctGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGt   >  6:146321/1‑149 (MQ=255)
                                                                                                                                              ctGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTc  <  6:385651/147‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                               
GACAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTC  >  NC_002947/122077‑122365

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: