Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,917,383 G→A 7.4% R17C (CGC→TGC)  ccmA ← cytochrome c biogenesis ABC transporter ATP‑binding protein CcmA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,917,3830GA 7.4% 196.4 / 4.1 67R17C (CGC→TGC) ccmAcytochrome c biogenesis ABC transporter ATP‑binding protein CcmA
Reads supporting (aligned to +/- strand):  ref base G (36/26);  new base A (2/3);  total (38/29)
Fisher's exact test for biased strand distribution p-value = 6.45e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

CCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTGGCAACCA  >  NC_002947/4917238‑4917531
                                                                                                                                                  |                                                                                                                                                      
ccACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGgcgc                                                                                                                                                       >  8:24303/1‑147 (MQ=255)
      agcagAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTcgcg                                                                                                                                               <  8:120343/149‑1 (MQ=255)
            aTCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCg                                                                                                                                           >  4:220626/1‑147 (MQ=255)
            aTCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCg                                                                                                                                           >  4:220628/1‑147 (MQ=255)
                            gCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCt                                                                                                                       >  7:7880/1‑149 (MQ=255)
                              tGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTg                                                                                                                     >  5:18278/1‑149 (MQ=255)
                               gCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTgg                                                                                                                    >  5:215425/1‑149 (MQ=255)
                                 aTCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCtt                                                                                                                      >  5:120438/1‑145 (MQ=255)
                                 aTCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCtt                                                                                                                      <  6:120438/145‑1 (MQ=255)
                                 aTCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAg                                                                                                                  <  6:215425/149‑1 (MQ=255)
                                  tCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAgg                                                                                                                 >  7:154299/1‑149 (MQ=255)
                                      cccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTg                                                                                                               >  2:90549/1‑147 (MQ=255)
                                       ccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAg                                                                                                            <  3:220628/149‑1 (MQ=255)
                                       ccGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAg                                                                                                            <  3:220626/149‑1 (MQ=255)
                                            cagcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGa                                                                                                                   >  4:125174/1‑137 (MQ=255)
                                            cagcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGa                                                                                                                   <  3:125174/137‑1 (MQ=255)
                                             agcagGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCAc                                                                                                      >  8:233183/1‑149 (MQ=255)
                                                  gcgcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGttt                                                                                                 <  7:233183/149‑1 (MQ=255)
                                                   cgcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTc                                                                                                <  7:284441/149‑1 (MQ=255)
                                                    gcAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCg                                                                                               >  7:109255/1‑149 (MQ=255)
                                                     cAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                                              <  3:173795/149‑1 (MQ=255)
                                                     cAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                                              >  1:7640/1‑149 (MQ=255)
                                                         aGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTa                                                                                          >  5:70675/1‑149 (MQ=255)
                                                             tGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCACCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCt                                                                                       >  1:258278/1‑149 (MQ=255)
                                                               gTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                                              <  4:20829/139‑1 (MQ=255)
                                                               gTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGc                                                                                              >  3:20829/1‑139 (MQ=255)
                                                                tCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAg                                                                                   >  7:230063/1‑149 (MQ=255)
                                                                  ttGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTa                                                                                 <  7:24303/149‑1 (MQ=255)
                                                                    gCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAg                                                                                   <  1:90549/146‑1 (MQ=255)
                                                                        cTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCACCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCtt                                                                           <  3:285691/149‑1 (MQ=255)
                                                                         tGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCttt                                                                          >  6:231198/1‑149 (MQ=255)
                                                                           ccGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCttt                                                                          <  4:20895/147‑1 (MQ=255)
                                                                            cGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTg                                                                         >  7:6352/1‑147 (MQ=255)
                                                                                   gCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                                                       >  6:341237/1‑142 (MQ=255)
                                                                                   gCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGaa                                                                       <  5:341237/142‑1 (MQ=255)
                                                                                    ccACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCaa                                                               <  7:177688/149‑1 (MQ=255)
                                                                                    ccACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCaa                                                               <  2:167161/149‑1 (MQ=255)
                                                                                       cTGATCTGCAGCGTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAacac                                                            >  3:282407/1‑149 (MQ=255)
                                                                                       cTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAacac                                                            >  3:151994/1‑149 (MQ=255)
                                                                                          aTCTGCAGCGTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCg                                                          <  4:282407/148‑1 (MQ=255)
                                                                                               cagcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCa                                                     <  2:258278/148‑1 (MQ=255)
                                                                                               cagcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAg                                                    <  6:70675/149‑1 (MQ=255)
                                                                                                 gcaTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCGTGGAGGTGAAGGGTCACGGATTCGCCGTATCTTCAGTAGCGCTTTGCAGCAGGCAGGCCCGCTCCAgca                                                  >  8:249025/1‑147 (MQ=255)
                                                                                                  caTGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAAc                                                 <  5:231198/149‑1 (MQ=255)
                                                                                                      tCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCGGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCtt                                             >  6:13514/1‑149 (MQ=255)
                                                                                                      tCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCtt                                             <  6:18278/149‑1 (MQ=255)
                                                                                                      tCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCtt                                             >  6:54294/1‑149 (MQ=255)
                                                                                                       cGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTtg                                            >  1:129772/1‑149 (MQ=255)
                                                                                                             cgcCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGtctc                                      <  8:6352/149‑1 (MQ=255)
                                                                                                                  cAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCGGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGt                                  <  5:13514/149‑1 (MQ=255)
                                                                                                                     cTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGc                               >  3:129275/1‑149 (MQ=255)
                                                                                                                      tCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcc                              >  3:275299/1‑149 (MQ=255)
                                                                                                                      tCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcc                              >  3:6440/1‑149 (MQ=255)
                                                                                                                        gAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGcccc                            <  8:7880/149‑1 (MQ=255)
                                                                                                                          aaTCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCcgc                          >  8:27142/1‑149 (MQ=255)
                                                                                                                          aaTCGAGCTGGTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGGCACGGGTTCGCCTTATCGTCAGTCGGGCTTTGAAGCAGCCAACACCGCTCCAGCCCCCTTGGGTCTCCAGTCGCCCcgc                          >  4:106098/1‑149 (MQ=255)
                                                                                                                            tCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCcgcgc                        >  1:207315/1‑149 (MQ=255)
                                                                                                                              gAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGtc                                        >  6:95416/1‑131 (MQ=255)
                                                                                                                              gAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGtc                                        <  5:95416/131‑1 (MQ=255)
                                                                                                                                  tGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGcc                  >  6:311795/1‑149 (MQ=255)
                                                                                                                                       gAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTta              >  2:37259/1‑148 (MQ=255)
                                                                                                                                       gAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTta              >  8:2482/1‑148 (MQ=255)
                                                                                                                                       gAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTta              >  1:132059/1‑148 (MQ=255)
                                                                                                                                          cagcagGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGTCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCAAGACCCCTTGTGTCTCACGTCGCCCCCCGCGCTGCCGttttcct          >  6:250173/1‑144 (MQ=255)
                                                                                                                                           agcagGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTg         <  7:2482/149‑1 (MQ=255)
                                                                                                                                            gcagGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTgg        <  4:151994/149‑1 (MQ=255)
                                                                                                                                             cagGCGCCAGTCGCGCTCGCAGGCCAGGCCCGGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTgg        >  1:215760/1‑149 (MQ=255)
                                                                                                                                               ggCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTGGCaa     <  2:207315/149‑1 (MQ=255)
                                                                                                                                                  gcCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTGGCAACca  >  8:100685/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                      
CCACCCAGCAGAATCTGCCCGGCCGTCGGCTGCATCAGCCCGGCCAGCAGGCGCAACAGGCTGGTCTTGCCGCTGCCGTTGGGGCCACTGATCTGCAGCATGTCGCCGGCGCCCAGCTCGAAATCGAGCTGTTCGAACAGCAGGCGCCAGTCGCGCTCGCAGGCCAGGCCCGCGGCTTGGAGGTGAAGGGTCACGGTTTCGCCTTATCTTCAGTAGGGCTTTGAAGCAGCCAACACCGCTCCAGAACCCTTGTGTCTCAAGTCGCCCCGCGCGCTGCCGTTATACTGGCAACCA  >  NC_002947/4917238‑4917531

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: