Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,567,880 A→G 5.7% T517A (ACC→GCC)  treZ → malto‑oligosyltrehalose trehalohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,567,8800AG 5.7% 198.4 / 2.0 71T517A (ACC→GCC) treZmalto‑oligosyltrehalose trehalohydrolase
Reads supporting (aligned to +/- strand):  ref base A (31/35);  new base G (2/2);  total (34/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 = 4.33e-01

GCCTCTTCACGCCCGCAGAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATACCCACCTGGCT  >  NC_002947/4567733‑4568023
                                                                                                                                                   |                                                                                                                                               
gCCTCTTCACGCCCGCAGAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGAc                                                                                                                                                >  8:152253/1‑149 (MQ=255)
gCCTCTTCACGCCCGCAGAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGAc                                                                                                                                                >  2:353303/1‑149 (MQ=255)
 ccTCTTCACGCCCGCAGAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGAcc                                                                                                                                               >  2:325566/1‑149 (MQ=255)
  ctctTCACGCCCGCAGAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACcg                                                                                                                                              >  3:100014/1‑149 (MQ=255)
                 gAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTgg                                                                                                                               <  1:325566/149‑1 (MQ=255)
                    cAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggc                                                                                                                             >  7:176973/1‑148 (MQ=255)
                      aGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggcaac                                                                                                                          <  6:32821/149‑1 (MQ=255)
                        gCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggcaacgg                                                                                                                        >  6:249038/1‑149 (MQ=255)
                           gTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGgggaacggcaa                                                                                                                     <  1:315610/149‑1 (MQ=255)
                              aTCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAaca                                                                                                                   <  8:58759/148‑1 (MQ=255)
                                  ccGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGGCCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTg                                                                                                              >  5:329112/1‑149 (MQ=255)
                                  ccGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGGCCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTg                                                                                                              >  5:334595/1‑149 (MQ=255)
                                   cGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCgcgcgc                                                                                                                                         >  4:231647/1‑121 (MQ=255)
                                   cGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTgc                                                                                                             <  1:315383/149‑1 (MQ=255)
                                   cGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTgc                                                                                                             >  2:10973/1‑149 (MQ=255)
                                         gCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCAt                                                                                                         <  8:167070/147‑1 (MQ=255)
                                         gCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATc                                                                                                        <  5:249038/148‑1 (MQ=255)
                                         gCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATc                                                                                                        >  7:16154/1‑148 (MQ=255)
                                           aCGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGAc                                                                                                     <  1:10973/149‑1 (MQ=255)
                                             ggCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGGCCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCt                                                                                                   <  6:334595/149‑1 (MQ=255)
                                             ggCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGGCCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCt                                                                                                   <  6:329112/149‑1 (MQ=255)
                                               ccTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCt                                                                                                                                 >  2:70788/1‑117 (MQ=255)
                                               ccTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCt                                                                                                                                 <  1:70788/117‑1 (MQ=255)
                                                 tttACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAAc                                                                                               >  7:166939/1‑149 (MQ=255)
                                                      ccgCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGc                                                                                          <  8:294952/147‑1 (MQ=255)
                                                      cAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGc                                                                                          <  4:100014/149‑1 (MQ=255)
                                                      cAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGc                                                                                          >  4:168654/1‑149 (MQ=255)
                                                         cAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCTCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGa                                                                                       <  8:130212/149‑1 (MQ=255)
                                                                  gACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggcaacggcaac                                                                                                                    >  6:294405/1‑111 (MQ=255)
                                                                  gACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggcaacggcaac                                                                                                                    <  5:294405/111‑1 (MQ=255)
                                                                  gACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGc                                                                               <  1:353303/148‑1 (MQ=255)
                                                                   aCCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAg                                                                             <  8:16154/149‑1 (MQ=255)
                                                                    ccTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAgg                                                                            <  3:49119/149‑1 (MQ=255)
                                                                     cTGCGCCGCCGCGACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggcaacggcaa                                                                                                                     >  8:306426/1‑107 (MQ=255)
                                                                     cTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGggcaacggcaa                                                                                                                     <  7:306426/107‑1 (MQ=255)
                                                                      tGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACAct                                                                                                                 <  1:201973/110‑1 (MQ=255)
                                                                      tGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACAct                                                                                                                 >  2:201973/1‑110 (MQ=255)
                                                                       gcgcCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCgg                                                                         >  6:199761/1‑149 (MQ=255)
                                                                       gcgcCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCgg                                                                         >  6:204776/1‑149 (MQ=255)
                                                                         gccgccgccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTc                                                                       <  8:176973/149‑1 (MQ=255)
                                                                         gccgccgccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTc                                                                       <  7:152253/149‑1 (MQ=255)
                                                                          ccgccgccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCg                                                                      <  3:168654/149‑1 (MQ=255)
                                                                           tcacgccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAATCGCTGACCTCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGACCGACACCCCGCAGGCGGTCGa                                                                     <  7:54789/146‑1 (MQ=255)
                                                                               gccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCAGTCGa                                                                     >  8:54789/1‑145 (MQ=255)
                                                                                ccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCTCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAAc                                                                                               >  3:129117/1‑118 (MQ=255)
                                                                                ccACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCTCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAAc                                                                                               <  4:129117/118‑1 (MQ=255)
                                                                                    atcatcCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGACCTGCCcgc                                                            >  8:327197/1‑149 (MQ=255)
                                                                                    atcatcCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCcgc                                                            >  5:83320/1‑149 (MQ=255)
                                                                                       atcCCGCCCTTGCCCGTCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCgcgc                                                          <  6:252575/148‑1 (MQ=255)
                                                                                       atcCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCgcgc                                                          <  2:77964/148‑1 (MQ=255)
                                                                                          ccGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATc                                                      >  7:11726/1‑149 (MQ=255)
                                                                                           cGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATc                                                      >  6:42300/1‑148 (MQ=255)
                                                                                                ttGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGCCCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCCGGCGGTCGAACTGCACGCGCCATCCACCCg                                                >  8:36894/1‑149 (MQ=255)
                                                                                                  gCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCATCACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGc                                              >  1:113257/1‑149 (MQ=255)
                                                                                                  gCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACgg                                               <  5:204776/148‑1 (MQ=255)
                                                                                                  gCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACgg                                               <  2:256581/148‑1 (MQ=255)
                                                                                                  gCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACgg                                               <  5:199761/148‑1 (MQ=255)
                                                                                                       gCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTc                                         <  2:300478/149‑1 (MQ=255)
                                                                                                            cGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGAcag                                    <  5:42300/149‑1 (MQ=255)
                                                                                                                              gTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCAccc                  <  8:11726/149‑1 (MQ=255)
                                                                                                                              gTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCAccc                  <  6:83320/149‑1 (MQ=255)
                                                                                                                                 cACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGa               >  2:348574/1‑149 (MQ=255)
                                                                                                                                   cGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCATCACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATa             <  2:113257/149‑1 (MQ=255)
                                                                                                                                     cTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGTCATCCGCACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATAcc           >  4:293929/1‑149 (MQ=255)
                                                                                                                                      tGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGc                                                                                          <  6:128717/69‑1 (MQ=255)
                                                                                                                                      tGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGc                                                                                          >  5:128717/1‑69 (MQ=255)
                                                                                                                                        aCAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATACCCAc        >  5:362381/1‑149 (MQ=255)
                                                                                                                                          aaaGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCAGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATACCCACCt      >  5:227509/1‑149 (MQ=255)
                                                                                                                                            aGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGCACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGCCCCCCGCCCCCCGGATACCCACCTgg    >  4:342103/1‑149 (MQ=255)
                                                                                                                                              cgcTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATACCCACCTGGCt  >  8:141723/1‑149 (MQ=255)
                                                                                                                                              cgcTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATACCCACCTGGCt  <  1:348574/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
GCCTCTTCACGCCCGCAGAACAAGGCGGTGATCGCCGGCTGGCACGGCCTTTACCAGCAACTGCTCGACCTGCGCCGCCGCCACATCATCCCGCACTTGCCCGGCAGCCGTGCGCTGGGCGCCGAGGTGCACGCTGACAAAGCGCTGACCGCGCGCTGGCAGCTGGGCAACGGCAACACTCTGCGCATCGACCTGAACCTGGCCGACACCCCGCAGGCGGTCGAACTGCCCGCGCCATCCACACGGCTGTTCGACAGCAGCGACACCCGCCACCCGGATACCCACCTGGCT  >  NC_002947/4567733‑4568023

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑