Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,326,822 T→C L38L (CTA→CTG ftsH ← ATP‑dependent zinc metalloprotease FtsH

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,326,8220TC100.0% 237.5 / NA 65L38L (CTA→CTGftsHATP‑dependent zinc metalloprotease FtsH
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (38/27);  total (38/27)

TTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTT  >  NC_000913/3326685‑3326958
                                                                                                                                         |                                                                                                                                        
ttGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAg                                                                                                                                        >  2:42839/1‑140 (MQ=255)
  ggTCAACAGGTTATCCATTAATGTCGGATCCTGAACCGGGATGTAAGTTGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGa                                                                                                                                      <  2:392426/140‑1 (MQ=255)
     cAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAgg                                                                                                                                   <  2:892082/140‑1 (MQ=255)
       aCAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTa                                                                                                                                 <  1:789264/140‑1 (MQ=255)
         aGGTTATCCAGTAATGTCGGATCCTTAACCGGAATGTGAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTaga                                                                                                                               <  1:499638/140‑1 (MQ=255)
                     aaTTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCtt                                                                                                                   >  2:217397/1‑140 (MQ=255)
                       tttCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTAc                                                                                                                 >  1:16893/1‑140 (MQ=255)
                         tCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACgg                                                                                                               >  2:524805/1‑140 (MQ=255)
                         tCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACATGGTCGTTATTGCCCTCTTGCAGGAAGGTAGAGTAATCCACCTTACgg                                                                                                               >  2:351863/1‑140 (MQ=255)
                              tCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGTAAGGTAGAGTAATACACCTTACGGCCAtt                                                                                                          >  1:277346/1‑140 (MQ=255)
                                  gAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGAc                                                                                                      >  2:48941/1‑140 (MQ=255)
                                   aaCCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACt                                                                                                     >  1:911094/1‑140 (MQ=255)
                                       ggAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCt                                                                                                 <  2:108959/140‑1 (MQ=255)
                                         aaTGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTgg                                                                                               >  2:529265/1‑140 (MQ=255)
                                             tAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGccc                                                                                           >  2:917486/1‑140 (MQ=255)
                                                   gTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCa                                                                                          >  2:885063/1‑135 (MQ=255)
                                                   gTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCa                                                                                          <  1:885063/135‑1 (MQ=255)
                                                      taACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTg                                                                                  >  1:255525/1‑140 (MQ=255)
                                                       aaCGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCATCTTACGGCCATTAGTCTCGCTGGGCCCAAAGCTCTgg                                                                                 >  1:463433/1‑140 (MQ=255)
                                                        aCGGTTACTATCTTTCTCGGTATCGTTGATTACACCTCCCTTGATACGTGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGa                                                                                <  2:277346/140‑1 (MQ=255)
                                                         cGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCCTACGGCCGTTAAACCCGCTGGTCCCAAAGCTCTGGaa                                                                               >  1:398745/1‑140 (MQ=255)
                                                              aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTg                                                                          >  1:15176/1‑140 (MQ=255)
                                                              aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTg                                                                          >  2:110615/1‑140 (MQ=255)
                                                              aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTg                                                                          >  1:995128/1‑140 (MQ=255)
                                                              aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTg                                                                          >  2:608378/1‑140 (MQ=255)
                                                              aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTg                                                                          >  1:88273/1‑140 (MQ=255)
                                                              aCTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTg                                                                          >  1:779542/1‑140 (MQ=255)
                                                                tATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGAc                                                                        <  2:936885/140‑1 (MQ=255)
                                                                tATCTTTCATGGTAACGTTGATTTCACGTCCGTCGCTACCCGCTTCACCAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGAc                                                                        <  2:485798/140‑1 (MQ=255)
                                                                          ggTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACaa                                                              <  2:16893/140‑1 (MQ=255)
                                                                          ggTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACaa                                                              >  2:733892/1‑140 (MQ=255)
                                                                              aCGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGc                                                          >  2:853323/1‑140 (MQ=255)
                                                                              aCGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGc                                                          >  2:83491/1‑140 (MQ=255)
                                                                              aCGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCGGGGGCCAAAGCTCTGGGATACAGACATCAGCACAACGGc                                                          >  2:277149/1‑140 (MQ=255)
                                                                               cGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCa                                                         >  2:527022/1‑140 (MQ=255)
                                                                                gTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCaa                                                        >  2:722139/1‑140 (MQ=255)
                                                                                   gATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGa                                                     <  2:15176/140‑1 (MQ=255)
                                                                                   gATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGa                                                     <  2:88273/140‑1 (MQ=255)
                                                                                   gATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGa                                                     <  2:730934/140‑1 (MQ=255)
                                                                                    aTTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTTACATCAGCACAACGGCAATGAc                                                    >  1:770453/1‑140 (MQ=255)
                                                                                      ttCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAcca                                                  >  1:664828/1‑140 (MQ=255)
                                                                                          cGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                              <  2:995128/140‑1 (MQ=255)
                                                                                          cGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                              <  2:255525/140‑1 (MQ=255)
                                                                                          cGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                              <  1:529265/140‑1 (MQ=255)
                                                                                              cGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGt                                          >  2:616443/1‑140 (MQ=255)
                                                                                                     aCGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                              <  1:510217/129‑1 (MQ=255)
                                                                                                     aCGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGAccagcca                                              >  2:510217/1‑129 (MQ=255)
                                                                                                     aCGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGt                                   <  1:547932/140‑1 (MQ=255)
                                                                                                                       gTCGTTATTGACCTCTTGCAGGAAGGTAGCGTAATCCCCCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCa                 <  1:277149/140‑1 (MQ=255)
                                                                                                                       gTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTTACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCa                 <  2:770453/140‑1 (MQ=255)
                                                                                                                       gTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCa                 >  1:947803/1‑140 (MQ=255)
                                                                                                                       gTCGTTATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCa                 <  1:83491/140‑1 (MQ=255)
                                                                                                                            tATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGa            >  1:981887/1‑140 (MQ=255)
                                                                                                                            tATTGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGa            >  2:642976/1‑140 (MQ=255)
                                                                                                                              ttGACCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGAtt          <  2:664828/140‑1 (MQ=255)
                                                                                                                                 aCCTCTTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGCTTTTTCGCCATGTCACTCAAGGGATTAAc       >  1:462850/1‑140 (MQ=255)
                                                                                                                                  ccTCTTGCAGGAAGTAAGAGTAATCCCCCTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct     <  1:841491/140‑1 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTAGAGTAATCCACCTTGCGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct     <  1:616443/140‑1 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGGCCAGCCAGGGTAGTAGCTTTTTCGCGATGTGACTCAAGGGATTAACct     >  2:275432/1‑140 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct     <  1:110615/140‑1 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct     <  1:42839/140‑1 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct     <  2:911094/140‑1 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCAATGTCACTCAAGGGATTAACct     >  2:316216/1‑140 (MQ=255)
                                                                                                                                   ctctTGCAGGAAGGTACGGTAATCCACCTTACGACGATTTGACTCGCTAGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACct     <  2:449299/140‑1 (MQ=255)
                                                                                                                                     ctTGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACctct   >  1:524337/1‑140 (MQ=255)
                                                                                                                                      ttGCAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCtt  >  2:252022/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                        
TTGGTCAACAGGTTATCCAGTAATTTCGGATCCTGAACCGGAATGTAAGTGGTATAACGGTTACTATCTTTCTTGGTAACGTTGATTTCACGTCCGTTGATACGCGCTTCACGAACCTGGTCGTTATTGACCTCTTGTAGGAAGGTAGAGTAATCCACCTTACGGCCATTAGACTCGCTGGGCCCAAAGCTCTGGAATACTGACATCAGCACAACGGCAATGACCAGCCAGAGTATTAGGTTTTTCGCCATGTCACTCAAGGGATTAACCTCTT  >  NC_000913/3326685‑3326958

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: