Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,583,280 A→C 6.2% Y50D (TAT→GAT)  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,2800AC 6.2% 198.4 / 4.6 65Y50D (TAT→GAT) treSBbifunctional trehalose synthase B/maltokinase
Reads supporting (aligned to +/- strand):  ref base A (26/35);  new base C (2/2);  total (28/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 = 7.83e-01

CCTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCATG  >  NC_002947/4583133‑4583428
                                                                                                                                                   |                                                                                                                                                    
ccTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAg                                                                                                                                                     >  1:30245/1‑149 (MQ=255)
    ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTc                                                                                                                                                   >  4:40541/1‑147 (MQ=255)
    ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTc                                                                                                                                                   <  4:324066/147‑1 (MQ=255)
    ggCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTc                                                                                                                                                   <  3:40541/147‑1 (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                                                                                                       >  1:30659/1‑149 (MQ=255)
                                              gTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGt                                                                                                       <  4:125474/149‑1 (MQ=255)
                                               tacACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGGCACGGCGTGGCGAGGGATAGAACGCCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGtt                                                                                                      <  8:271170/147‑1 (MQ=255)
                                               tGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGtt                                                                                                      <  7:338136/149‑1 (MQ=255)
                                               tGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGtt                                                                                                      <  2:168562/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                                                                                         <  2:30659/148‑1 (MQ=255)
                                                             ttCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAAc                                                                                         <  7:66845/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:193757/1‑134 (MQ=255)
                                                                                   tcACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGAt                                                                                 >  1:166723/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                             >  3:250429/1‑134 (MQ=255)
                                                                                                                                       gCTCGACTATATCGCCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAg                             <  4:250429/134‑1 (MQ=255)
                                                                                                                                         tCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACg            >  2:297114/1‑149 (MQ=255)
                                                                                                                                            gCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCtt         <  5:357385/149‑1 (MQ=255)
                                                                                                                                                tataGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGcc     <  1:297114/149‑1 (MQ=255)
                                                                                                                                                 ataGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCa    >  4:37154/1‑149 (MQ=255)
                                                                                                                                                   aGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTAGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCAtg  >  7:140640/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CCTCGGCGATGAAACGGCGGGCGTCGGCCATGCTGCCGTAGTCAGGGTGCACGGCCTTGTATTCGGCGATGTCGTAGCCGTCGTCACGGCGTGGCGAGGGATAGAACGGCAGCAGCCACAAGGTGTTGACACCCAGCTCGGCGATATAGTCGAGCTTGCTGATCAGCCCGGCGAAATCGCCGATGCCGTCATTGTTCGAGTCGAAGAACGACTTGATGTGCAGCTGGTAGATGACGGCGTCCTTGTACCACAGCGGGTCGTCGATGAAGGCTGCCGGGCGGGAACGCTTGGCCATG  >  NC_002947/4583133‑4583428

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