Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,682,703 C→G 5.9% intergenic (+24/+20) crtZ → / ← PP_3247 beta‑carotene hydroxylase/bile acid/Na+ symporter family transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,682,7030CG 5.9% 622.7 / 21.2 187intergenic (+24/+20)crtZ/PP_3247beta‑carotene hydroxylase/bile acid/Na+ symporter family transporter
Reads supporting (aligned to +/- strand):  ref base C (82/94);  new base G (6/5);  total (88/99)
Fisher's exact test for biased strand distribution p-value = 7.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.68e-01

Only 100 of 189 total aligned reads displayed.

CCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCGTTGTTC  >  NC_002947/3682556‑3682852
                                                                                                                                                   |                                                                                                                                                     
CCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTg                                                                                                                                                     <  2:1429369‑M1/150‑1 (MQ=255)
 CGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGc                                                                                                                                                    <  2:878479‑M1/150‑1 (MQ=255)
   cgcgcACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCg                                                                                                                                                  <  2:487212/150‑1 (MQ=255)
    gcgcACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCgg                                                                                                                                                 <  2:3311125/150‑1 (MQ=255)
    gcgcACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCgg                                                                                                                                                 <  2:4564898/150‑1 (MQ=255)
    gcgcACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCgg                                                                                                                                                 <  2:3078412/150‑1 (MQ=255)
       cACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCAGGTGCCGGGcg                                                                                                                                              <  2:5150359/150‑2 (MQ=255)
        aCCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGACGCAGCAATTGCGCAGCAGGCGCAGCCCATCGTGACGCATCCCTGGCACGCCCGGCCCCTGCCGGGCag                                                                                                                                             >  2:4115867/1‑149 (MQ=255)
            cTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCggg                                                                                                                                         <  2:1882494/150‑1 (MQ=255)
                  ccAAGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAAt                                                                                                                                   <  2:163322/150‑1 (MQ=255)
                    aTGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCa                                                                                                                                 <  2:2699349/150‑1 (MQ=255)
                    aTGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCa                                                                                                                                 >  2:5331506/1‑150 (MQ=255)
                           gAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACg                                                                                                                          <  1:1845551/150‑1 (MQ=255)
                           gAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACg                                                                                                                          >  1:1092810/1‑150 (MQ=255)
                              ggggCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGtt                                                                                                                       >  2:1840632/1‑150 (MQ=255)
                                 gcgcACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCa                                                                                                                    <  1:4336231/150‑1 (MQ=255)
                                   gcACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCAcc                                                                                                                  >  2:652397/1‑150 (MQ=255)
                                    cACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCAccc                                                                                                                 >  1:698589/1‑150 (MQ=255)
                                     aCGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCt                                                                                                                >  2:2148475/1‑150 (MQ=255)
                                       ggGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCGTg                                                                                                              >  1:4092224‑M2/3‑150 (MQ=255)
                                        gGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCACTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTgg                                                                                                             <  1:3478032‑M2/149‑1 (MQ=255)
                                        gGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTgg                                                                                                             <  1:1429369‑M2/149‑1 (MQ=255)
                                         ggAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTggg                                                                                                            <  2:3679223/150‑1 (MQ=255)
                                          gAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGCAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGc                                                                                                           <  2:4061034/150‑1 (MQ=255)
                                             cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGtt                                                                                                                              <  1:1116992/128‑1 (MQ=255)
                                             cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCa                                                                                                        <  2:1944524/150‑1 (MQ=255)
                                             cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCa                                                                                                        >  1:5052458/1‑150 (MQ=255)
                                             cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTACCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCa                                                                                                        <  2:1325149/150‑1 (MQ=255)
                                                gtCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgcc                                                                                                     <  1:2917519/150‑1 (MQ=255)
                                                  gcggagggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgccgc                                                                                                   <  1:1951901/145‑1 (MQ=255)
                                                    ttCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCca                                                                                                 <  2:4455962/150‑1 (MQ=255)
                                                        ggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCTACATg                                                                                             >  1:2939506/1‑150 (MQ=255)
                                                        ggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATg                                                                                             >  2:128691/1‑150 (MQ=255)
                                                        ggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATg                                                                                             >  1:4244730/1‑150 (MQ=255)
                                                             ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcgc                                                                                        >  2:907514/1‑150 (MQ=255)
                                                             ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcgc                                                                                        >  2:2171904/1‑150 (MQ=255)
                                                                tATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTgg                                                                                                             >  1:9382/1‑126 (MQ=255)
                                                                tATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGc                                                                                     >  1:1092068/1‑150 (MQ=255)
                                                                tATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGc                                                                                     >  1:4027315/1‑150 (MQ=255)
                                                                  tGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTc                                                                                   >  2:2034887/1‑150 (MQ=255)
                                                                    cgccgccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTcac                                                                                 <  2:4410844/150‑1 (MQ=255)
                                                                       cgccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTcaccca                                                                              >  2:363207/1‑150 (MQ=255)
                                                                        gccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTggg                                                                                                            >  2:425767/1‑119 (MQ=255)
                                                                        gccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTcacccac                                                                             <  1:3758255/150‑1 (MQ=255)
                                                                          cgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACccacca                                                                           <  2:2399245/150‑1 (MQ=255)
                                                                            cTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGc                                                                         >  2:2527740/1‑150 (MQ=255)
                                                                               aaGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGGTCACGGTTTCACCCTAGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCaa                                                                      >  2:2991275/1‑150 (MQ=255)
                                                                                     ttGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCgg                                                                <  2:4432831/150‑1 (MQ=255)
                                                                                       gAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGcc                                                              >  1:3657280/1‑150 (MQ=255)
                                                                                        aaGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTcacccac                                                                             <  2:5210459/134‑1 (MQ=255)
                                                                                         agcagcAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCcgc                                                            <  2:4642834/150‑1 (MQ=255)
                                                                                          gcagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCATCCAAACCCGGCcgcg                                                           >  2:5317078/1‑150 (MQ=255)
                                                                                          gcagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCcgcg                                                           >  2:1147563/1‑150 (MQ=255)
                                                                                          gcagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCcgcg                                                           <  1:4177758/150‑1 (MQ=255)
                                                                                           cagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCgg                                                          >  1:5107551/1‑150 (MQ=255)
                                                                                            agcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGt                                                         >  1:339939/1‑150 (MQ=255)
                                                                                               aaTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTaaa                                                      <  2:992146/150‑1 (MQ=255)
                                                                                               aaTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTaaa                                                      <  2:339939/150‑1 (MQ=255)
                                                                                                aTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAg                                                     <  1:1049917/150‑1 (MQ=255)
                                                                                                  tGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGaa                                                   <  1:749031/150‑1 (MQ=255)
                                                                                                    cgcAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACa                                                 <  1:2867594/150‑1 (MQ=255)
                                                                                                    cgcAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACa                                                 <  1:1735300/150‑1 (MQ=255)
                                                                                                       agcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGa                                              >  1:266145/1‑150 (MQ=255)
                                                                                                       agcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGa                                              <  2:3145052/150‑1 (MQ=255)
                                                                                                        gcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCaaa                                                                     <  2:938619/126‑1 (MQ=255)
                                                                                                        gcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCaaa                                                                     >  1:938619/1‑126 (MQ=255)
                                                                                                         cagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATc                                            >  1:1068913/1‑150 (MQ=255)
                                                                                                           ggCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCaa                                          <  1:536390/150‑1 (MQ=255)
                                                                                                            gcgcAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAg                                         >  1:5249050/1‑150 (MQ=255)
                                                                                                            gcgcAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAg                                         <  1:2486468/150‑1 (MQ=255)
                                                                                                                aGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGt                                     <  2:3390607/150‑1 (MQ=255)
                                                                                                                aGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGt                                     <  1:198208/150‑1 (MQ=255)
                                                                                                                    aaTCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACaa                                 >  1:402749/1‑150 (MQ=255)
                                                                                                                     aTCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCAACGGGTGAAATCAGGTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCACACCCGGCCGCGGTACAGAACCTGATCAAGCTGTACAAc                                >  2:2865292/1‑150 (MQ=255)
                                                                                                                        gTGACGCATCCCTGGCTCGCCCGGCACCTGCCGGGCATCTGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCa                             <  2:3946064/150‑1 (MQ=255)
                                                                                                                        gTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCa                             <  2:837124/150‑1 (MQ=255)
                                                                                                                           aCGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGcc                          >  1:3289133/1‑150 (MQ=255)
                                                                                                                           aCGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGcc                          <  1:2395400/150‑1 (MQ=255)
                                                                                                                             gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCgg                        <  2:5107551/150‑1 (MQ=255)
                                                                                                                             gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCACCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCgg                        >  1:1837797/1‑150 (MQ=255)
                                                                                                                                tCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCAc                     <  1:172266/150‑1 (MQ=255)
                                                                                                                                 cccTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCAcc                    <  2:261295/150‑1 (MQ=255)
                                                                                                                                   acccgatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGACTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGc                  >  1:5229813/8‑150 (MQ=255)
                                                                                                                                    tGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGAt                                             >  1:4355585/1‑122 (MQ=255)
                                                                                                                                    tGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGcc                 >  2:5393005/1‑150 (MQ=255)
                                                                                                                                    cccgatgCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGcc                 >  1:91680/7‑150 (MQ=255)
                                                                                                                                     ccgatgCCCGGCAGGTGCCGGGCGTCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACg                               <  1:3257355/130‑1 (MQ=255)
                                                                                                                                     ccgatgCCCGGCAGGTGCCGGGCGTCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCa                >  1:2230215/6‑150 (MQ=255)
                                                                                                                                     ccgatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCa                >  2:1837977/6‑150 (MQ=255)
                                                                                                                                      gCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCTGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCAt               <  2:4271637/150‑1 (MQ=255)
                                                                                                                                      cgatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCAt               >  2:4353608/5‑150 (MQ=255)
                                                                                                                                        aCGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATgg             >  1:3361105/1‑150 (MQ=255)
                                                                                                                                        aCGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATgg             >  1:635162/1‑150 (MQ=255)
                                                                                                                                          gCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTc           <  1:5146757/150‑1 (MQ=255)
                                                                                                                                          gCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTc           <  1:128818/150‑1 (MQ=255)
                                                                                                                                            ccGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTAc                                   <  2:4765381/124‑1 (MQ=255)
                                                                                                                                            ccGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTAc                                   >  1:4765381/1‑124 (MQ=255)
                                                                                                                                              ggCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCgt       <  2:1357039/150‑1 (MQ=255)
                                                                                                                                               gCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAAAATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCgtt      <  1:941178/150‑1 (MQ=255)
                                                                                                                                                 aCCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCgttgt    <  2:1092193/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCGTTGTTC  >  NC_002947/3682556‑3682852

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