Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,502,238 C→A 5.1% intergenic (+30/‑134) yccA → / → xanP FtsH‑mediated proteolysis substrate/modulator/xanthine permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,502,2380CA 5.1% 252.3 / 2.2 79intergenic (+30/‑134)yccA/xanPFtsH‑mediated proteolysis substrate/modulator/xanthine permease
Reads supporting (aligned to +/- strand):  ref base C (39/36);  new base A (2/2);  total (41/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 = 7.96e-01

GATCATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCA  >  NC_002947/4502091‑4502380
                                                                                                                                                    |                                                                                                                                               
gATCATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCCGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGGAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTcg                                                                                                                                                >  2:116567/1‑149 (MQ=255)
gATCATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTcg                                                                                                                                                >  7:224054/1‑149 (MQ=255)
gATCATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTcg                                                                                                                                                <  2:8465/149‑1 (MQ=255)
          ggtggCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACCAGGGGATTTATTGCTGTACCGCAATAAGCCCGCATCGGCGGGCttt                                                                                                                                      >  1:146940/1‑149 (MQ=255)
           gtggCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCtttt                                                                                                                                     <  5:64292/148‑1 (MQ=255)
            tggCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCtttttt                                                                                                                                   <  8:224054/149‑1 (MQ=255)
            tggCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCtttttt                                                                                                                                   <  2:111850/149‑1 (MQ=255)
                    gTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGgcgc                                                                                                                           >  4:92901/1‑149 (MQ=255)
                       aCTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTg                                                                                                                         <  6:212294/148‑1 (MQ=255)
                             tcatGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATgg                                                                                                                  <  6:332306/149‑1 (MQ=255)
                                tGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGtt                                                                                                               >  5:305441/1‑149 (MQ=255)
                                   cAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAg                                                                                                            <  6:305441/149‑1 (MQ=255)
                                    aaCCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGc                                                                                                           >  5:252321/1‑149 (MQ=255)
                                    aaCCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAAAAGCCCGCCGAAGCGGGcttcttgcg                                                                                                                                 <  5:154334/127‑4 (MQ=255)
                                    aaCCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAAAAGCCCGCCGAAGCGGGcttcttgcg                                                                                                                                 >  6:154334/1‑124 (MQ=255)
                                       cATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGc                                                                                                        <  5:51536/149‑1 (MQ=255)
                                        aTCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCa                                                                                                       <  3:266779/149‑1 (MQ=255)
                                         tCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAg                                                                                                      <  4:68639/149‑1 (MQ=255)
                                         tCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAg                                                                                                      <  3:84003/149‑1 (MQ=255)
                                            gCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATg                                                                                                                   <  8:252636/133‑1 (MQ=255)
                                            gCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATg                                                                                                                   >  7:252636/1‑133 (MQ=255)
                                            gCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGGACCGCAAGAAGCCCGGTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGGTGCAGGag                                                                                                   >  8:68912/1‑148 (MQ=255)
                                              cTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGc                                                                                                 <  8:330414/149‑1 (MQ=255)
                                                 tATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTc                                                                                               <  3:427098/148‑1 (MQ=255)
                                                  aTGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGt                                                                                             >  6:152422/1‑149 (MQ=255)
                                                  aTGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGt                                                                                             >  2:158156/1‑149 (MQ=255)
                                                     tggcGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGc                                                                                           <  2:132515/145‑1 (MQ=255)
                                                      gtCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTgg                                                                                                                               <  3:242930/111‑1 (MQ=255)
                                                      gtCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTgg                                                                                                                               >  4:242930/1‑111 (MQ=255)
                                                         gATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCa                                                                                      >  4:343741/1‑149 (MQ=255)
                                                           tCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATc                                                                                    >  6:404203/1‑149 (MQ=255)
                                                              aCAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCg                                                                                 >  5:277978/1‑149 (MQ=255)
                                                                aaCCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGtt                                                                               >  6:181543/1‑149 (MQ=255)
                                                                aaCCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAAAAGCCCGCCGAAGCGGGcttcttgcgg                                                                                                                                <  4:129175/100‑5 (MQ=37)
                                                                aaCCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAAAAGCCCGCCGAAGCGGGcttcttgcgg                                                                                                                                >  3:129175/1‑96 (MQ=37)
                                                                   gtGTTTGTCAGGCTGCTGCAATTGTTTGGCATCATGGGTCGGGACGACTGATTGTTTGCTGTAGCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGTATATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTtc                                                                            <  8:182645/148‑1 (MQ=255)
                                                                   cTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTtc                                                                            >  2:151302/1‑149 (MQ=255)
                                                                    tGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTc                                                                                       >  5:48923/1‑137 (MQ=255)
                                                                    tGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTc                                                                                       <  6:48923/137‑1 (MQ=255)
                                                                      tttGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTtct                                                                           <  5:404203/147‑1 (MQ=255)
                                                                          tCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTAcc                                                                      <  6:277978/148‑1 (MQ=255)
                                                                          tCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCg                                                                     >  4:200752/1‑149 (MQ=255)
                                                                           cAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGGGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGc                                                                    >  8:312684/1‑149 (MQ=255)
                                                                            aGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGc                                                                    >  7:259536/1‑148 (MQ=255)
                                                                              ccTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTAc                                                                       <  5:8704/143‑1 (MQ=255)
                                                                              ccTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTAc                                                                       >  6:8704/1‑143 (MQ=255)
                                                                                    gcAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCAt                                                                                     >  1:167672/1‑123 (MQ=255)
                                                                                    gcAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCAt                                                                                     <  2:167672/123‑1 (MQ=255)
                                                                                    gcAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCa                                                           >  7:7101/1‑149 (MQ=255)
                                                                                    gcAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCa                                                           >  5:56601/1‑149 (MQ=255)
                                                                                     cAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTAATCCCGTTTTCTCTACCGCTATCCGCCAt                                                          >  4:52553/1‑149 (MQ=255)
                                                                                      aaTTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCAt                                                          >  5:295840/1‑148 (MQ=255)
                                                                                          gTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCAccc                                                     >  8:61483/1‑149 (MQ=255)
                                                                                           tttGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTa                                                                        <  2:434070/129‑1 (MQ=255)
                                                                                           tttGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTa                                                                        >  1:434070/1‑129 (MQ=255)
                                                                                             tGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCg                                                   >  4:128346/1‑148 (MQ=255)
                                                                                                catcatGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCa                                                           >  8:225597/1‑137 (MQ=255)
                                                                                                catcatGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCa                                                           <  7:225597/137‑1 (MQ=255)
                                                                                                catcatGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGa                                               <  7:61483/149‑1 (MQ=255)
                                                                                                  tcatGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGGGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAAt                                             <  1:116567/149‑1 (MQ=255)
                                                                                                    atGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATgc                                           >  5:340919/1‑149 (MQ=255)
                                                                                                     tGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATgcg                                          <  5:152422/149‑1 (MQ=255)
                                                                                                     tGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGGTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATgcg                                          <  2:146940/149‑1 (MQ=255)
                                                                                                      gggTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGGTATCCGCCATCACCCCGTAGAATgcgc                                         >  6:45106/1‑149 (MQ=255)
                                                                                                                    tgattgatTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGGTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTttcttc                            >  2:400285/1‑149 (MQ=255)
                                                                                                                       ttgattgCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCgg                         >  8:127313/1‑149 (MQ=255)
                                                                                                                       ttgattgCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGATGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCgg                         >  6:413467/1‑149 (MQ=255)
                                                                                                                         gattgCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCggg                        <  8:353048/148‑1 (MQ=255)
                                                                                                                              cTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGggcggc                   <  6:340919/148‑1 (MQ=255)
                                                                                                                               tGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCtt                 >  1:133667/1‑149 (MQ=255)
                                                                                                                                 tACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGTTGTAGGATGCTCGTGTCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCttt                <  1:140717/148‑1 (MQ=255)
                                                                                                                                 tACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCttt                <  8:259536/148‑1 (MQ=255)
                                                                                                                                   ccGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCa              <  5:119968/148‑1 (MQ=255)
                                                                                                                                    cGCAAGAAGCCCGCTTCGGCGGGCTTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCttt                >  4:200154/1‑146 (MQ=255)
                                                                                                                                    cGCAAGAAGCCCGCTTCGGCGGGCTTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCa              <  3:200154/148‑1 (MQ=255)
                                                                                                                                    cGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAg             >  6:415477/1‑148 (MQ=255)
                                                                                                                                      cAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCaa          <  6:295840/149‑1 (MQ=255)
                                                                                                                                        agaagCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATg        <  3:128346/149‑1 (MQ=255)
                                                                                                                                            gCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGc     >  4:122822/1‑148 (MQ=255)
                                                                                                                                             cccGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTc   <  5:415477/149‑1 (MQ=255)
                                                                                                                                              ccGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGCCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTca  >  6:379375/1‑149 (MQ=255)
                                                                                                                                               cGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTca  >  4:388933/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
GATCATTCATGGTGGCGAACGTAACTACATCATGGCAACCATCAGCCTGTATGTGTCGATCTACAACCTGTTTGTCAGCCTGCTGCAATTGTTTGGCATCATGGGTCGTGACGACTGATTGATTGCTGTACCGCAAGAAGCCCGCTTCGGCGGGCTTTTTTGTGGGCGCTGGAAATGGGTTGAGCTGCAGGATGCTCGTGCCGTCATCCCGTTTTCTCTACCGCTATCCGCCATCACCCCGTAGAATGCGCTCCTTTTTTCTTCCGGGGCGGCTTTCAGCAATGAGCTCA  >  NC_002947/4502091‑4502380

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