Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,669,674 C→A 11.7% intergenic (+41/+34) prpD → / ← acnB 2‑methylcitrate dehydratase/aconitate hydratase B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,669,6740CA11.7% 223.2 / 16.6 77intergenic (+41/+34)prpD/acnB2‑methylcitrate dehydratase/aconitate hydratase B
Reads supporting (aligned to +/- strand):  ref base C (38/30);  new base A (4/5);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 7.24e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

CGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTCGCG  >  NC_002947/2669529‑2669822
                                                                                                                                                  |                                                                                                                                                    
cGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGc                                                                                                                                                   >  2:333160/1‑149 (MQ=255)
  cTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCg                                                                                                                                                 >  5:282934/1‑149 (MQ=255)
     tGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCggg                                                                                                                                               <  4:209047/148‑1 (MQ=255)
         cGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCAGCGATGCCGGGGCttt                                                                                                                                          <  7:117331/149‑1 (MQ=255)
                  tgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCggg                                                                                                                                               >  1:92352/1‑135 (MQ=255)
                  tgtCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCggg                                                                                                                                               <  2:92352/135‑1 (MQ=255)
                      agcagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTg                                                                                                                              >  4:176512/1‑148 (MQ=255)
                        cagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                           <  1:318845/149‑1 (MQ=255)
                        cagGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTaa                                                                                                                           <  3:188644/149‑1 (MQ=255)
                           gTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTc                                                                                                                        <  6:282934/149‑1 (MQ=255)
                            ttttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                             >  7:214760/1‑127 (MQ=255)
                            ttttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                             <  8:214760/127‑1 (MQ=255)
                            ttttCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACt                                                                                                                         >  3:140744/1‑147 (MQ=255)
                                 gAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTa                                                                                                                  >  3:125211/1‑149 (MQ=255)
                                 gAACTGTGAAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTa                                                                                                                  >  4:215659/1‑149 (MQ=255)
                                       tgCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCt                                                                                                             <  5:351121/148‑1 (MQ=255)
                                       tgCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCt                                                                                                             >  6:351121/1‑148 (MQ=255)
                                       tacAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                             <  2:44280/114‑1 (MQ=255)
                                       tacAAGGACCAGCAGAAGCTGGAAGGCACGGCAGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGc                                                                                                                                             >  1:44280/3‑116 (MQ=255)
                                         cAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGaa                                                                                                          >  2:48202/1‑149 (MQ=255)
                                               ccAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCATGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCg                                                                                                    <  6:115689/149‑1 (MQ=255)
                                                cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGtt                                                                                                                                       <  6:171708/113‑1 (MQ=255)
                                                cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGtt                                                                                                                                       >  5:171708/1‑113 (MQ=255)
                                                cagcagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCgg                                                                                                   >  3:142620/1‑149 (MQ=255)
                                                   cagAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGAt                                                                                                >  8:175516/1‑149 (MQ=255)
                                                    agAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATc                                                                                               >  6:231367/1‑149 (MQ=255)
                                                       aGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTg                                                                                            <  3:215659/149‑1 (MQ=255)
                                                            gAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGt                                                                                        <  2:220048/148‑1 (MQ=255)
                                                                  aCGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACCACCTGTAACTCCGCTTACGCCTGACCCACCGGGATCTTGGCGTTGGCTGc                                                                                 >  8:74559/1‑149 (MQ=255)
                                                                  aCGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCa                                                                                                                       >  5:365998/1‑111 (MQ=255)
                                                                  aCGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCa                                                                                                                       <  6:365998/111‑1 (MQ=255)
                                                                  aCGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCCGAACCACCGGGATCTTGGCGTTGGCTGc                                                                                 >  5:120217/1‑149 (MQ=255)
                                                                    ggCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTgg                                                                                     >  4:158044/1‑143 (MQ=255)
                                                                    ggCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTgg                                                                                     <  3:158044/143‑1 (MQ=255)
                                                                       tgtgCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCtt                                                                            >  6:146836/1‑149 (MQ=255)
                                                                          gCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgc                                                                         >  2:174941/1‑149 (MQ=255)
                                                                           cacaGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTcgcg                                                                        <  3:176512/149‑1 (MQ=255)
                                                                                   cGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCg                                                                 <  3:65621/148‑1 (MQ=255)
                                                                                    gTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGc                                                               <  5:146836/149‑1 (MQ=255)
                                                                                         cTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCAGCGATCt                                                          >  7:305257/1‑149 (MQ=255)
                                                                                          ttgtttgtGATATGAAGCGGCCGGCACAGGCCCAATACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCt                                                          >  1:234540/1‑149 (MQ=255)
                                                                                            gtttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGTGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTg                                                         <  1:48202/147‑1 (MQ=255)
                                                                                             tttgtGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTc                                                      >  8:333193/1‑149 (MQ=255)
                                                                                               tgtgATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGAAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTCGCTGCCACATCGCGGAACTCGGCGATATGGTCg                                                     >  2:278903/1‑148 (MQ=255)
                                                                                                 tgATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAg                                                  <  4:125211/149‑1 (MQ=255)
                                                                                                  gATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                       >  5:254450/1‑143 (MQ=255)
                                                                                                  gATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGt                                                       <  6:254450/143‑1 (MQ=255)
                                                                                                     atGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCt                                                >  6:210410/1‑147 (MQ=255)
                                                                                                       gAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAgg                                            >  7:313535/1‑149 (MQ=255)
                                                                                                         aGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTa                                          <  1:333160/149‑1 (MQ=255)
                                                                                                             gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCtt                                                                                             >  8:87946/1‑94 (MQ=255)
                                                                                                             gccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCtt                                                                                             <  7:87946/94‑1 (MQ=255)
                                                                                                              ccggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCCGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                     <  6:120217/149‑1 (MQ=255)
                                                                                                               cggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGt                                     >  7:287354/1‑148 (MQ=255)
                                                                                                                ggcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTaa                                   <  4:142620/149‑1 (MQ=255)
                                                                                                                 gcACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTaaa                                  >  4:139882/1‑149 (MQ=255)
                                                                                                                   acaGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACg                                >  4:6124/1‑149 (MQ=255)
                                                                                                                     aGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTc                              >  4:26750/1‑149 (MQ=255)
                                                                                                                       gCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCgg                            >  8:69731/1‑149 (MQ=255)
                                                                                                                              taaacaaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCAt                     >  6:27935/3‑149 (MQ=255)
                                                                                                                               aaacaaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCAt                     >  5:274322/2‑149 (MQ=255)
                                                                                                                               aaacaaAGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCAt                     >  8:158615/2‑149 (MQ=255)
                                                                                                                                aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATg                    <  5:210410/149‑1 (MQ=255)
                                                                                                                                aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATg                    >  4:81306/1‑149 (MQ=255)
                                                                                                                                aacaaAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCAGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATg                    <  8:305257/149‑1 (MQ=255)
                                                                                                                                aacaaAGCCCCGGCAGCGATGCCGGTGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATg                    <  6:274322/149‑1 (MQ=255)
                                                                                                                                   aaaGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCt                                                          >  6:252153/1‑108 (MQ=255)
                                                                                                                                   aaaGCCCCGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCt                                                          <  5:252153/108‑1 (MQ=255)
                                                                                                                                      gCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                            >  1:33343/1‑103 (MQ=255)
                                                                                                                                      gCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGAt                                                            <  2:33343/103‑1 (MQ=255)
                                                                                                                                       ccccGGCAGCGATGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGa             <  7:158615/149‑1 (MQ=255)
                                                                                                                                        cccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGAt            <  4:304923/149‑1 (MQ=255)
                                                                                                                                         ccGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATg           >  2:357373/1‑149 (MQ=255)
                                                                                                                                          cGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt          <  7:333193/149‑1 (MQ=255)
                                                                                                                                          cGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt          >  8:255087/1‑149 (MQ=255)
                                                                                                                                          cGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGt          <  3:26750/149‑1 (MQ=255)
                                                                                                                                                 gCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTc     <  7:255087/147‑1 (MQ=255)
                                                                                                                                                  cTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTcgcg  >  3:316580/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
CGCTTTGCCCGGCAGCGGTGTCAGCAGGTTTTCGAACTGTGCAAGGACCAGCAGAAGCTGGAAGGCACGGCTGTGCACAGGTTCGTGGACTTGTTTGTGATATGAAGCGGCCGGCACAGGCCCAAACAACAAAGCCCCGGCATCGCTGCCGGGGCTTTGTTTACAACCTGTAACTCAGCTTACGCCTGAACCACCGGGATCTTGGCGTTGGCTGCCGCTTCGCGGAACTCGGCGATCTGGTCGAAGCTCAGGTAGCGGTAAACGTCGGCAGCCATGCTGTCGATGTCTTTCGCG  >  NC_002947/2669529‑2669822

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: