Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,306,072 C→T 5.7% V176I (GTT→ATT)  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,0720CT 5.7% 225.0 / 2.4 70V176I (GTT→ATT) yjcSmetallo‑beta‑lactamase superfamily protein
Reads supporting (aligned to +/- strand):  ref base C (30/36);  new base T (2/2);  total (32/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.60e-01

ATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGAT  >  NC_000913/4305924‑4306201
                                                                                                                                                    |                                                                                                                                 
atgTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAAc                                                                                                                                   >  4:150667/1‑149 (MQ=255)
   ttACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACgg                                                                                                                                 <  5:278047/148‑1 (MQ=255)
   ttACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACgg                                                                                                                                 <  2:363653/148‑1 (MQ=255)
   ttACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACgg                                                                                                                                 <  7:140166/148‑1 (MQ=255)
      cccGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGcaa                                                                                                                              >  4:8713/1‑148 (MQ=255)
            aGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCgg                                                                                                                       <  1:91991/149‑1 (MQ=255)
                 gTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGttttt                                                                                                                  >  6:335252/1‑149 (MQ=255)
                 gTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGttttt                                                                                                                  >  7:256089/1‑149 (MQ=255)
                    ttCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGc                                                                                                                <  3:150667/148‑1 (MQ=255)
                     tCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCg                                                                                                               <  3:115654/148‑1 (MQ=255)
                      cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATc                                                                                                                         >  7:63089/1‑137 (MQ=255)
                      cGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATc                                                                                                                         <  8:63089/137‑1 (MQ=255)
                          gATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGAt                                                                                                         >  3:29603/1‑149 (MQ=255)
                           aTGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTAAGAGAGAATGCCTTTCACGCCACCATAGTGGTCGGTGTGCCTGTGAGTGTAGCTAACGGCAACAATCTGTTTATGCgg                                                                                                              >  5:371830/1‑143 (MQ=255)
                              gCTTCGTCCATAAAGCCTGCAGTGGCAATCACCTGAACTTTGCCGGATTTAAGATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTTTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTg                                                                                                     <  3:45676/149‑1 (MQ=255)
                               cTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTgg                                                                                                    <  4:332/149‑1 (MQ=255)
                                   gTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCTGACGATGCTGGAAg                                                                                                 >  6:193804/1‑148 (MQ=255)
                                    tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGt                                                                                                >  6:362143/1‑148 (MQ=255)
                                    tCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGt                                                                                                >  7:6084/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                                                               >  8:394856/1‑149 (MQ=255)
                                                                    tttGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGgcggcg                                                               >  5:340661/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                                                                      <  4:33586/121‑1 (MQ=255)
                                                                                         ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                      <  4:134494/121‑1 (MQ=255)
                                                                                         ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                      >  3:33586/1‑121 (MQ=255)
                                                                                         ttCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCa                                                                      >  3:134494/1‑121 (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           >  6:82185/1‑130 (MQ=255)
                                                                                                                                           tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt           <  5:82185/130‑1 (MQ=255)
                                                                                                                                           tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt           >  4:149157/1‑130 (MQ=255)
                                                                                                                                           tgtAGATAATGGCAACAATCGGTTTTTGCGGATGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAAt           <  3:149157/130‑1 (MQ=255)
                                                                                                                                           tgtAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                       <  4:14106/102‑1 (MQ=255)
                                                                                                                                              aGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                       >  4:284737/1‑99 (MQ=255)
                                                                                                                                              aGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACaataa                                       <  3:284737/99‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                 
ATGTTACCCGCCAGCACGTTTTCGCTGATGGCTTCGTCCATAAAGCCTGCAGGGGCAATCACCTGAACTTTGCCGGATTTAACATCGGCTTCAGAGATAATGCCTTTCACGCCACCATAGTGGTCGGTGTGGCTGTGAGTGTAGATAACGGCAACAATCGGTTTTTGCGGACGATGCTGGAAGTAAAGGTCAAGTGCGGCTTTTGCGGCAGGCGGCGTCACCAGCGGGTCGATAACAATAATGCCTTTCTCGCCCTCAACGAACGTAATGTTAGAGAT  >  NC_000913/4305924‑4306201

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: ‑