Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,071,211 T→A 7.8% C22S (TGC→AGC)  yiaY → Fe‑containing alcohol dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,071,2110TA 7.8% 172.1 / 3.5 60C22S (TGC→AGC) yiaYFe‑containing alcohol dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (25/30);  new base A (2/3);  total (27/33)
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.53e-01

CCCCGCTTTGGTCTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGTATCGACTACTGCCTGT  >  NC_002947/3071066‑3071352
                                                                                                                                                  |                                                                                                                                             
ccccGCTTTGGTCTGGTTGCTTTCGTGGCGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGtgccggctgcc                                                                                                                                            >  2:276945/1‑149 (MQ=255)
  ccGCTTTGGTCTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGtgccggctgccgg                                                                                                                                          >  4:173492/1‑149 (MQ=255)
     cTTTGGTCTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCAGCCGGcac                                                                                                                                       <  5:78234/149‑1 (MQ=255)
         ggTCTGGTTGCTTTCGTGGTGGAATGCTTTACCGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCATGCGCAAGTTCGTATCGCCTGAATTCATATTTGTTGCCGGATGCCGGCACAAt                                                                                                                                    >  1:231986/1‑149 (MQ=255)
          gTCTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAAtgt                                                                                                                                  <  4:189983/149‑1 (MQ=255)
          gTCTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCAGCCGGcaccaaa                                                                                                                                   <  2:378248/148‑5 (MQ=255)
          gTCTGGTTGCTTTCGTGGGGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAAtgt                                                                                                                                  <  1:291108/149‑1 (MQ=255)
            cTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTCTCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTgg                                                                                                                                <  5:119715/149‑1 (MQ=255)
               ggtGCTCTCGTGGGGGGATGCTTTA‑CGCTGGCGGTTATCTTCCACAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGcc                                                                                                                              <  6:209840/146‑1 (MQ=255)
               gTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGTCGGCTGCCGGCACAAtgtcacgc                                                                                                                             >  1:324448/1‑144 (MQ=255)
               gTTGCTTTCGTGGTGGGATGCTTTA‑CGCTCGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATATTTCGTGCCGGCTGCCGGCACAAAGTGGcc                                                                                                                              >  5:80537/1‑148 (MQ=255)
                tttcTTTCGGGGTGGGATGCTTTA‑CTCTGGCTGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCTCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCaa                                                                                                                            <  6:80537/146‑1 (MQ=255)
                      tCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCCGCCCCCTTCGCAAGTTCGTATAGCGTGAAATCATCTTTGGTGCCGGCAGCGGGcac                                                                                                                                       >  6:280273/1‑132 (MQ=255)
                      tCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACg                                                                                                                       <  3:173492/148‑1 (MQ=255)
                      tCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATCAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCAGCCGGcac                                                                                                                                       <  5:280273/132‑1 (MQ=255)
                           gtggGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCaaaa                                                                                                                 >  3:158989/1‑149 (MQ=255)
                            tggGATGCTGTA‑CGCTGGCGCTTATCATCCAGAACAATAAGAACCGCCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAAc                                                                                                                <  8:101393/149‑1 (MQ=255)
                                aTGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTc                                                                                                            >  3:39454/1‑149 (MQ=255)
                                   cTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCgg                                                                                                          <  2:231986/148‑1 (MQ=255)
                                   cTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGt                                                                                                         <  4:39454/149‑1 (MQ=255)
                                          gCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGcc                                                                                                                     >  7:271611/1‑131 (MQ=255)
                                          gCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGcc                                                                                                                     <  8:271611/131‑1 (MQ=255)
                                           cTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTgcgcgc                                                                                                   >  1:3176/1‑148 (MQ=255)
                                              gCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAgg                                                                                               <  2:3176/149‑1 (MQ=255)
                                                ggTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCATATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTa                                                                                             >  5:68360/1‑149 (MQ=255)
                                                    aTCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACtg                                                                                          <  4:158989/148‑1 (MQ=255)
                                                     tCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACtgg                                                                                         <  5:311880/148‑1 (MQ=255)
                                                        tCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTggtggt                                                                                     >  1:363690/1‑149 (MQ=255)
                                                         ccAGAACAATAAGAACCGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACtggtg                                                                                       <  6:60721/146‑1 (MQ=255)
                                                          cAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCTGTGCGCGCAAGGTACTGGTGGTCa                                                                                   <  5:371351/149‑1 (MQ=255)
                                                          cAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTggtggt                                                                                     >  1:37827/1‑147 (MQ=255)
                                                           aGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAg                                                                                  >  6:1280/1‑149 (MQ=255)
                                                           aGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAAACATCTTTGGTGCCGGCTGCCGCCACAATGTGGCCAATTACGCCAAAACCTTCGCTGCGCGCAAGGAACTGGTGGTCAg                                                                                  >  8:151596/1‑149 (MQ=255)
                                                            gAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACtg                                                                                          >  5:300317/1‑140 (MQ=255)
                                                            gAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACtg                                                                                          <  6:300317/140‑1 (MQ=255)
                                                               cAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGAc                                                                              >  6:177277/1‑149 (MQ=255)
                                                                 aTAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGAccc                                                                            >  7:30171/1‑149 (MQ=255)
                                                                   aaGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCCGTGCGCGCAAGGTACtgg                                                                                         <  8:157042/134‑1 (MQ=255)
                                                                   aaGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCAGTGCGCGCAAGGTACtgg                                                                                         >  7:157042/1‑134 (MQ=255)
                                                                    aGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTgg                                                                         >  8:50583/1‑149 (MQ=255)
                                                                    aGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTgg                                                                         <  5:107025/149‑1 (MQ=255)
                                                                           gCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATc                                                                  <  1:276945/149‑1 (MQ=255)
                                                                           gCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATc                                                                  >  1:281603/1‑149 (MQ=255)
                                                                                        ccAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTggg                                                     >  4:194012/1‑149 (MQ=255)
                                                                                        ccAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTggg                                                     >  8:130216/1‑149 (MQ=255)
                                                                                              tttCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCg                                                >  8:185985/1‑148 (MQ=255)
                                                                                                tCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGa                                              <  8:30171/148‑1 (MQ=255)
                                                                                                 cAGCCCCCTTCGCAAGTTCATATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGAt                                             <  6:68360/148‑1 (MQ=255)
                                                                                                   gCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGAtgtg                                          <  7:97536/149‑1 (MQ=255)
                                                                                                    cccccTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGAtgtg                                          >  8:159939/1‑148 (MQ=255)
                                                                                                      cccTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAg                                       <  5:1280/149‑1 (MQ=255)
                                                                                                           cGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCa                                   <  7:185985/148‑1 (MQ=255)
                                                                                                             cAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGcc                                >  6:243453/1‑149 (MQ=255)
                                                                                                              aaGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCt                               >  7:183115/1‑149 (MQ=255)
                                                                                                                         gCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGTTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAgg                    <  2:281603/149‑1 (MQ=255)
                                                                                                                         gCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGctggcc                    <  2:37827/149‑6 (MQ=255)
                                                                                                                         gCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAgg                    <  7:50583/149‑1 (MQ=255)
                                                                                                                           cTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAggg                   >  4:133229/1‑148 (MQ=255)
                                                                                                                           cTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCGGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGt                  <  8:183115/149‑1 (MQ=255)
                                                                                                                                        ttgTGCCGGCAGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGTATCGACTACTGcc     >  3:37513/3‑149 (MQ=255)
                                                                                                                                        tGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGTATCGactact        <  3:133229/146‑1 (MQ=255)
                                                                                                                                        tGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGTATCGACTACTGc      <  7:130216/148‑1 (MQ=255)
                                                                                                                                           tgccggctgccggcACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGTATCGACTACTGCCTGt  <  2:363690/149‑1 (MQ=255)
                                                                                                                                             ccggctgccggcACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAgg                    >  5:105459/1‑129 (MQ=255)
                                                                                                                                             ccggctgccggcACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAgg                    <  6:105459/129‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                             
CCCCGCTTTGGTCTGGTTGCTTTCGTGGTGGGATGCTTTA‑CGCTGGCGGTTATCTTCCAGAACAATAAGAACAGGCTTACCGATGAGCCAGAGTTTCAGCCCCCTTCGCAAGTTCGTATCGCCTGAAATCATCTTTGGTGCCGGCTGCCGGCACAATGTGGCCAATTACGCCAAAACCTTCGGTGCGCGCAAGGTACTGGTGGTCAGCGACCCTGGCGTGATCGCCGCCGGCTGGGTGGCGGATGTGGAGGCCAGCCTGCAGGCCCAGGGTATCGACTACTGCCTGT  >  NC_002947/3071066‑3071352

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑