Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,682,702 C→G 5.9% intergenic (+23/+21) 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,7020CG 5.9% 621.7 / 21.2 187intergenic (+23/+21)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.47e-01

Only 100 of 190 total aligned reads displayed.

GCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCGTTGTT  >  NC_002947/3682553‑3682851
                                                                                                                                                     |                                                                                                                                                     
gCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCAc                                                                                                                                                       <  2:1399498/150‑1 (MQ=255)
    CGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGc                                                                                                                                                   <  2:878479‑M1/150‑1 (MQ=255)
      cgcgcACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCg                                                                                                                                                 <  2:487212/150‑1 (MQ=255)
       gcgcACCGCTTGCACCATGCGGTGTAGGGGCGCACGGGGAGAGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCgg                                                                                                                                                <  2:3311134/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)
           aCCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCATGCCCGGCAGGTGCCGGGcgt                                                                                                                                            <  1:3331233/150‑3 (MQ=255)
            ccGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCAGGTGCCGGGcgtg                                                                                                                                           <  1:947384/150‑4 (MQ=255)
            ccGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGACACGCCCGGCACCTGCCGGGCATc                                                                                                                                           >  1:5040752/1‑150 (MQ=255)
             cGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCg                                                                                                                                          >  1:1663777/1‑150 (MQ=255)
               cTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCggg                                                                                                                                        <  2:1882494/150‑1 (MQ=255)
                ttGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGt                                                                                                                                       >  1:5366324/1‑150 (MQ=255)
                ttGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGt                                                                                                                                       <  1:2363324/150‑1 (MQ=255)
                   cACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGaa                                                                                                                                    <  1:221417/150‑1 (MQ=255)
                       aTGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCa                                                                                                                                <  2:2699349/150‑1 (MQ=255)
                       aTGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCa                                                                                                                                >  2:5331506/1‑150 (MQ=255)
                         gCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGt                                                                                                                              <  1:3750074/150‑1 (MQ=255)
                             tGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCAc                                                                                                                          >  2:1247096/1‑150 (MQ=255)
                              gAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACg                                                                                                                         <  2:1548971/150‑1 (MQ=255)
                              gAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACg                                                                                                                         >  1:1092810/1‑150 (MQ=255)
                              gAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACg                                                                                                                         <  1:1845551/150‑1 (MQ=255)
                                 ggggCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGtt                                                                                                                      >  2:1840632/1‑150 (MQ=255)
                                      gcACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCAcc                                                                                                                 >  2:2456419/1‑150 (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)
                                           gGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTgg                                                                                                            <  1:1429369‑M2/149‑1 (MQ=255)
                                            ggAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTggg                                                                                                           <  2:3679223/150‑1 (MQ=255)
                                             gAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGCAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGc                                                                                                          <  2:4061034/150‑1 (MQ=255)
                                                cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGtt                                                                                                                             >  2:1116992/1‑128 (MQ=255)
                                                cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGtt                                                                                                                             <  1:1116992/128‑1 (MQ=255)
                                                cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCa                                                                                                       <  2:1944524/150‑1 (MQ=255)
                                                cGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTACCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCa                                                                                                       <  2:1325149/150‑1 (MQ=255)
                                                      gTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgccgcc                                                                                                 >  1:3465213/1‑150 (MQ=255)
                                                       ttCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCca                                                                                                <  2:4455962/150‑1 (MQ=255)
                                                           ggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCTACATg                                                                                            >  1:2939506/1‑150 (MQ=255)
                                                           ggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATg                                                                                            >  2:128818/1‑150 (MQ=255)
                                                           ggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATg                                                                                            >  1:4244701/1‑150 (MQ=255)
                                                            gTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGc                                                                                           <  1:1975230/150‑1 (MQ=255)
                                                                ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcgc                                                                                       >  2:2312928/1‑150 (MQ=255)
                                                                ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcgc                                                                                       >  2:2171708/1‑150 (MQ=255)
                                                                ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcgc                                                                                       >  2:907514/1‑150 (MQ=255)
                                                                   tATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTgg                                                                                                            >  1:9382/1‑126 (MQ=255)
                                                                   tATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGc                                                                                    >  1:1092068/1‑150 (MQ=255)
                                                                     tGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTc                                                                                  >  2:2034887/1‑150 (MQ=255)
                                                                          cgccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTcaccca                                                                             >  2:363207/1‑150 (MQ=255)
                                                                           gccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTggg                                                                                                           >  2:425767/1‑119 (MQ=255)
                                                                           gccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTggg                                                                                                           <  1:425767/119‑1 (MQ=255)
                                                                             cgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACccacca                                                                          <  1:2875428/150‑1 (MQ=255)
                                                                             cgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACccacca                                                                          <  1:1734221/150‑1 (MQ=255)
                                                                                  aaGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGGTCACGGTTTCACCCTAGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCaa                                                                     >  2:2991275/1‑150 (MQ=255)
                                                                                        ttGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCgg                                                               <  1:50645/150‑1 (MQ=255)
                                                                                        ttGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCgg                                                               <  2:4432831/150‑1 (MQ=255)
                                                                                           aaGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTcacccac                                                                            <  2:5210459/134‑1 (MQ=255)
                                                                                            agcagcAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCcgc                                                           <  2:4642834/150‑1 (MQ=255)
                                                                                             gcagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCATCCAAACCCGGCcgcg                                                          >  2:5317078/1‑150 (MQ=255)
                                                                                             gcagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCcgcg                                                          <  1:4177758/150‑1 (MQ=255)
                                                                                               agcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGt                                                        >  1:339938/1‑150 (MQ=255)
                                                                                               agcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGt                                                        >  1:1847141/1‑150 (MQ=255)
                                                                                                  aaTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTaaa                                                     <  2:992146/150‑1 (MQ=255)
                                                                                                  aaTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTaaa                                                     <  2:339938/150‑1 (MQ=255)
                                                                                                   aTTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAg                                                    <  1:1049917/150‑1 (MQ=255)
                                                                                                     tGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGaa                                                  <  1:1771057/150‑1 (MQ=255)
                                                                                                     tGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGaa                                                  <  1:749031/150‑1 (MQ=255)
                                                                                                        gcagcaGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACAt                                               >  1:402573/1‑150 (MQ=255)
                                                                                                          agcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGa                                             <  2:3145052/150‑1 (MQ=255)
                                                                                                          agcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGa                                             >  1:266146/1‑150 (MQ=255)
                                                                                                           gcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCaaa                                                                    <  2:938619/126‑1 (MQ=255)
                                                                                                           gcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCaaa                                                                    >  1:939024/1‑126 (MQ=255)
                                                                                                            cagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATc                                           >  1:1068913/1‑150 (MQ=255)
                                                                                                               gcgcAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAg                                        <  1:2486691/150‑1 (MQ=255)
                                                                                                               gcgcAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAg                                        >  1:5249050/1‑150 (MQ=255)
                                                                                                                       aaTCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACaa                                >  1:402749/1‑150 (MQ=255)
                                                                                                                        aTCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAAc                               <  1:363207/150‑1 (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:2395562/150‑1 (MQ=255)
                                                                                                                              aCGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGcc                         >  1:3289111/1‑150 (MQ=255)
                                                                                                                                gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCgg                       <  2:5107551/150‑1 (MQ=255)
                                                                                                                                    cccTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCAcc                   <  2:261295/150‑1 (MQ=255)
                                                                                                                                      acccgatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGACTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGc                 >  1:5229813/8‑150 (MQ=255)
                                                                                                                                       tGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGAt                                            <  2:4355585/122‑1 (MQ=255)
                                                                                                                                       tGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGcc                >  2:5393005/1‑150 (MQ=255)
                                                                                                                                       cccgatgCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGcc                >  1:91680/7‑150 (MQ=255)
                                                                                                                                        ccgatgCCCGGCAGGTGCCGGGCGTCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACg                              <  1:3257355/130‑1 (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:635162/1‑150 (MQ=255)
                                                                                                                                             gCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTc          <  1:5146757/150‑1 (MQ=255)
                                                                                                                                             gCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTc          <  1:128818/150‑1 (MQ=255)
                                                                                                                                               ccGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTcgc        <  1:4353608/150‑1 (MQ=255)
                                                                                                                                               ccGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTAc                                  <  2:4765381/124‑1 (MQ=255)
                                                                                                                                               ccGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTAc                                  >  1:4765270/1‑124 (MQ=255)
                                                                                                                                                 ggCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCgt      <  2:1357039/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCACCGCGCGCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCGTTGTT  >  NC_002947/3682553‑3682851

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: