Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,578,135 T→G 8.5% E659A (GAG→GCG) ‡ 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,1350TG 8.5% 211.0 / 10.0 71E659A (GAG→GCG) ‡glgB1,4‑alpha‑glucan branching enzyme
Reads supporting (aligned to +/- strand):  ref base T (34/31);  new base G (3/3);  total (37/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.47e-01

GCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGGCGCT  >  NC_002947/4577988‑4578282
                                                                                                                                                   |                                                                                                                                                   
gCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTc                                                                                                                                                    >  7:107731/1‑149 (MQ=255)
       ccccGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCCCCACCAGCACCGGCTCGCCgctg                                                                                                                                             >  8:218203/1‑149 (MQ=255)
        cccGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctgc                                                                                                                                            >  5:92783/1‑149 (MQ=255)
          cGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCgctgctg                                                                                                                                          >  2:112482/1‑149 (MQ=255)
                   cGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgcg                                                                                                                                  >  4:134805/1‑148 (MQ=255)
                   cGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAgcg                                                                                                                                  >  1:146282/1‑148 (MQ=255)
                          cGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAg                                                                                                                           <  2:262554/148‑1 (MQ=255)
                          cGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAg                                                                                                                           <  6:92783/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                                                                                                        >  7:35231/1‑149 (MQ=255)
                                            ggCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTgcg                                                                                                        >  8:233985/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                                                                              >  4:209046/1‑149 (MQ=255)
                                                                      aaCGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGaaa                                                                              >  1:163494/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)
                                                                                            aaTGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTcgc                                                                                                                                                  >  8:164447/1‑59 (MQ=255)
                                                                                            aaTGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTcgc                                                                                                                                                  <  7:164447/59‑1 (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                                                                                                 >  4:166799/1‑97 (MQ=255)
                                                                                                       ccTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCgtcgtc                                                                                                 <  3:166799/97‑1 (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                                                    <  2:288207/117‑1 (MQ=255)
                                                                                                                                ggCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTg                                                    >  1:288207/1‑117 (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                   >  6:138340/1‑141 (MQ=255)
                                                                                                                                         cagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg                   <  4:226826/141‑1 (MQ=255)
                                                                                                                                         cagcagCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAg                   <  5:138340/141‑1 (MQ=255)
                                                                                                                                           gCAGCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGtt                      <  6:224527/136‑1 (MQ=255)
                                                                                                                                           gCAGCGGCGAGCCGGTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGtt                      >  5:224527/1‑136 (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)
                                                                                                                                                   |                                                                                                                                                   
GCCACCGCCCCGAGATTGCCGACATTCGACCCGGCATACAGCTCGGCATCACTGTTCAGAAGCTCCTGCCAACGCTCGCCGAACGGCACGCCAATGCGGTAACCCTCACGGGGTACCGGGGTGAAGTTGGCCACCACCAGCACCGGCTCGCCGCTGCTGCTCCAGCGCAGCCAGGCGTACACGCTGTTCTGCGCGTCGTCACCGATCAACCACTGGAAACCCTGCGGCTGGCAATCCTGCTCGTGCAACGCCGGCTCTTCGCGGTACAGGCGGTTGAGGTCGCCCACCAGGCGCT  >  NC_002947/4577988‑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: