Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,326,401 T→C 9.5% I385T (ATC→ACC) ‡ mqo‑III → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,326,4010TC 9.5% 174.8 / 7.8 63I385T (ATC→ACC) ‡mqo‑IIImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (35/22);  new base C (3/3);  total (38/25)
Fisher's exact test for biased strand distribution p-value = 6.74e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.07e-01

TGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTC  >  NC_002947/3326260‑3326546
                                                                                                                                             |                                                                                                                                                 
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                                                                                                                                      <  3:351476/141‑1 (MQ=255)
              gCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                                      >  4:351476/1‑141 (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                                                                                                                   >  6:236151/1‑149 (MQ=255)
                         gACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCt                                                                                                                   >  1:341891/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                                                                                                                          >  5:272120/1‑110 (MQ=255)
                                                         aCCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGAt                                                                                                                          <  6:272120/110‑1 (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                                                                 >  1:310960/1‑149 (MQ=255)
                                                                           gCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTggcgg                                                                 >  5:368463/1‑149 (MQ=255)
                                                                             gTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCt                                                               <  5:235229/149‑1 (MQ=255)
                                                                             gTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCt                                                               >  5:359076/1‑149 (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                          <  3:104839/134‑1 (MQ=255)
                                                                                                                                ttGCAGCACCTGCCCGATCAGGTACTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCaagaaa                          >  4:104839/2‑135 (MQ=21)
                                                                                                                                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)
                                                                                                                                             |                                                                                                                                                 
TGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGGCGTGCTC  >  NC_002947/3326260‑3326546

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