Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,578,141 C→G 9.5% S657T (AGC→ACC)  glgB ← 1,4‑alpha‑glucan branching enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,578,1410CG 9.5% 196.1 / 9.8 63S657T (AGC→ACC) glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base C (28/29);  new base G (3/3);  total (31/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.70e-01

CCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGGCGCT  >  NC_002947/4577995‑4578282
                                                                                                                                                  |                                                                                                                                             
ccccGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCCCCACCAGCACCGGCTCGCCgctg                                                                                                                                             >  8:218203/1‑149 (MQ=255)
 cccGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctgc                                                                                                                                            >  5:92783/1‑149 (MQ=255)
   cGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctgctg                                                                                                                                          >  2:112482/1‑149 (MQ=255)
            cGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgcg                                                                                                                                  >  1:146282/1‑148 (MQ=255)
            cGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgcg                                                                                                                                  >  4:134805/1‑148 (MQ=255)
                   cGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAg                                                                                                                           <  6:92783/148‑1 (MQ=255)
                   cGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAg                                                                                                                           <  2:262554/148‑1 (MQ=255)
                        cGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTac                                                                                                                     <  5:46386/149‑1 (MQ=255)
                                  cTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCt                                                                                                           <  8:107731/149‑1 (MQ=255)
                                    ctacATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTgc                                                                                                         <  1:112482/146‑1 (MQ=255)
                                    cGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTgc                                                                                                         <  7:150344/149‑1 (MQ=255)
                                     ggCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTgcg                                                                                                        >  8:233985/1‑149 (MQ=255)
                                     ggCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTgcg                                                                                                        >  7:35231/1‑149 (MQ=255)
                                     ggCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTgc                                                                                                         <  3:134805/148‑1 (MQ=255)
                                       cATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGcgcg                                                                                                      >  2:58001/1‑149 (MQ=255)
                                        aTCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGcgt                                                                                                     >  3:302123/1‑149 (MQ=255)
                                         tCATTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACTAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGcgtc                                                                                                    >  6:239260/1‑149 (MQ=255)
                                         tCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGcgtc                                                                                                    <  8:41678/149‑1 (MQ=255)
                                           aCTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGcgtcgt                                                                                                  >  3:87494/1‑149 (MQ=255)
                                            cTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCgtcgtc                                                                                                 <  6:42531/149‑1 (MQ=255)
                                              gTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctgc                                                                                                                                            >  7:18431/1‑104 (MQ=255)
                                              gTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctgc                                                                                                                                            <  8:18431/104‑1 (MQ=255)
                                               ttCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCa                                                                                                <  4:302123/147‑1 (MQ=255)
                                               ttCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCAc                                                                                               <  3:131238/148‑1 (MQ=255)
                                                      gCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAAc                                                                                       <  8:35231/149‑1 (MQ=255)
                                                          cTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACt                                                                                   >  8:117385/1‑149 (MQ=255)
                                                               aaCGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGaaa                                                                              >  1:163494/1‑149 (MQ=255)
                                                               aaCGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGaaa                                                                              >  4:209046/1‑149 (MQ=255)
                                                                  gCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAAccc                                                                           >  2:56548/1‑149 (MQ=255)
                                                                      gCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCg                                                                       >  5:69237/1‑149 (MQ=255)
                                                                              gCATGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTg                                                                   <  3:65092/145‑1 (MQ=255)
                                                                              gCATGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTg                                                                   >  4:65092/1‑145 (MQ=255)
                                                                               cACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAAt                                                              >  7:242729/1‑149 (MQ=255)
                                                                                  gCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTg                                                                   <  1:104681/141‑1 (MQ=255)
                                                                                  gCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTg                                                                   >  2:104681/1‑141 (MQ=255)
                                                                                     aaTGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCt                                                        >  3:249712/1‑149 (MQ=255)
                                                                                      aTGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTc                                                       >  7:192582/1‑149 (MQ=255)
                                                                                          ggTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGc                                                   <  7:233985/149‑1 (MQ=255)
                                                                                               cccTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGc                                              >  8:285446/1‑149 (MQ=255)
                                                                                                ccTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCgtcgtc                                                                                                 <  3:166799/97‑1 (MQ=255)
                                                                                                ccTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCgtcgtc                                                                                                 >  4:166799/1‑97 (MQ=255)
                                                                                                    aCGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCgg                                           <  6:69237/147‑1 (MQ=255)
                                                                                                     cGGGGTACCGGTGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGctc                                        <  8:192582/149‑1 (MQ=255)
                                                                                                        ggTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGc                                                         <  6:71071/129‑1 (MQ=255)
                                                                                                        ggTACCGGGGTGAAGTTGGCCACCACCAGCACAGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGc                                                         >  5:71071/1‑129 (MQ=255)
                                                                                                              ggggTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTa                               >  6:8590/1‑149 (MQ=255)
                                                                                                                 gTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAg                            <  2:163494/149‑1 (MQ=255)
                                                                                                                  tGAAGTTGGCCACCACTAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAgg                           <  5:239260/149‑1 (MQ=255)
                                                                                                                   gAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGc                          >  5:283122/1‑149 (MQ=255)
                                                                                                                    aaGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCg                         <  8:242729/149‑1 (MQ=255)
                                                                                                                         ggCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTg                                                    >  1:288207/1‑117 (MQ=255)
                                                                                                                         ggCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTg                                                    <  2:288207/117‑1 (MQ=255)
                                                                                                                          gCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGa                    <  6:283122/148‑1 (MQ=255)
                                                                                                                           ccaccaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAgg                  <  5:8590/149‑1 (MQ=255)
                                                                                                                           ccaccaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAgg                  <  4:249712/149‑1 (MQ=255)
                                                                                                                            caccaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGt                 >  1:3006/1‑149 (MQ=255)
                                                                                                                            caccaccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGt                 >  4:235690/1‑149 (MQ=255)
                                                                                                                                 ccaGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGccc            >  3:36328/1‑149 (MQ=255)
                                                                                                                                  cagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg                   >  3:226826/1‑141 (MQ=255)
                                                                                                                                  cagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg                   <  4:226826/141‑1 (MQ=255)
                                                                                                                                  cagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg                   <  5:138340/141‑1 (MQ=255)
                                                                                                                                  cagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg                   >  6:138340/1‑141 (MQ=255)
                                                                                                                                    gCAGCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGtt                      >  5:224527/1‑136 (MQ=255)
                                                                                                                                    gCAGCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGtt                      <  6:224527/136‑1 (MQ=255)
                                                                                                                                    gCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCccac          >  1:233839/1‑148 (MQ=255)
                                                                                                                                     cACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCccacca        <  2:3006/149‑1 (MQ=255)
                                                                                                                                         cgCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGTTCGCCCACCAGgcg    <  1:58001/148‑1 (MQ=255)
                                                                                                                                           cTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCGCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGGCGCt  <  1:248026/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGGCGCT  >  NC_002947/4577995‑4578282

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: