Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,326,394 T→C 9.5% Y383H (TAC→CAC)  mqo‑III → malate:quinone oxidoreductase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,326,3940TC 9.5% 173.9 / 7.5 63Y383H (TAC→CAC) mqo‑IIImalate:quinone oxidoreductase
Reads supporting (aligned to +/- strand):  ref base T (33/24);  new base C (3/3);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

AAGCGCCTGATCCTGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGG  >  NC_002947/3326247‑3326539
                                                                                                                                                   |                                                                                                                                                 
aaGCGCCTGATCCTGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTa                                                                                                                                                  >  1:338672/1‑149 (MQ=255)
 aGCGCCTGATCCTGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTAc                                                                                                                                                 <  7:282850/149‑1 (MQ=255)
  gcgcCTGATCCTGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTTCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTAcc                                                                                                                                                <  5:309692/149‑1 (MQ=255)
      cTGATCCTGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGCACCTgcc                                                                                                                                            >  6:194281/1‑147 (MQ=255)
             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                                                                                                                               >  2:21815/1‑138 (MQ=255)
                              aCCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGtgctg                                                                                                                               <  1:21815/138‑1 (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                                                                                                                   >  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                   >  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)
                                                                                                                                                   |                                                                                                                                                 
AAGCGCCTGATCCTGTTTGGGCCGTTTGCCACCTTCTCGACCAAGTTCCTCAAGCAGGGCTCATTGCTTGACCTGCTCGCCAGTACCAGCGTGCACAATGCCTGGCCGATGGTGCGGGTAGGGGTACGCGAGTTCGACCTGGTGCAGTACCTGATCGGGCAGGTGCTGCAATCGGACGATGACCGCTTCGAAGCGCTGCGCACCTACTTCCCCGAGGCGAAGAAAGAAGACTGGCGGCTGTGGCAGGCGGGCCAGCGGGTGCAAATCATCAAGAAAGACGAAGCGTTGGGCGG  >  NC_002947/3326247‑3326539

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