Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,583,273 G→A 6.0% A52V (GCC→GTC)  treSB ← bifunctional trehalose synthase B/maltokinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,583,2730GA 6.0% 195.0 / 2.1 66A52V (GCC→GTC) treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base G (27/35);  new base A (2/2);  total (29/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.40e-01

GTGGGCCTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTT  >  NC_002947/4583128‑4583421
                                                                                                                                                 |                                                                                                                                                    
gtgGGCCTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGc                                                                                                                                                     <  5:83068/147‑1 (MQ=255)
gtgGGCCTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGa                                                                                                                                                   >  3:73743/1‑149 (MQ=255)
     ccTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAg                                                                                                                                              >  1:30245/1‑149 (MQ=255)
         ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTc                                                                                                                                            <  3:40541/147‑1 (MQ=255)
         ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTc                                                                                                                                            <  4:324066/147‑1 (MQ=255)
         ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTc                                                                                                                                            >  4:40541/1‑147 (MQ=255)
         ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGa                                                                                                                                          <  4:341907/149‑1 (MQ=255)
              tGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTg                                                                                                                                     <  6:55204/149‑1 (MQ=255)
                 aaCGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGc                                                                                                                                    <  6:99837/147‑1 (MQ=255)
                 aaCGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGc                                                                                                                                    >  5:99837/1‑147 (MQ=255)
                         gcgTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGTGATATAGTCGAGCTTGCTGATCAGccc                                                                                                                          >  8:266417/1‑149 (MQ=255)
                         gcgTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGAAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGccc                                                                                                                          >  1:255887/1‑149 (MQ=255)
                           gTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCgg                                                                                                                        <  6:177074/149‑1 (MQ=255)
                            tCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGc                                                                                                                       <  7:131404/149‑1 (MQ=255)
                              ggCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCg                                                                                                                      >  7:168550/1‑148 (MQ=255)
                               gCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGa                                                                                                                     >  1:86374/1‑148 (MQ=255)
                               gCCATGCTACCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGCCGa                                                                                                                     >  4:33164/1‑148 (MQ=255)
                                 cATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAAt                                                                                                                  >  6:163385/1‑149 (MQ=255)
                                 cATGCTGCCGTAGTCAGGGTGCACGGCCTTCTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGCGCTTGCTGATCAGCCCCGCGAAAt                                                                                                                  >  5:51901/1‑149 (MQ=255)
                                       gcCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCg                                                                                                             >  8:338136/1‑148 (MQ=255)
                                        ccGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGTCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCg                                                                                                             <  4:165760/147‑1 (MQ=255)
                                        ccGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCg                                                                                                             <  4:160108/147‑1 (MQ=255)
                                          gTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGAACCGCCCGGCGAAATCGCCGATGc                                                                                                         >  7:233223/1‑149 (MQ=255)
                                             gTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGt                                                                                                      >  3:125474/1‑149 (MQ=255)
                                                   gTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGt                                                                                                <  4:125474/149‑1 (MQ=255)
                                                   gTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGt                                                                                                >  1:30659/1‑149 (MQ=255)
                                                    tacACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGGCACGGCGTGGCGAGGGATAGAACGCCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGtt                                                                                               <  8:271170/147‑1 (MQ=255)
                                                    tGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGtt                                                                                               <  2:168562/149‑1 (MQ=255)
                                                    tGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGtt                                                                                               <  7:338136/149‑1 (MQ=255)
                                                       aCGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGa                                                                                            <  8:168550/149‑1 (MQ=255)
                                                        cGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAg                                                                                           >  1:356095/1‑149 (MQ=255)
                                                             ttGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCg                                                                                        >  6:357385/1‑147 (MQ=255)
                                                                tATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCgaagaa                                                                                   >  6:201495/1‑149 (MQ=255)
                                                                  ttCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAAc                                                                                  <  7:66845/148‑1 (MQ=255)
                                                                  ttCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAAc                                                                                  <  2:30659/148‑1 (MQ=255)
                                                                       cGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGACGTCATTGTTCGAGTCGAAGAACGACtt                                                                             >  2:287910/1‑148 (MQ=255)
                                                                         aTGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAt                                                                          <  2:107149/149‑1 (MQ=255)
                                                                            tCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGTGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAtgt                                                                        <  7:266417/148‑1 (MQ=255)
                                                                                   cgtcgtCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTgg                                                                <  5:240807/149‑1 (MQ=255)
                                                                                    gtcgtcACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGt                                                               <  1:287910/149‑1 (MQ=255)
                                                                                        tcACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAt                                                                          >  1:166723/1‑134 (MQ=255)
                                                                                        tcACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAt                                                                          >  1:193757/1‑134 (MQ=255)
                                                                                        tcACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAt                                                                          <  2:166723/134‑1 (MQ=255)
                                                                                        tcACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAt                                                                          <  2:193757/134‑1 (MQ=255)
                                                                                                  gCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTcc                                                 <  5:201495/149‑1 (MQ=255)
                                                                                                  gCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTc                                                  <  2:106272/148‑1 (MQ=255)
                                                                                                     aGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTg                                              <  4:162555/149‑1 (MQ=255)
                                                                                                       gcaTAGAACGGCAGCAGCCACAAGGTGATGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTa                                            <  4:73743/147‑1 (MQ=255)
                                                                                                            gAACGGCAGAAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACcaca                                       <  2:255887/149‑1 (MQ=255)
                                                                                                             aaCGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACcaca                                       <  3:33164/148‑1 (MQ=255)
                                                                                                                     cagcCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGgtcgt                              <  1:264999/149‑1 (MQ=255)
                                                                                                                          acaAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCTGGTCGTCGATg                         <  4:330776/149‑1 (MQ=255)
                                                                                                                            aaGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATg                         >  6:342513/1‑147 (MQ=255)
                                                                                                                                tgtTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCt                   >  2:247788/1‑149 (MQ=255)
                                                                                                                                  ttGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGc                 >  6:83277/1‑149 (MQ=255)
                                                                                                                                     acacCCAGCTGGGCGATTTAGTCGAGATTGCTGTTCTGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGATTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgg              <  2:356095/149‑1 (MQ=255)
                                                                                                                                        gcaaGCTCGACTATATCGCCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcg           >  3:160178/4‑149 (MQ=255)
                                                                                                                                        ccgaGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcg           <  5:83277/146‑1 (MQ=255)
                                                                                                                                         ccAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcgg          >  6:56053/1‑149 (MQ=255)
                                                                                                                                         ccAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcgg          >  6:321041/1‑149 (MQ=255)
                                                                                                                                         caaGCTCGACTATATCGCCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcgg          <  4:160178/147‑1 (MQ=255)
                                                                                                                                          cAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcggg         <  2:86374/149‑1 (MQ=255)
                                                                                                                                          cAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCcgggcggg         <  5:163385/149‑1 (MQ=255)
                                                                                                                                            gCTCGACTATATCGCCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAg                      <  4:250429/134‑1 (MQ=255)
                                                                                                                                            gCTCGACTATATCGCCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAg                      >  3:250429/1‑134 (MQ=255)
                                                                                                                                              tCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACg     >  2:297114/1‑149 (MQ=255)
                                                                                                                                                 gCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCtt  <  5:357385/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                    
GTGGGCCTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTT  >  NC_002947/4583128‑4583421

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: