Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,984,043 A→C 10.2% intergenic (+37/+153) ilvE → / ← PP_3512 branched‑chain‑amino acid aminotransferase/membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,984,0430AC10.2% 167.2 / 9.9 59intergenic (+37/+153)ilvE/PP_3512branched‑chain‑amino acid aminotransferase/membrane protein
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base C (3/3);  total (31/28)
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.09e-01

CACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGA  >  NC_002947/3983896‑3984190
                                                                                                                                                   |                                                                                                                                                   
cACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAAt                                                                                                                                                    >  7:65848/1‑149 (MQ=255)
cACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAAt                                                                                                                                                    >  5:232696/1‑149 (MQ=255)
     ccTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCAt                                                                                                                                               <  6:232696/149‑1 (MQ=255)
         gAAAAGGTCGGCCCGGTTACCCAGCAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGGCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGccc                                                                                                                                           <  3:252721/149‑1 (MQ=255)
         gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTGCAAGGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGcc                                                                                                                                            <  7:69972/148‑1 (MQ=255)
         gAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGc                                                                                                                                             <  8:65848/147‑1 (MQ=255)
                   gCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCa                                                                                                                                 >  5:197640/1‑149 (MQ=255)
                     ccGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCAtg                                                                                                                               <  4:278787/149‑1 (MQ=255)
                            cccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTa                                                                                                                         <  8:304079/148‑1 (MQ=255)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                       >  1:106866/1‑149 (MQ=255)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                       >  5:177041/1‑149 (MQ=255)
                             ccAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTaca                                                                                                                       >  2:79225/1‑149 (MQ=255)
                                  aGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGcc                                                                                                                  >  4:254743/1‑149 (MQ=255)
                                       tACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTGTTTTCATGTGTTTAGACGgccgccgc                                                                                                             >  4:276381/1‑147 (MQ=255)
                                         cAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                            >  2:13922/1‑148 (MQ=255)
                                          aaCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                            >  8:27829/1‑147 (MQ=255)
                                          aaCGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCg                                                                                                            >  7:60169/1‑147 (MQ=255)
                                                    ccGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGtt                                                                                                >  8:2500/1‑149 (MQ=255)
                                                       gTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTAc                                                                                             >  4:279875/1‑149 (MQ=255)
                                                          tCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcg                                                                                          >  3:328618/1‑149 (MQ=255)
                                                           ccAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                         >  5:359024/1‑149 (MQ=255)
                                                             agagCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgcg                                                                                        >  3:157275/1‑148 (MQ=255)
                                                                 cggcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAg                                                                                   >  8:159966/1‑149 (MQ=255)
                                                                   gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                 <  6:177041/149‑1 (MQ=255)
                                                                   gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                 >  3:170082/1‑149 (MQ=255)
                                                                   gcgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                 <  6:359024/149‑1 (MQ=255)
                                                                    cgATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGc                                                                                 >  8:262193/1‑148 (MQ=255)
                                                                        gTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCacac                                                                            >  6:360525/1‑149 (MQ=255)
                                                                          agagGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGttt                                                                                                                                      >  2:103379/1‑89 (MQ=37)
                                                                          agagGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCca                                                                               <  2:158873/144‑1 (MQ=255)
                                                                          agagGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCca                                                                               >  1:158873/1‑144 (MQ=255)
                                                                          agagGCACCGGCAGGCTGGATCGTCAAGGTCGCATGAAACGCATCAGCTGTATGAAAACGGGGCATGCCATTGCTGGCATGCCCCGttt                                                                                                                                      <  1:103379/89‑1 (MQ=25)
                                                                                       ggCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCaataat                                                             <  6:147801/149‑1 (MQ=255)
                                                                                         cTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGc                                                           <  3:148734/149‑1 (MQ=255)
                                                                                          tGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                         >  7:158160/1‑118 (MQ=255)
                                                                                          tGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCcgc                                                                                         <  8:158160/118‑1 (MQ=255)
                                                                                          tGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGcc                                                          >  7:186626/1‑149 (MQ=255)
                                                                                            gATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTg                                                        <  2:234303/149‑1 (MQ=255)
                                                                                             aTCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTgg                                                       <  8:60169/149‑1 (MQ=255)
                                                                                              tCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACAGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTggg                                                      <  3:95544/149‑1 (MQ=255)
                                                                                                gTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGga                                                     <  7:159966/148‑1 (MQ=255)
                                                                                                       tCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCg                                             >  4:306658/1‑149 (MQ=255)
                                                                                                       tCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCg                                             >  3:320898/1‑149 (MQ=255)
                                                                                                         gCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGtt                                           <  4:328618/149‑1 (MQ=255)
                                                                                                          ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAACCAATAATGCCTGGGAGATCACCg                                             >  2:340829/1‑146 (MQ=255)
                                                                                                          ccTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAACCAATAATGCCTGGGAGATCACCGtt                                           <  1:340829/148‑1 (MQ=255)
                                                                                                             gAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGt                                         >  7:201932/1‑147 (MQ=255)
                                                                                                              aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCa                                       <  7:2500/148‑1 (MQ=255)
                                                                                                              aaaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCa                                       <  1:13922/148‑1 (MQ=255)
                                                                                                               aaCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGt                                     <  3:254743/149‑1 (MQ=255)
                                                                                                                aCGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTg                                    <  4:320898/149‑1 (MQ=255)
                                                                                                                           atgaaaaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCa                            <  1:24491/141‑1 (MQ=38)
                                                                                                                           atgaaaaaaCGGGGCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCa                            >  2:24491/6‑146 (MQ=38)
                                                                                                                              tGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCg                      <  2:106866/149‑1 (MQ=255)
                                                                                                                                 aaaCGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTACCggg                    >  8:143001/1‑148 (MQ=255)
                                                                                                                                     ggggCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGc               <  6:197640/149‑1 (MQ=255)
                                                                                                                                       ggCATGCCATTGCTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGccc             >  2:89447/1‑149 (MQ=255)
                                                                                                                                               aGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATAACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCggg     <  1:79225/149‑1 (MQ=255)
                                                                                                                                                gCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGt    <  5:334374/149‑1 (MQ=255)
                                                                                                                                                  aaTGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGa  >  6:126889/1‑149 (MQ=255)
                                                                                                                                                  aaTGGCATGCACCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGACTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGa  >  6:193816/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CACGACCTGGAAAAGGTCGGCCCGGTTACCCAGAAGCTCTACAACGAGCTGACCGGTATCCAGAGCGGCGATGTAGAGGCACCGGCAGGCTGGATCGTCAAGGTCGCCTGAAACGCATCAGCTGTATGAAAACGGGGCATGCCAGCAATGGCATGCCCCGTTTTTTCATGTGTTTACACGGCCCCCTCGTAGGAGCGGGTTTACCCGCGAATAGGCCACACCCGCCAGCCAATAATGCCTGGGAGATCACCGTTGTCAGTGCAGCCGCATTAGCGGGCAAGCCCGCGCCGGGTGA  >  NC_002947/3983896‑3984190

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: