Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,734,028 T→G 5.1% E536A (GAG→GCG)  sucA ← 2‑oxoglutarate dehydrogenase subunit E1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,734,0280TG 5.1% 251.8 / 4.2 78E536A (GAG→GCG) sucA2‑oxoglutarate dehydrogenase subunit E1
Reads supporting (aligned to +/- strand):  ref base T (39/35);  new base G (2/2);  total (41/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.19e-01

CTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTGC  >  NC_002947/4733883‑4734175
                                                                                                                                                  |                                                                                                                                                   
cTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCGTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgc                                                                                                                                                  >  3:101371/1‑149 (MQ=255)
 ttCGACACCTGACGCTGGACCACGGAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgcg                                                                                                                                                 >  8:167454/1‑149 (MQ=255)
  tCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCgg                                                                                                                                                >  7:204847/1‑149 (MQ=255)
     acacCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTc                                                                                                                                             <  2:289083/149‑1 (MQ=255)
        ccTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGc                                                                                                                                          <  6:164364/149‑1 (MQ=255)
             cGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                     <  7:88844/149‑1 (MQ=255)
              gCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                     >  1:170213/1‑148 (MQ=255)
              gCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCtt                                                                                                                                     >  7:17803/1‑148 (MQ=255)
                  gACCACGGAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                                  <  7:167454/147‑1 (MQ=255)
                    ccACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgc                                                                                                                                           <  1:103842/136‑2 (MQ=255)
                    ccACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACCGCGAGCTGTTCgc                                                                                                                                           >  2:103842/1‑135 (MQ=255)
                    ccACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAgg                                                                                                                              >  3:66210/1‑149 (MQ=255)
                     cACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcgc                                                                                                                                                  >  8:141296/1‑128 (MQ=255)
                      aCGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCt                                                                                                                            >  3:56108/1‑149 (MQ=255)
                       cGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCtt                                                                                                                           >  5:413002/1‑149 (MQ=255)
                             cTTCCGGCAGCTCGAGCAGCTTGGCCGAGCAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTAc                                                                                                                      >  1:24083/1‑149 (MQ=255)
                                ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCAc                                                                                                                   <  2:170213/148‑1 (MQ=255)
                                ccGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCAc                                                                                                                   <  8:17803/148‑1 (MQ=255)
                                 cGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGt                                                                                                                 >  1:321149/1‑149 (MQ=255)
                                       cTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcg                                                                                                                                                   >  2:160644/1‑109 (MQ=255)
                                       cTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTcg                                                                                                                                                   <  1:160644/109‑1 (MQ=255)
                                            gcagcTTGGCCGACAGATCCTGCAGGGTCTGGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTg                                                                                                      <  4:138399/149‑1 (MQ=255)
                                              agcTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                    <  4:65392/149‑1 (MQ=255)
                                              agcTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                    <  8:92842/149‑1 (MQ=255)
                                               gcTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                    <  3:95033/148‑1 (MQ=255)
                                                cTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTc                                                                                                    >  4:95033/1‑147 (MQ=255)
                                                     ccGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGcgcg                                                                                             >  5:53418/1‑149 (MQ=255)
                                                      cGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                                  <  3:178546/111‑1 (MQ=255)
                                                      cGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGAc                                                                                                                                  >  4:178546/1‑111 (MQ=255)
                                                      cGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGt                                                                                            >  5:187946/1‑149 (MQ=255)
                                                       gACAGATGATGCAGTGTCTTGTGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGACAGTCGACGAAAAGCTCGCGGTTCGGCTACTTTACCAGGCTTTTTACCAGGTGCAGGCCTTTGTCCAGAGCGtt                                                                                           <  4:101371/149‑1 (MQ=255)
                                                       gACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGtt                                                                                           <  2:24083/149‑1 (MQ=255)
                                                           ataTCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTc                                                                                                                                             <  8:362203/94‑1 (MQ=255)
                                                           aGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCg                                                                                        >  4:299882/1‑149 (MQ=255)
                                                             aTCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCgg                                                                                       >  7:131098/1‑148 (MQ=255)
                                                                cTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACt                                                                                   >  2:21281/1‑149 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                <  6:53418/149‑1 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                <  7:209099/149‑1 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                <  4:56108/149‑1 (MQ=255)
                                                                   cAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                <  6:187946/149‑1 (MQ=255)
                                                                     gggTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgtcg                                                                              <  8:131098/149‑1 (MQ=255)
                                                                       gTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGAt                                                                            <  8:204847/149‑1 (MQ=255)
                                                                        tCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTAACCACGTGCAGGCCATTGTCCAGCGCGTTGCTGTACTCGTCGAt                                                                            >  4:16195/1‑148 (MQ=255)
                                                                              ggTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGcc                                                                     >  3:310596/1‑149 (MQ=255)
                                                                              ggTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGc                                                                      <  3:299882/148‑1 (MQ=255)
                                                                              ggTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCAAGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCAATGTCCAGCGCGTTGCGGTACTCGTCGATCTTGTcc                                                                     >  3:305057/1‑149 (MQ=255)
                                                                               gTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                    >  6:199087/1‑149 (MQ=255)
                                                                                 cGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                    >  1:294840/1‑147 (MQ=255)
                                                                                    aaCGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTAc                                                                                    >  3:252391/1‑128 (MQ=255)
                                                                                    aaCGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTAc                                                                                    <  4:252391/128‑1 (MQ=255)
                                                                                           gtCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCgg                                                        <  4:66210/149‑1 (MQ=255)
                                                                                              gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCCGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGAGGTACTCGTCGATCTTGGCCTGAGAGCTCTCGGCGt                                                     >  8:294986/1‑149 (MQ=255)
                                                                                              gTGACGCGCGGTCCAGGCATGGCCCAGGTACGGAAGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTAATCGTCGATCTTGGCCTGAGCGCGCTCGGCGt                                                     >  7:72359/1‑149 (MQ=255)
                                                                                                  cgcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGAt                                                 <  8:302502/149‑1 (MQ=255)
                                                                                                   gcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACtcgt                                                                                >  4:342862/1‑117 (MQ=255)
                                                                                                   gcgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATg                                                >  1:266554/1‑149 (MQ=255)
                                                                                                    cgcgGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGc                                               <  6:413002/149‑1 (MQ=255)
                                                                                                         tCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCGATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGAGCGCTCGGCGTCGATGCGACCg                                          >  5:301773/1‑149 (MQ=255)
                                                                                                            aGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGcc                                       <  5:199087/149‑1 (MQ=255)
                                                                                                                    cccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAg                               >  1:371361/1‑149 (MQ=255)
                                                                                                                    cccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAg                               <  2:321149/149‑1 (MQ=255)
                                                                                                                     ccAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAgc                              >  6:420018/1‑149 (MQ=255)
                                                                                                                       aGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                    >  1:78092/1‑109 (MQ=255)
                                                                                                                       aGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCt                                                                    <  2:78092/109‑1 (MQ=255)
                                                                                                                             ggACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGc                                                                      <  3:15015/101‑1 (MQ=255)
                                                                                                                             ggACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGc                                                                      >  4:15015/1‑101 (MQ=255)
                                                                                                                             ggACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGc                      <  7:252609/149‑1 (MQ=255)
                                                                                                                                        gccgaACCGCGAGCTGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCt              >  1:374351/3‑146 (MQ=255)
                                                                                                                                        gccgaACCGCGAGCTGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCt              <  2:374351/144‑1 (MQ=255)
                                                                                                                                          cgaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACgg         <  2:266554/149‑1 (MQ=255)
                                                                                                                                          cgaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACgg         <  2:371361/149‑1 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCt                                      <  6:127735/119‑1 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCt                                      >  5:127735/1‑119 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATAc                  <  6:377562/139‑1 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATAc                  >  5:377562/1‑139 (MQ=255)
                                                                                                                                           gaACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTAGGCATACAGCTCACgg         >  3:48983/1‑148 (MQ=255)
                                                                                                                                              cAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAg     >  1:234036/1‑149 (MQ=255)
                                                                                                                                                 cTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTgc  >  7:167772/1‑149 (MQ=255)
                                                                                                                                                  tCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCa                                                                                                  <  4:310529/52‑1 (MQ=255)
                                                                                                                                                  tCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCa                                                                                                  >  3:310529/1‑52 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CTTCGACACCTGACGCTGGACCACGAAACCTTCCGGCAGCTCGAGCAGCTTGGCCGACAGATCCTGCAGGGTCTTGAGGTCGAAACGGGTGTCGTGACGCGCGGTCCAGGCATGGCCCAGGTACGGACGCCAGTCGACGAACAGCTCGCGGTTCGGCTCCTTGACCAGGCTTTTTACCACGTGCAGGCCATTGTCCAGCGCGTTGCGGTACTCGTCGATCTTGGCCTGAGCGCGCTCGGCGTCGATGCGACCGGCCTGGATCAGCTGCTCGGCATACAGCTCACGGGTAGTGC  >  NC_002947/4733883‑4734175

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: