Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,892 C→G 6.3% G228R (GGC→CGC)  ldcC ← lysine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,677,8920CG 6.3% 196.9 / 2.2 61G228R (GGC→CGC) ldcClysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (30/27);  new base G (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACCCGTCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTA  >  NC_002947/4677749‑4678040
                                                                                                                                               |                                                                                                                                                    
aCCCGTCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGt                                                                                                                                                 <  3:851/149‑1 (MQ=255)
     tCATGATGCGCGCGTGCACCACCGACTTGGGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTgcgc                                                                                                                                            <  8:77394/149‑1 (MQ=255)
        tgatgaTCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGcgcc                                                                                                                                         >  6:292089/1‑149 (MQ=255)
            gatCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgat                                                                                                                                      >  4:128534/1‑146 (MQ=255)
             atCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgatc                                                                                                                                     <  3:128534/148‑4 (MQ=255)
              tCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgatca                                                                                                                                    >  8:247140/1‑144 (MQ=255)
              tCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTgc                                                                                                                                              <  1:288737/138‑1 (MQ=255)
              tCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTgc                                                                                                                                              >  2:288737/1‑138 (MQ=255)
                gcgTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCGCAACTTCGGcgccgatcaca                                                                                                                                  <  7:247140/148‑7 (MQ=255)
                gcgTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGc                                                                                                                                 <  3:148524/149‑1 (MQ=255)
                   tGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTc                                                                                                                              <  5:243496/149‑1 (MQ=255)
                     caccacCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                            >  8:218683/1‑149 (MQ=255)
                          ccGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCg                                                                                                                       <  5:120733/149‑1 (MQ=255)
                          ccGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCg                                                                                                                       <  8:46139/149‑1 (MQ=255)
                            gACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCGATTAATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                            >  2:156995/1‑142 (MQ=255)
                            gACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTAATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                            <  1:156995/142‑1 (MQ=255)
                                       gTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCtg                                                                                                          >  2:186623/1‑149 (MQ=255)
                                            ggTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTgg                                                                                                     >  2:324120/1‑149 (MQ=255)
                                                caccaggaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCgg                                                                                                                            <  8:270295/122‑1 (MQ=255)
                                                caccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGt                                                                                                    >  7:304872/1‑146 (MQ=255)
                                                caccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGt                                                                                                    <  8:304872/146‑1 (MQ=255)
                                                 accagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGag                                                                                                >  2:211052/1‑149 (MQ=255)
                                                  ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCt                                                                                                           >  1:294929/1‑137 (MQ=255)
                                                  ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCt                                                                                                           <  2:294929/137‑1 (MQ=255)
                                                  ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGag                                                                                                >  6:224639/1‑148 (MQ=255)
                                                   cagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTc                                                                                                   >  3:55150/1‑144 (MQ=255)
                                                   cagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTc                                                                                                   <  4:55150/144‑1 (MQ=255)
                                                            gtcgtcACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAAcc                                                                                       >  8:20254/1‑147 (MQ=255)
                                                                tcACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTc                                                                                 <  5:224639/149‑1 (MQ=255)
                                                                   cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGc                                                                                                             >  2:99359/1‑118 (MQ=255)
                                                                   cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGc                                                                                                             <  1:99359/118‑1 (MQ=255)
                                                                   cGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCgg                                                                              <  1:324120/149‑1 (MQ=255)
                                                                        gACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACgg                                                                         <  7:20254/149‑1 (MQ=255)
                                                                            ccgatctGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGaaa                                                                      <  1:59693/142‑1 (MQ=255)
                                                                                  tGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAg                                                               <  4:169276/149‑1 (MQ=255)
                                                                                    ccAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGAGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGaaa                                                                      >  2:59693/1‑140 (MQ=255)
                                                                                     cAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTc                                                            >  7:85797/1‑149 (MQ=255)
                                                                                      aGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGAGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTcc                                                           >  7:306949/1‑149 (MQ=255)
                                                                                         cGATCATGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGACTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAg                                                        >  6:324310/1‑149 (MQ=255)
                                                                                              ttgttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAg                                                   >  1:84260/1‑149 (MQ=255)
                                                                                                gttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAgg                                                  <  1:186623/148‑1 (MQ=255)
                                                                                                 ttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACgg                                                                   <  7:203253/130‑1 (MQ=255)
                                                                                                 ttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACgg                                                                   >  8:203253/1‑130 (MQ=255)
                                                                                                 ttgGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGgt                                                >  1:26205/1‑149 (MQ=255)
                                                                                                   ggCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGc                                                                             >  2:166461/1‑118 (MQ=255)
                                                                                                      cGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGtttt                                           >  7:18561/1‑149 (MQ=255)
                                                                                                          gAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                        <  2:84260/148‑1 (MQ=255)
                                                                                                                 cATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGGCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAAc                                <  7:218683/149‑1 (MQ=255)
                                                                                                                  aTTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACt                               >  4:22347/1‑149 (MQ=255)
                                                                                                                      atgaCGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGt                            >  5:315599/1‑148 (MQ=255)
                                                                                                                      atgaCGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGt                            <  3:22347/148‑1 (MQ=255)
                                                                                                                        gaCGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCAc                                                                           <  2:169939/99‑1 (MQ=255)
                                                                                                                           gAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAgg                      >  1:160099/1‑149 (MQ=255)
                                                                                                                           gAAAAAGGTGTGATCGGCACCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAg                       >  5:279098/1‑148 (MQ=255)
                                                                                                                              aaaGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCtt                   <  8:18561/149‑1 (MQ=255)
                                                                                                                               aaGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTg                  >  1:187810/1‑149 (MQ=255)
                                                                                                                                 ggTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGc                <  2:26205/149‑1 (MQ=255)
                                                                                                                                   tgtgATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGcc               >  5:174818/1‑148 (MQ=255)
                                                                                                                                           gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt      >  2:217090/1‑149 (MQ=255)
                                                                                                                                           gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt      <  6:355555/149‑1 (MQ=255)
                                                                                                                                           gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt      >  5:355555/1‑149 (MQ=255)
                                                                                                                                               cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCGTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTa  <  5:324310/149‑1 (MQ=255)
                                                                                                                                               cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                        <  3:198416/111‑1 (MQ=255)
                                                                                                                                               cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                        >  4:198416/1‑111 (MQ=255)
                                                                                                                                               cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTa  <  1:217090/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
ACCCGTCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTA  >  NC_002947/4677749‑4678040

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: