Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,306,095 C→T 6.0% R168H (CGT→CAT)  yjcS ← metallo‑beta‑lactamase superfamily protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,306,0950CT 6.0% 202.4 / 2.5 69R168H (CGT→CAT) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (30/32);  new base T (2/2);  total (35/34)
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.72e-01

GATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTG  >  NC_000913/4305950‑4306241
                                                                                                                                                 |                                                                                                                                                  
gATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGAt                                                                                                                                                 >  3:29603/1‑149 (MQ=255)
    gCTTCGTCCATAAAGCCTGCAGTGGCAATCACCTGAACTTTGCCGGATTTAAGATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTTTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg                                                                                                                                             <  3:45676/149‑1 (MQ=255)
     cTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTgg                                                                                                                                            <  4:332/149‑1 (MQ=255)
         gTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCTGACGATGCTGGAAg                                                                                                                                         >  6:193804/1‑148 (MQ=255)
          tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGt                                                                                                                                        >  7:6084/1‑148 (MQ=255)
          tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGt                                                                                                                                        >  6:362143/1‑148 (MQ=255)
                 aGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGATTTTTCTTACTATACTGTAAGTAAAGGTc                                                                                                                                 >  2:382677/1‑148 (MQ=255)
                  gCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCaa                                                                                                                               <  8:256089/149‑1 (MQ=255)
                            gCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACCAACGGATAATTAGTACGATCCTGGAAGTAAAGGTCAAGTGCGGCttt                                                                                                                     >  8:233049/1‑149 (MQ=255)
                               aTCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTATGATGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGc                                                                                                                  >  4:107696/1‑149 (MQ=255)
                                tCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCg                                                                                                                 <  5:335252/149‑1 (MQ=255)
                                      gAACTTTGACGGAATTAACATCGGCTTCAGAGATAAAGCCTTTCACGCCACCATAGTGGTCGCTGTCGCTGTGAGTGTAGAAAACGGCAACAATCGCTTTATGCAGACGATGCTCGAAGTAAAGGTCAAGTCCGGCGTTAGCGGAAggc                                                                                                           >  8:385567/1‑149 (MQ=255)
                                       aaCTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAggcg                                                                                                          <  6:358307/149‑1 (MQ=255)
                                          tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAg                                                                                                                                    >  7:351739/1‑120 (MQ=255)
                                          tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                                       >  5:340661/1‑149 (MQ=255)
                                          tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                                                                       >  8:394856/1‑149 (MQ=255)
                                          tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGAGGACGATGCTGGAAGTAAAg                                                                                                                                    <  8:351739/120‑1 (MQ=255)
                                             gCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTc                                                                                                     <  7:394856/148‑1 (MQ=255)
                                              ccGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTTACGACGATCCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAc                                                                                                   >  8:179250/1‑149 (MQ=255)
                                               cGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCAc                                                                                                   >  4:120707/1‑148 (MQ=255)
                                                  aTTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGc                                                                                               <  3:8237/149‑1 (MQ=255)
                                                      aaCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGAAACAAACGGATATTTAAAACAAAGCTGGAAGTAAAGGTAAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGt                                                                                           >  4:91529/1‑149 (MQ=255)
                                                       agaTCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACTATCGGTTTTTGGGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGt                                                                                           <  6:371830/146‑1 (MQ=255)
                                                       aCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGCAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTc                                                                                          <  5:193804/149‑1 (MQ=255)
                                                       aCATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGGGTGGCTGTGAGTGTAGATAACGGCAACAATCGGCTTTTGCGGACGATGCGGGAAGTAAAGGTCAAGGGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGt                                                                                           <  7:179250/148‑1 (MQ=255)
                                                          tCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGAt                                                                                       <  1:195988/149‑1 (MQ=255)
                                                            ggCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATaa                                                                                     <  5:362143/149‑1 (MQ=255)
                                                             gCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAAACGGTTTTTTATGAAGATGATGGAAGTAAAGGTCAAGTGCTGCTTTTGCGTCAGGCGGCGTCACCAGCGGGTCGATaa                                                                                     >  2:125969/1‑148 (MQ=255)
                                                               ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                                                              >  3:134494/1‑121 (MQ=255)
                                                               ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                                                              <  4:134494/121‑1 (MQ=255)
                                                               ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                                                              >  3:33586/1‑121 (MQ=255)
                                                               ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                                                              <  4:33586/121‑1 (MQ=255)
                                                                 cAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaa                                                                                  <  3:8713/147‑1 (MQ=255)
                                                                  agagATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                               <  3:107696/149‑1 (MQ=255)
                                                                    agaTAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTATAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATg                                                                             <  3:316534/149‑1 (MQ=255)
                                                                        aaTGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCaa                                                                                                                               <  4:310467/95‑1 (MQ=255)
                                                                        aaTGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCaa                                                                                                                               >  3:310467/1‑95 (MQ=255)
                                                                         aTGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCttt                                                                        <  4:29603/149‑1 (MQ=255)
                                                                             cTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTGTCCAGCGGGTCGATAACAATAATGCCTTTCTCg                                                                    <  6:340661/149‑1 (MQ=255)
                                                                                    gCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaa                                                             >  7:169694/1‑149 (MQ=255)
                                                                                      caccaTAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTGAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaac                                                            >  7:80520/1‑148 (MQ=255)
                                                                                        ccaTAGTGGTCGGTGTGGCTGTGAGTGTAGATCACGGAAAACGACGTATAATTAAAAAAATACTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGCGGGGCGATAACAATAATGCCTTTCTCTCCCTCaacgac                                                         >  8:137893/1‑148 (MQ=255)
                                                                                         caTAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAAAAAACGGTTTTTGAGAAAAATTATTGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCATGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCaacgaac                                                        >  6:299098/1‑149 (MQ=255)
                                                                                             gTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTa                                                     <  8:6084/148‑1 (MQ=255)
                                                                                              tGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt                                                   <  8:169694/149‑1 (MQ=255)
                                                                                               ggTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATg                                                  >  5:270739/1‑149 (MQ=255)
                                                                                                gTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGt                                                 <  1:382677/149‑1 (MQ=255)
                                                                                                    gtgtGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCCACGAACGTAATGTTaga                                             <  7:233049/149‑1 (MQ=255)
                                                                                                       tgGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTGAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGat                                          <  8:80520/149‑1 (MQ=255)
                                                                                                                gtgtAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTaa                                                                                                                                      <  7:393055/48‑1 (MQ=255)
                                                                                                                gtgtAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTaa                                                                                                                                      >  8:393055/1‑48 (MQ=255)
                                                                                                                 tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt                                                   <  3:149157/130‑1 (MQ=255)
                                                                                                                 tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt                                                   >  6:82185/1‑130 (MQ=255)
                                                                                                                 tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt                                                   <  5:82185/130‑1 (MQ=255)
                                                                                                                 tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt                                                   >  4:149157/1‑130 (MQ=255)
                                                                                                                 tgtAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                               <  4:14106/102‑1 (MQ=255)
                                                                                                                    aGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                               >  4:284737/1‑99 (MQ=255)
                                                                                                                    aGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                                                               <  3:284737/99‑1 (MQ=255)
                                                                                                                            gCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATa                     >  7:397383/1‑149 (MQ=255)
                                                                                                                            gCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATa                     >  8:243684/1‑149 (MQ=255)
                                                                                                                             caacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCCCCAGGGGGTCGATTACAATAATGCCGTTCTCGCCCTCAACGAGCGTAATGTTAGAGATATCCTGGCCGCGCACCTGATAc                    >  8:390233/1‑149 (MQ=255)
                                                                                                                              aacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACa                   >  8:236633/1‑149 (MQ=255)
                                                                                                                              aacaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACa                   >  5:245890/1‑149 (MQ=255)
                                                                                                                               acaaTCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGcc                                                                  >  2:213243/1‑101 (MQ=255)
                                                                                                                                 aaTCGGTTTTTGCGGACGATTCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACAttt                <  7:385567/149‑1 (MQ=255)
                                                                                                                                   tCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTa              >  8:358102/1‑149 (MQ=255)
                                                                                                                                      gTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCg           >  7:348113/1‑149 (MQ=255)
                                                                                                                                              ggACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACttt   <  6:270739/149‑1 (MQ=255)
                                                                                                                                               gACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTg  <  7:243684/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
GATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGATATCCTGGCCGCGCACCTGATACATTTTATCGGTGACTTTG  >  NC_000913/4305950‑4306241

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