Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 122,225 T→C 5.1% intergenic (+35/+26) katE → / ← PP_0116 hydroperoxidase/putative lipoprotein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947122,2250TC 5.1% 235.9 / 2.1 78intergenic (+35/+26)katE/PP_0116hydroperoxidase/putative lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (37/37);  new base C (2/2);  total (39/39)
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.34e-01

CAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCAGCC  >  NC_002947/122079‑122372
                                                                                                                                                  |                                                                                                                                                   
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                                                                                                              >  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                             >  3:123199/1‑149 (MQ=255)
                                                                                                                      cTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCa                             >  7:31751/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                                                <  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)
                                                                                                                                                tGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCAGc   <  2:144266/149‑1 (MQ=255)
                                                                                                                                                 gTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCAGc   >  2:235947/1‑148 (MQ=255)
                                                                                                                                                  tAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCAGcc  <  5:150351/148‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CAAGGGGCTGTTGCTCGGGGATGAGCAGAAGGCAGTGGATGCCTTTGTCAAAGCGGTAGAAGGCCACCGGGTGTGGGAGCGGGAGGCGGCGGCTGAAGCTGTTCCTGCCTGACAAGGCCTTTTCGCGGGTAAACCCCCTCCTGATGTAGGAGCGGGTTTACCCGCGAAGCGTTTACTGCTGATGCGGGACCAGCACCACCTGTGCCGGCAGCGAGCGCTGGATCTGATGCTTCTGCTTCAGCTCGAACCCACTGTCAATCTTGCGCACGCGCTTTTCCAGCAAAGTCTTCAGCC  >  NC_002947/122079‑122372

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: