Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,677,900 C→G 6.2% R225P (CGC→CCC) ‡ ldcC ← lysine decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,677,9000CG 6.2% 203.8 / 4.0 65R225P (CGC→CCC) ‡ldcClysine decarboxylase
Reads supporting (aligned to +/- strand):  ref base C (30/30);  new base G (2/2);  total (33/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.03e-01

TCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTATGATAGGC  >  NC_002947/4677754‑4678048
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                                               <  8:46139/149‑1 (MQ=255)
                     ccGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCg                                                                                                                               <  5:120733/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                                                                                                                   <  2:294929/137‑1 (MQ=255)
                                             ccagcaccagGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCt                                                                                                                   >  1:294929/1‑137 (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              <  6:355555/149‑1 (MQ=255)
                                                                                                                                      gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt              >  2:217090/1‑149 (MQ=255)
                                                                                                                                      gcgcCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCt              >  5:355555/1‑149 (MQ=255)
                                                                                                                                          cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCGTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTa          <  5:324310/149‑1 (MQ=255)
                                                                                                                                          cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                                >  4:198416/1‑111 (MQ=255)
                                                                                                                                          cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTccc                                                <  3:198416/111‑1 (MQ=255)
                                                                                                                                          cGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTa          <  1:217090/149‑1 (MQ=255)
                                                                                                                                           gAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAt         <  6:315599/149‑1 (MQ=255)
                                                                                                                                           gAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTAt         <  8:85797/149‑1 (MQ=255)
                                                                                                                                            aaGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCggg                >  8:19110/1‑141 (MQ=255)
                                                                                                                                            aaGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCggg                <  7:19110/141‑1 (MQ=255)
                                                                                                                                                  cgcgcgGCCCTGGGCTCGGCTTTGTCTAGCGGGCCTGTGTGGTCGAGGAGCGAACCCAGCTCGGGGACGGAAACGGACAGGTCCGAGCGCAGGTTGTTTTCCCAGAAGAACTGGTGAAAGGCTTGGCGTACCGGGCTTTTATGATAGGc  <  8:306949/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCATGATGATCGCGTGCACCACCGACTTGTGGCAGTTGCGGTCCACCAGCACCAGGTCGTCACGGCCGACCATGGCATGCCAGACGATCTTGTTGGCCGTGGAGGTGCCATTGATGACGAAAAAGGTGTGATCGGCGCCGAAGTTGCGCGCCGCCCTGGCCTCGGCTTCGGCCAGCGGGCCTGTGTGGTCGAGCAGCGAACCCAGCTCCGGCACGGAAACGGACAGGTCCGAGCGCAGGGTGTTTTCCCCGAAGAACTGGTGAAAGGCTTGGCCTACCGGGCTTTTATGATAGGC  >  NC_002947/4677754‑4678048

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: