Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,830,229 T→A 6.3% D384V (GAC→GTC)  PP_3383 ← gluconate 2‑dehydrogenase flavoprotein subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,830,2290TA 6.3% 184.1 / 3.0 63D384V (GAC→GTC) PP_3383gluconate 2‑dehydrogenase flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base A (2/2);  total (30/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 = 4.60e-01

AGTAGTCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGAAGTTACGCCCGACCACACCCTCGTTCTTTACCGGGTC  >  NC_002947/3830083‑3830377
                                                                                                                                                  |                                                                                                                                                     
agtagtCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCg                                                                                                                                                     >  7:169975/1‑149 (MQ=255)
 gtagtCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCgg                                                                                                                                                    >  4:274803/1‑149 (MQ=255)
 gtagtCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCg                                                                                                                                                     <  6:274257/148‑1 (MQ=255)
 gtagtCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCg                                                                                                                                                     <  2:68234/148‑1 (MQ=255)
    gtCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGt                                                                                                                                                 >  1:202765/1‑149 (MQ=255)
    gtCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGt                                                                                                                                                 >  8:196504/1‑149 (MQ=255)
               ccgccctcCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAAtcgt                                                                                                                                        <  1:305722/141‑1 (MQ=255)
                tttCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGa                                                                                                                                    >  1:161208/1‑149 (MQ=255)
                tttCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGACGGCGTTGAAATCGTCGa                                                                                                                                    >  7:289863/1‑149 (MQ=255)
                tttCTTCCAGGCACTGCCACAGGCCGGAGTGACCGGCGGGTTGGAGGTACCTGCAATCGGCCGGCTGCCGGACTCGTTCACCCACATTGGCGAACCGCCGACGAAGACGTGCGGGACATGGTCGAAGTTGTCGGCCTTGAAATCGTCGa                                                                                                                                    >  7:201944/1‑149 (MQ=255)
                  tcttcCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGa                                                                                                                                    <  5:32320/147‑1 (MQ=255)
                      ccAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATgg                                                                                                                                 <  7:196504/146‑1 (MQ=255)
                       cAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAAt                                                                                                                             >  7:266592/1‑149 (MQ=255)
                        aGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATg                                                                                                                            >  3:20766/1‑149 (MQ=255)
                         ggCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGc                                                                                                                           <  2:161208/149‑1 (MQ=255)
                           cACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCa                                                                                                                         >  4:270371/1‑149 (MQ=255)
                           cACTGCCCCAGGCCAGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCa                                                                                                                         <  5:109874/149‑1 (MQ=255)
                                  ccAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGACGGCGTTGAAATCGTCGATGGCAATGCCATTGCCAc                                                                                                                  >  7:134284/1‑149 (MQ=255)
                                   cAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCAcc                                                                                                                 <  2:160235/149‑1 (MQ=255)
                                     ggCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGACGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCcg                                                                                                               <  8:134284/149‑1 (MQ=255)
                                             tGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGAt                                                                                                        >  1:270097/1‑148 (MQ=255)
                                               cccGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGaa                                                                                                     >  4:181600/1‑149 (MQ=255)
                                                    cggGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAgttgt                                                                                                <  2:202765/149‑1 (MQ=255)
                                                     gggTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAgttgtt                                                                                               <  8:266592/149‑1 (MQ=255)
                                                      ggTTTGAGGTACCTGCAATAGTCCGGGTGCCGGCCTGTTTCACCCAAATTGGCGAACCGCAGACGAAGCCGTGCTGGGCATGTTCGAAGTTTTCGGCTTTGAAATCGTCTATGGCAATGCCATTGCCACCCGGGCCGATGAAGttgttg                                                                                              <  8:201944/149‑1 (MQ=255)
                                                      ggTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGttgttg                                                                                              <  3:181600/149‑1 (MQ=255)
                                                          ggAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGgtg                                                                                           <  2:270097/148‑1 (MQ=255)
                                                           gAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAAtcgtc                                                                                                                                       <  2:73350/104‑1 (MQ=255)
                                                           gAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAAtcgtc                                                                                                                                       >  1:73350/1‑104 (MQ=255)
                                                           gAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGACGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTa                                                                                         <  8:289863/149‑1 (MQ=255)
                                                            aGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAg                                                                                        <  4:158469/149‑1 (MQ=255)
                                                            aGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAg                                                                                        <  4:219960/149‑1 (MQ=255)
                                                              gTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGgt                                                                                      >  4:268427/1‑149 (MQ=255)
                                                                 ccGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGgtgt                                                                                    <  3:268427/148‑1 (MQ=255)
                                                                  cGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGttgttg                                                                                              >  5:24856/1‑137 (MQ=255)
                                                                  cGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGttgttg                                                                                              <  6:24856/137‑1 (MQ=255)
                                                                   ggCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAgttgtt                                                                                               <  6:60816/135‑1 (MQ=255)
                                                                   ggCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAgttgtt                                                                                               >  5:60816/1‑135 (MQ=255)
                                                                          ggccggcTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCgaa                                                                          <  3:270371/149‑1 (MQ=255)
                                                                            ccggcTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCgaaga                                                                        <  5:87814/149‑1 (MQ=255)
                                                                            ccggcTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCgaag                                                                         >  8:59754/1‑148 (MQ=255)
                                                                              ggcTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAAc                                                                      >  8:296700/1‑149 (MQ=255)
                                                                              ggcTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAAc                                                                      <  7:296700/149‑1 (MQ=255)
                                                                                   ccGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTgg                                                                                             <  4:321180/121‑1 (MQ=255)
                                                                                   ccGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTgg                                                                                             >  3:321180/1‑121 (MQ=255)
                                                                                           gTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGcc                                                                   <  2:81958/139‑1 (MQ=255)
                                                                                           gTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGcc                                                                   >  1:81958/1‑139 (MQ=255)
                                                                                           gTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTAc                                                         >  2:45114/1‑149 (MQ=255)
                                                                                               aCCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGc                                                                    <  4:153442/134‑1 (MQ=255)
                                                                                               aCCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGc                                                                    >  3:153442/1‑134 (MQ=255)
                                                                                                 ccACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGt                                                   >  2:137283/1‑149 (MQ=255)
                                                                                                 ccACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGt                                                   >  5:261289/1‑149 (MQ=255)
                                                                                                  cacaTTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGtt                                                  >  7:162757/1‑149 (MQ=255)
                                                                                                      ttGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCACTTCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGa                                              <  7:230102/149‑1 (MQ=255)
                                                                                                        ggCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGAt                                             <  3:274803/148‑1 (MQ=255)
                                                                                                           gAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGc                                         <  6:261289/149‑1 (MQ=255)
                                                                                                              ccgccgACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGaa                                      <  4:20766/149‑1 (MQ=255)
                                                                                                              ccgccgACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGaa                                      <  8:169975/149‑1 (MQ=255)
                                                                                                                    acgaAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCAAGTTCTGATAGGCGAAGTTACg                                >  6:372571/1‑149 (MQ=255)
                                                                                                                     cgaAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGAAGTTACGc                               <  8:162757/149‑1 (MQ=255)
                                                                                                                      gaAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGAATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGt                                                           >  1:78207/1‑121 (MQ=255)
                                                                                                                        aGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGAAGTTACGccc                             >  2:304900/1‑148 (MQ=255)
                                                                                                                                           gAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGAAGTTACGCCCGACCACACCCTCGTTCTTTa         >  1:5624/1‑149 (MQ=255)
                                                                                                                                                  tCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGAAGTTACGCCCGACCACACCCTCGTTCTTTACCGGGTc  >  2:182877/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
AGTAGTCGGCCGTGGCTTTCTTCCAGGCACTGCCCCAGGCCGGAGTGCCCGGCGGGTTGGAGGTACCGGCAATCGGCCGGCTGCCGGCCTGGTTCACCCACATTGGCGAACCGCCGACGAAGCCGTGCGGGCCATGGTCGAAGTTGTCGGCGTTGAAATCGTCGATGGCAATGCCATTGCCACCCGCGCCGATGAAGTTGTTGGTGTAGGTGTCCTTGTCGAAGAACGCCTTCACTGTACCCATGTTCTGATAGGCGAAGTTACGCCCGACCACACCCTCGTTCTTTACCGGGTC  >  NC_002947/3830083‑3830377

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: