Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,134,511 A→T 5.4% Y226F (TAC→TTC)  katG → catalase‑peroxidase HPI, heme b‑containing

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,134,5110AT 5.4% 240.7 / 3.5 74Y226F (TAC→TTC) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base A (38/32);  new base T (2/2);  total (40/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 = 9.27e-01

GTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTAAAACCCACGG  >  NC_000913/4134364‑4134655
                                                                                                                                                   |                                                                                                                                                
gTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTAc                                                                                                                                                 >  3:132438/1‑149 (MQ=255)
gTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTAc                                                                                                                                                 <  6:276575/149‑1 (MQ=255)
gTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTAACCGTCATACGTAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCAGATTTAc                                                                                                                                                 >  1:216645/1‑149 (MQ=255)
        ggTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAAcc                                                                                                                                          <  4:316129/148‑1 (MQ=255)
                   cGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCg                                                                                                                               >  5:126375/1‑148 (MQ=255)
                                gTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGTGATGGGTCTGATTTACGTTAACCCGGAAGGCCCgg                                                                                                                              <  1:69935/136‑1 (MQ=255)
                                gTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGTGATGGGTCTGATTTACGTTAACCCGGAAGGCCCgg                                                                                                                              >  2:69935/1‑136 (MQ=255)
                                gTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcggcg                                                                                                                  >  7:56231/1‑148 (MQ=255)
                                   tGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCTCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAAcc                                                                                                              >  5:76024/1‑149 (MQ=255)
                                    gggAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCg                                                                                                             >  4:31101/1‑149 (MQ=255)
                                     ggAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGc                                                                                                            >  4:184588/1‑149 (MQ=255)
                                      gAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGt                                                                                                                                               >  2:87956/1‑113 (MQ=255)
                                      gAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGt                                                                                                                                               <  1:87956/113‑1 (MQ=255)
                                           ggATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCTAACCGCTTTCTg                                                                                                      <  6:164473/149‑1 (MQ=255)
                                            gATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGc                                                                                                     <  4:8215/149‑1 (MQ=255)
                                                tGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCgg                                                                                                   >  3:1415/1‑147 (MQ=255)
                                                  gATGTTCACTGGGGTGCTGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTACCCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCgg                                                                                                   <  8:122715/145‑1 (MQ=255)
                                                  gATGTTAACTGGGGTGATGAACAAGCCTGTCTGACTCACCGTCATCCGTAATCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                                                                                <  4:294850/148‑1 (MQ=255)
                                                  gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGAAGGGTCTGATTTTCGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagc                                                                                               <  4:300596/149‑1 (MQ=255)
                                                  gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGAAGGGTCTGATTTTCGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcagc                                                                                               <  4:305403/149‑1 (MQ=255)
                                                  gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGAAGGGTCTGATTTTCGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                                                                                >  3:300596/1‑148 (MQ=255)
                                                  gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGAAGGGTCTGATTTTCGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGgcag                                                                                                >  3:305403/1‑148 (MQ=255)
                                                  gATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGGGCAACCGAGATGGGGCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCg                                                                                                    <  7:171989/144‑1 (MQ=255)
                                                         aCTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCAGAAGGCCCTGATCACAGCGGCGAACCGCTTTCTGCGGCAGAAGCTAt                                                                                         <  1:107342/148‑1 (MQ=255)
                                                          cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTg                                                                                                      >  5:136384/1‑134 (MQ=255)
                                                          cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTg                                                                                                      <  6:136384/134‑1 (MQ=255)
                                                          cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                                                                         >  6:81439/1‑147 (MQ=255)
                                                          cTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTAt                                                                                         <  5:81439/147‑1 (MQ=255)
                                                             gggTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGGCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCcg                                                                                      <  1:284455/147‑1 (MQ=255)
                                                                 gatgaAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGAGGCAGCAGCTATCCGCGCGAc                                                                                >  3:274002/1‑149 (MQ=255)
                                                                   tgaAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCt                                                                              >  5:339282/1‑149 (MQ=255)
                                                                     aaaaaGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCtt                                                                             <  7:953/148‑1 (MQ=255)
                                                                     aaaaaGCCTGGCGGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCTAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCTAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCtt                                                                             <  5:84606/148‑1 (MQ=255)
                                                                      aaaaGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGCCCttgg                                                                           >  5:135115/1‑147 (MQ=255)
                                                                           ccTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCaa                                                                       >  1:92096/1‑148 (MQ=255)
                                                                            cTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa                                                                     <  4:51626/149‑1 (MQ=255)
                                                                            cTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACa                                                                     <  3:184588/149‑1 (MQ=255)
                                                                               gCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATgg                                                                  >  3:88737/1‑149 (MQ=255)
                                                                                cTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCATCAGCTATCCGCGAGACCTTCGGCAACATggg                                                                 <  2:216645/149‑1 (MQ=255)
                                                                                     tCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATg                                                             >  4:1960/1‑148 (MQ=255)
                                                                                     tCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATg                                                             >  2:267528/1‑148 (MQ=255)
                                                                                       aCCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac                                                          >  3:148049/1‑149 (MQ=255)
                                                                                       aCCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac                                                          <  3:1960/149‑1 (MQ=255)
                                                                                       aCCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCAGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAac                                                          >  8:221607/1‑149 (MQ=255)
                                                                                         cGTCATCCGGAAGCGCTGGCGAAACCACCTCTGGGTGCACCACAGAAGGGTCAGAATAACGTGAACCCTGAAGGCCCTGATCACACCGCACAACCGCTTTCTGCTGCAGCAGCTATCCTCGCGACCATCGGCAACAAGGGCATGCacga                                                        >  4:157410/1‑149 (MQ=255)
                                                                                           tCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACTGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacgac                                                       >  6:135873/1‑148 (MQ=255)
                                                                                             aTCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACgaa                                                    <  8:56231/149‑1 (MQ=255)
                                                                                              tCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAacgacg                                                      >  8:91561/1‑146 (MQ=255)
                                                                                                cGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACgaaga                                                  >  1:344825/1‑148 (MQ=255)
                                                                                                cGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACgaaga                                                  <  6:339282/148‑1 (MQ=255)
                                                                                                   aaGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAAcc                                              >  8:68127/1‑149 (MQ=255)
                                                                                                      cgcTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTg                                           <  3:200332/149‑1 (MQ=255)
                                                                                                       gcTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTgg                                          >  3:132817/1‑149 (MQ=255)
                                                                                                             gAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGa                                    <  1:267528/149‑1 (MQ=255)
                                                                                                             gAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGa                                    >  6:189802/1‑149 (MQ=255)
                                                                                                                      gCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGggtg                           >  3:305427/1‑149 (MQ=255)
                                                                                                                      gCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGggtg                           <  4:59381/149‑1 (MQ=255)
                                                                                                                       cTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGggtgg                          <  4:132817/149‑1 (MQ=255)
                                                                                                                          ggTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCa                       >  8:285302/1‑149 (MQ=255)
                                                                                                                          ggTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCa                       >  8:143575/1‑149 (MQ=255)
                                                                                                                                  cGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTgg               <  2:19300/149‑1 (MQ=255)
                                                                                                                                   gagaTGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACCGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTggg              >  5:149904/1‑149 (MQ=255)
                                                                                                                                   gagaTGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTggg              <  4:148049/149‑1 (MQ=255)
                                                                                                                                     gaTGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGt             >  3:105375/1‑148 (MQ=255)
                                                                                                                                      aTGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTg                                           >  7:332514/1‑117 (MQ=255)
                                                                                                                                      aTGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTg                                           <  8:332514/117‑1 (MQ=255)
                                                                                                                                         ggTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTaaaa         <  3:129516/148‑1 (MQ=255)
                                                                                                                                         ggTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTAAAAc        <  4:88737/149‑1 (MQ=255)
                                                                                                                                           tCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGACCGACGAAGAACCCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTAAAAcc       >  5:195631/1‑148 (MQ=255)
                                                                                                                                            cTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTggg              >  3:180938/1‑140 (MQ=255)
                                                                                                                                            cTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTggg              <  4:180938/140‑1 (MQ=255)
                                                                                                                                             tGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTAAAACCCAc    >  5:346067/1‑149 (MQ=255)
                                                                                                                                               aTTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCTACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGCGTGGTCATACGCTGGGTAAAAccaccgg  >  2:315911/1‑144 (MQ=255)
                                                                                                                                               aTTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTAAAACCCAc    <  7:68127/147‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
GTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGGCGAACCGCTTTCTGCGGCAGCAGCTATCCGCGCGACCTTCGGCAACATGGGCATGAACGACGAAGAAACCGTGGCGCTGATTGCGGGTGGTCATACGCTGGGTAAAACCCACGG  >  NC_000913/4134364‑4134655

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: