Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,326,404 G→A 9.4% G386E (GGG→GAG) ‡ mqo‑III → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,326,4040GA 9.4% 174.0 / 7.7 64G386E (GGG→GAG) ‡mqo‑IIImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base G (36/22);  new base A (3/3);  total (39/25)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.54e-01

TGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTC  >  NC_002947/3326260‑3326552
                                                                                                                                                |                                                                                                                                                    
tGTTTGGGCCGTTTGCCACCTTCTCGACAAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAg                                                                                                                                                  >  4:205983/1‑149 (MQ=255)
      ggCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                                            >  6:235229/1‑149 (MQ=255)
      ggCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgct                                                                                                                                             <  3:46285/148‑1 (MQ=255)
      ggCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGCACCTGCCCGATCAGGtactg                                                                                                                                            <  2:134475/149‑1 (MQ=255)
       gCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGCACCTGCCCGATCAGGtactgc                                                                                                                                           <  5:194281/149‑1 (MQ=255)
              gCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                                            >  4:351476/1‑141 (MQ=255)
              gCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                                            <  3:351476/141‑1 (MQ=255)
               ccACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGAc                                                                                                                                   >  2:12405/1‑149 (MQ=255)
                 aCCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                                            <  1:21815/138‑1 (MQ=255)
                 aCCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                                            >  2:21815/1‑138 (MQ=255)
                 aCCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGa                                                                                                                                 >  3:151215/1‑149 (MQ=255)
                         gACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCt                                                                                                                         >  1:341891/1‑149 (MQ=255)
                         gACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCt                                                                                                                         >  6:236151/1‑149 (MQ=255)
                            cAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTc                                                                                                                       >  4:73954/1‑148 (MQ=255)
                                  ccTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAgcg                                                                                                                 <  3:205983/148‑1 (MQ=255)
                                   cTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCt                                                                                                               >  4:146816/1‑149 (MQ=255)
                                       aGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTgcg                                                                                                            >  1:164934/1‑148 (MQ=255)
                                           gggCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCAcc                                                                                                        >  4:180267/1‑148 (MQ=255)
                                                        gACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAg                                                                                            <  2:164934/147‑1 (MQ=255)
                                                         aCCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGAt                                                                                                                                <  6:272120/110‑1 (MQ=255)
                                                         aCCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGAt                                                                                                                                >  5:272120/1‑110 (MQ=255)
                                                              cTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCgaaga                                                                                     >  5:3961/1‑148 (MQ=255)
                                                                 gCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGaagaaa                                                                                   <  6:325702/147‑1 (MQ=255)
                                                                  ccAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGaagatag                                                                                  >  8:24696/1‑147 (MQ=255)
                                                                  ccAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAagaaagaa                                                                                <  1:12405/149‑1 (MQ=255)
                                                                   cAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAaagaag                                                                               <  2:341891/149‑1 (MQ=255)
                                                                         cAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggc                                                                         <  6:3961/149‑1 (MQ=255)
                                                                           gCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggcgg                                                                       >  5:368463/1‑149 (MQ=255)
                                                                           gCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggcgg                                                                       >  1:310960/1‑149 (MQ=255)
                                                                             gTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCt                                                                     >  5:359076/1‑149 (MQ=255)
                                                                             gTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCt                                                                     <  5:235229/149‑1 (MQ=255)
                                                                                cacaATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCtgtg                                                                  >  6:125958/1‑149 (MQ=255)
                                                                                 acaATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTgg                                                                 >  6:293297/1‑149 (MQ=255)
                                                                                    aTGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACATGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAg                                                              >  7:212835/1‑149 (MQ=255)
                                                                                     tGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAgg                                                             <  7:330500/149‑1 (MQ=255)
                                                                                        cTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAgg                                                             >  4:128175/1‑146 (MQ=255)
                                                                                        cTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAgg                                                             <  3:128175/146‑1 (MQ=255)
                                                                                        cTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAAGCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTGCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCgg                                                          >  2:315963/1‑149 (MQ=255)
                                                                                         tGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCggg                                                         >  3:147621/1‑149 (MQ=255)
                                                                                          ggCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGc                                                        >  8:268311/1‑149 (MQ=255)
                                                                                           gCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGcc                                                       <  4:151215/149‑1 (MQ=255)
                                                                                               aTGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCg                                                   >  2:118035/1‑149 (MQ=255)
                                                                                                tGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCgg                                                  <  5:236151/149‑1 (MQ=255)
                                                                                                     cGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGc                                              <  6:359076/148‑1 (MQ=255)
                                                                                                           ggggTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAatcat                                       <  4:147621/149‑1 (MQ=255)
                                                                                                             ggTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAatcatc                                      >  8:44224/1‑148 (MQ=255)
                                                                                                             ggTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAatcatc                                      >  3:344780/1‑148 (MQ=255)
                                                                                                               tACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCaag                                   <  3:73954/149‑1 (MQ=255)
                                                                                                                 cgcgAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCaagaa                                 <  1:207810/149‑1 (MQ=255)
                                                                                                                          gACCTGGTGCAGTACGTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGc                        <  8:212835/149‑1 (MQ=255)
                                                                                                                          gACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGc                        >  3:77270/1‑149 (MQ=255)
                                                                                                                           tccgattGCAGCACCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCg                       >  5:56516/7‑149 (MQ=255)
                                                                                                                             cTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCggt                     >  2:139684/1‑147 (MQ=255)
                                                                                                                                ttGCAGCACCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCaagaaa                                >  4:104839/2‑135 (MQ=21)
                                                                                                                                ttGCAGCACCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCaagaaa                                <  3:104839/134‑1 (MQ=255)
                                                                                                                                gTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTgg                   <  3:146816/148‑1 (MQ=255)
                                                                                                                                  gcagcaCCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGGTGggcg                >  5:237779/1‑149 (MQ=25)
                                                                                                                                  gCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCgg                                                  <  6:5173/115‑1 (MQ=255)
                                                                                                                                  gCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCgg                                                  >  5:5173/1‑115 (MQ=255)
                                                                                                                                   cAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGggcgg               >  3:98232/1‑149 (MQ=255)
                                                                                                                                        ccTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGc          >  8:274226/1‑149 (MQ=255)
                                                                                                                                         cTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCt         >  6:121555/1‑149 (MQ=255)
                                                                                                                                           gATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTc        >  7:352463/1‑148 (MQ=255)
                                                                                                                                              cGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGc    >  1:26264/1‑149 (MQ=255)
                                                                                                                                               gggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTg           <  2:356598/141‑1 (MQ=255)
                                                                                                                                               gggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTg           >  1:356598/1‑141 (MQ=255)
                                                                                                                                                ggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGCGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTc  <  2:338672/149‑1 (MQ=255)
                                                                                                                                                ggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTc  <  7:24696/149‑1 (MQ=255)
                                                                                                                                                ggCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTc  <  5:293297/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTCAAGCTC  >  NC_002947/3326260‑3326552

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: ‑