Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,168,391 C→T 8.5% G501D (GGC→GAC)  katG ← catalase‑peroxidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,168,3910CT 8.5% 163.2 / 3.2 59G501D (GGC→GAC) katGcatalase‑peroxidase
Reads supporting (aligned to +/- strand):  ref base C (25/28);  new base T (2/3);  total (27/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.81e-01

AGGGAAACCTTCTTGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGCTTTGGGCAACGG  >  NC_002947/4168244‑4168524
                                                                                                                                                   |                                                                                                                                     
aGGGAAACCTTCTTGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGcc                                                                                                                                      >  1:13774/1‑149 (MQ=255)
        cttcttGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAgg                                                                                                                              <  3:104874/149‑1 (MQ=255)
        cttcttGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTATGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAgg                                                                                                                              >  7:129154/1‑149 (MQ=255)
           cttGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGt                                                                                                                             <  1:8911/147‑1 (MQ=255)
           cttGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTATGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGt                                                                                                                             <  8:129154/147‑1 (MQ=255)
               ccGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGAcgc                                                                                                                       >  4:308487/1‑149 (MQ=255)
               ccGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGACCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGAcgc                                                                                                                       >  1:88571/1‑149 (MQ=255)
                cGCTGCAGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGACCACGCCCTGGTTTGCAGCCAAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCCCGGAAGGTCGAcgcg                                                                                                                      >  4:28600/1‑149 (MQ=255)
                     ctgTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGAt                                                                                                                   <  4:221539/147‑1 (MQ=255)
                       gTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGcc                                                                                                                >  1:133543/1‑148 (MQ=255)
                       gTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGcc                                                                                                                >  2:85946/1‑148 (MQ=255)
                          gAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAg                                                                                                             <  3:308487/148‑1 (MQ=255)
                           aaCTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAgg                                                                                                            <  5:21563/148‑1 (MQ=255)
                             cTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCgg                                                                                                         >  5:7184/1‑149 (MQ=255)
                                gCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTg                                                                                                       >  1:17907/1‑148 (MQ=255)
                                 cccTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGTACCACCGCGCTTGTCGGAGTCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGga                                                                                                     <  2:14428/149‑1 (MQ=255)
                                      gATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCa                                                                                                <  7:173050/149‑1 (MQ=255)
                                      gATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCa                                                                                                <  2:133543/149‑1 (MQ=255)
                                        tcttctCCAATGCCGCCAGCACCTTCTCCACCCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAACCGTGCGCCGTTGGCACCACCGCGCTTGTAGGAGCCGAGGAAGGTCGACGCGGATGCCCAGCCGGAGGAGACCAGt                                                                                              >  3:220686/1‑149 (MQ=255)
                                        tattctCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTACTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGt                                                                                              >  5:279804/3‑149 (MQ=255)
                                               cAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGTCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCa                                                                                       >  7:160400/1‑149 (MQ=255)
                                               cAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCa                                                                                                >  1:95358/1‑140 (MQ=255)
                                               cAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCa                                                                                                <  2:95358/140‑1 (MQ=255)
                                               cAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGc                                                                                          <  1:44615/146‑1 (MQ=255)
                                                aaTGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGTCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGccctc                                                                                      >  6:172415/1‑147 (MQ=255)
                                                      gccAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCa                                                                                 >  2:16939/1‑148 (MQ=255)
                                                            aCCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAAccc                                                                            >  5:207977/1‑147 (MQ=255)
                                                             ccTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAg                                                                         <  6:279804/149‑1 (MQ=255)
                                                              cTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAgg                                                                        >  2:138419/1‑149 (MQ=255)
                                                                    cACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGa                                                                  >  7:50695/1‑149 (MQ=255)
                                                                             gTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGa                                                           <  8:50695/147‑1 (MQ=255)
                                                                                 gCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGa                                                      <  2:17907/148‑1 (MQ=255)
                                                                                  cAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAgg                                                    <  2:13774/149‑1 (MQ=255)
                                                                                       cAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTcgc                                               <  1:138419/149‑1 (MQ=255)
                                                                                       cAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTcgc                                               >  3:275617/1‑149 (MQ=255)
                                                                                        aGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTcgc                                               >  1:158273/1‑148 (MQ=255)
                                                                                        aGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTcgc                                               >  5:37946/1‑148 (MQ=255)
                                                                                          tCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCACAAACCCGAGTCCAGGATCCTGGACTTGAGGGTCGCGAt                                            <  3:28600/149‑1 (MQ=255)
                                                                                               cTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTc                                         <  1:16939/147‑1 (MQ=255)
                                                                                                tgtgGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTg                                      >  4:304/1‑149 (MQ=255)
                                                                                                tgtgGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTg                                      >  6:142801/1‑149 (MQ=255)
                                                                                                  tgGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTg                                      <  3:304/147‑1 (MQ=255)
                                                                                                    gcgcCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTc                                   <  6:7184/148‑1 (MQ=255)
                                                                                                      gcCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGACCAGGCGGTGGAGACCAGTTCTCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGAGATGTCTTGCTCGCt                                >  2:83010/1‑149 (MQ=255)
                                                                                                         aaGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCAAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGAt                             >  5:22886/1‑149 (MQ=255)
                                                                                                          aGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATc                            <  7:98903/149‑1 (MQ=255)
                                                                                                          aGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATc                            <  2:88571/149‑1 (MQ=255)
                                                                                                          aGCGCAAACATGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATc                            >  6:45502/1‑149 (MQ=255)
                                                                                                               aaaCGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGc                                 <  5:141434/139‑1 (MQ=255)
                                                                                                               aaaCGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGc                                 >  6:141434/1‑139 (MQ=255)
                                                                                                                aaCGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGAt                      <  4:275617/149‑1 (MQ=255)
                                                                                                                  cGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATgg                    <  1:178691/149‑1 (MQ=255)
                                                                                                                        ccGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGc              <  1:85946/149‑1 (MQ=255)
                                                                                                                           ttGGCACCACCGCGCTTGTCGGAGGCGCGGAAGGTCGACGCTGATGCCCAGGCTGTGGAGACCAGTTGGGCCTCGCTCAAACCCGAGGCCTGGATGCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGGCAGCttt           <  4:220686/149‑1 (MQ=255)
                                                                                                                              gCACCACCGCGCTTGTCGGAGTCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGCTTTgg         <  8:160400/148‑1 (MQ=255)
                                                                                                                              gCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGCTTTgg         <  2:158273/148‑1 (MQ=255)
                                                                                                                                 ccaccGCGCTTGTCGGAGTCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCTCGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGCTTTGGGCaa     <  5:172415/149‑1 (MQ=255)
                                                                                                                                    ccGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGCTTTGGGCAACgg  <  2:265540/149‑1 (MQ=255)
                                                                                                                                                gAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTg                                                       <  7:96296/84‑1 (MQ=255)
                                                                                                                                                gAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTg                                                       >  8:96296/1‑84 (MQ=255)
                                                                                                                                                   |                                                                                                                                     
AGGGAAACCTTCTTGCCGCTGCTGTTGAACTCGCCCTGGATCTTCTCCAATGCCGCCAGCACCTTGTCCACGCCCTGGTTTGCAGCCCAGTCCTTCTGTGGCGCCAAGCGCAAACGTGCGCCGTTGGCACCACCGCGCTTGTCGGAGCCGCGGAAGGTCGACGCGGATGCCCAGGCGGTGGAGACCAGTTCGCCCACGCTCAAACCCGAGGCCAGGATCCTGGACTTGAGGGTCGCGATGTCTTGCTCGCTGATCGTGGATGGGTCAGCTTTGGGCAACGG  >  NC_002947/4168244‑4168524

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: