Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,680,738 C→G 8.3% 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,680,7380CG 8.3% 220.7 / 9.8 72intergenic (+24/+20)crtZ/PP_3247beta‑carotene hydroxylase/bile acid/Na+ symporter family transporter
Reads supporting (aligned to +/- strand):  ref base C (36/30);  new base G (4/2);  total (40/32)
Fisher's exact test for biased strand distribution p-value = 6.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.59e-01

GCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCGTTG  >  NC_002947/3680597‑3680884
                                                                                                                                              |                                                                                                                                                  
gcACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCggg                                                                                                                                             <  5:313128/149‑1 (MQ=255)
   ccGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCAGGTGCCGGGcgt                                                                                                                                          <  3:20370/149‑3 (MQ=255)
       ttGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCggg                                                                                                                                      >  7:84224/1‑149 (MQ=255)
         gCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCATGCCCGGCAGGTGCCGGGcgtgcca                                                                                                                                      <  8:228512/147‑7 (MQ=255)
         gCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGt                                                                                                                                     <  4:262759/148‑1 (MQ=255)
           aCCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAACTCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGa                                                                                                                                   <  1:75718/149‑1 (MQ=255)
           aCCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGa                                                                                                                                   >  8:75688/1‑148 (MQ=255)
           aCCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGa                                                                                                                                   >  5:36886/1‑148 (MQ=255)
            ccATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGaaa                                                                                                                                 <  8:84224/149‑1 (MQ=255)
                   gTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCACAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGtt                                                                                                                           >  4:144872/1‑148 (MQ=255)
                       aGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCCTCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACgg                                                                                                                      >  8:83128/1‑149 (MQ=255)
                                  gggAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTg                                                                                                           <  2:245172/149‑1 (MQ=255)
                                   ggAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTgg                                                                                                          <  6:36886/149‑1 (MQ=255)
                                     aGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGc                                                                                                        <  3:295102/149‑1 (MQ=255)
                                        gtgtCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCt                                                                                                       >  3:223799/1‑147 (MQ=255)
                                         tgtCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAg                                                                                                    >  3:15789/1‑149 (MQ=255)
                                             gTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCgg                                                                                                                                              >  4:286094/1‑103 (MQ=255)
                                             gTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCgg                                                                                                                                              <  3:286094/103‑1 (MQ=255)
                                                cGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgc                                                                                                   >  7:274592/1‑143 (MQ=255)
                                                cGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgc                                                                                                   <  8:274592/143‑1 (MQ=255)
                                                 gggTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCcaca                                                                                            >  2:272410/1‑149 (MQ=255)
                                                   gTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCa                                                                                                                              <  7:157975/113‑1 (MQ=255)
                                                   gTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCa                                                                                                                              >  8:157975/1‑113 (MQ=255)
                                                   gTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATg                                                                                          <  2:156153/149‑1 (MQ=255)
                                                       ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCCCCCTTGGGCTCAGCCGCCACATGCCcg                                                                                      >  5:64525/1‑149 (MQ=255)
                                                       ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcg                                                                                      >  5:269774/1‑149 (MQ=255)
                                                       ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcg                                                                                      >  5:82225/1‑149 (MQ=255)
                                                       ttctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcg                                                                                      >  5:59621/1‑149 (MQ=255)
                                                         ctATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCcgcg                                                                                    >  1:197959/1‑149 (MQ=255)
                                                              cgccgccgcTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTc                                                                                <  4:223799/148‑1 (MQ=255)
                                                              cgccgccgcTGAAGGTGTTGAAGCAGCAATTGCACAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTc                                                                                <  3:144872/148‑1 (MQ=255)
                                                                              gtTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAAccc                                                               >  3:140185/1‑149 (MQ=255)
                                                                               ttGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCg                                                              >  1:270500/1‑149 (MQ=255)
                                                                                 gAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGc                                                            >  1:143997/1‑149 (MQ=255)
                                                                                    gcagcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCcgc                                                         >  6:97484/1‑149 (MQ=255)
                                                                                       gcaATTGCGCAGCAGGCGCAGCCAATCGTGACGCATACCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCgg                                                       >  8:34881/1‑148 (MQ=255)
                                                                                           ttGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAg                                                  <  4:15789/149‑1 (MQ=255)
                                                                                              cgcAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGaa                                                <  7:34881/148‑1 (MQ=255)
                                                                                                 agcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATg                                            <  6:82225/149‑1 (MQ=255)
                                                                                                 agcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATg                                            <  8:150269/149‑1 (MQ=255)
                                                                                                  gcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGa                                           <  1:272410/149‑1 (MQ=255)
                                                                                                  gcagGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGa                                           <  1:182650/149‑1 (MQ=255)
                                                                                                    agGAGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGTCAAACCCGGCCGCGGTAAAGAACATGAt                                          >  8:36524/1‑148 (MQ=255)
                                                                                                        gcAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGC                                     >  8:188904‑M1/1‑149 (MQ=255)
                                                                                                        gcAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGC                                     <  2:270500‑M1/149‑1 (MQ=255)
                                                                                                          aGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTG                                   <  4:140185‑M1/149‑1 (MQ=255)
                                                                                                              aaTCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATCATCAAGCTGTAc                                >  2:75778/1‑148 (MQ=255)
                                                                                                                   gTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACg                            <  5:97484/148‑1 (MQ=255)
                                                                                                                    tGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCa                          <  7:140413/149‑1 (MQ=255)
                                                                                                                     gACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAg                         >  3:212850/1‑149 (MQ=255)
                                                                                                                       cGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGcc                       <  6:269774/149‑1 (MQ=255)
                                                                                                                        gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGaacaacgcgacca                      >  1:314051‑M1/1‑136 (MQ=255)
                                                                                                                        gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCg                      >  4:23468/1‑149 (MQ=255)
                                                                                                                        gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTG                                   <  8:162595‑M1/136‑1 (MQ=255)
                                                                                                                        gCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTG                                   >  7:162595‑M1/1‑136 (MQ=255)
                                                                                                                            cccTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCCGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAATCCCGGCCGCGGTAAAGATCATGAt                                          <  8:209293/125‑1 (MQ=255)
                                                                                                                             ccTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCAcc                 >  4:206698/1‑149 (MQ=255)
                                                                                                                              cTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGaacaacgcgaccatggcgg                >  4:182875‑M1/1‑130 (MQ=255)
                                                                                                                                ggCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGcc                       >  3:302823/1‑140 (MQ=255)
                                                                                                                                ggCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGcc                       <  4:302823/140‑1 (MQ=255)
                                                                                                                                ccgatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGc               >  2:331505/6‑148 (MQ=255)
                                                                                                                                ccgatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGc               >  7:106345/6‑148 (MQ=255)
                                                                                                                                 gCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgccgcc                                                                                               >  7:256954/1‑67 (MQ=255)
                                                                                                                                 gCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAgccgcc                                                                                               <  8:256954/67‑1 (MQ=255)
                                                                                                                                  gatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCAt            >  8:197051/4‑149 (MQ=255)
                                                                                                                                  gatgCCCGGCAGGTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCAt            >  8:202386/4‑149 (MQ=255)
                                                                                                                                   aCGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATg           >  6:290628/1‑149 (MQ=255)
                                                                                                                                       ccGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAACCAGtacaacgcgaccatcgcggtgccggctg       >  8:155917‑M1/1‑121 (MQ=255)
                                                                                                                                       ccGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCGCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCGATGGTcg       <  7:83128/149‑1 (MQ=255)
                                                                                                                                        cGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTcgc      >  8:55409/1‑149 (MQ=255)
                                                                                                                                            aCCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAATAACATGATCAAGCTGtacaacgcgaccatggcggtgccggctgcgttg  <  7:155917‑M1/149‑34 (MQ=255)
                                                                                                                                            aCCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCgttg  >  3:223566/1‑149 (MQ=255)
                                                                                                                                              |                                                                                                                                                  
GCACCGCTTGCACCATGCGGTGAAGGGGCGCACGGGGAGCGTGTCGTTCGGGTTCTTCTATGCGCCGCCGCTGAAGGTGTTGAAGCAGCAATTGCGCAGCAGGCGCAGCCAATCGTGACGCATCCCTGGCACGCCCGGCACCTGCCGGGCATCGGGTGAAATCAGTTCACGGTTTCACCCTTGGGCTCAGCCGCCACATGCCCGCGGCTCACCCACCAGCCAAACCCGGCCGCGGTAAAGAACATGATCAAGCTGTACAACGCAGCCGGCACCGCCATGGTCGCGTTG  >  NC_002947/3680597‑3680884

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: