Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,136,037 C→T intergenic (+22/‑71) katG → / → yijE catalase‑peroxidase HPI, heme b‑containing/EamA‑like transporter family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,136,0370CT97.8% 599.1 / 0.9 179intergenic (+22/‑71)katG/yijEcatalase‑peroxidase HPI, heme b‑containing/EamA‑like transporter family protein
Reads supporting (aligned to +/- strand):  ref base C (2/2);  new base T (68/106);  total (71/108)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.72e-01

Only 100 of 191 total aligned reads displayed.

GCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCC  >  NC_000913/4135938‑4136133
                                                                                                   |                                                                                                
  ccACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgttt                                                                                                <  1:1972574/100‑1 (MQ=255)
   cACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgtttg                                                                                               <  1:3452095/100‑1 (MQ=255)
   cACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgtttg                                                                                               >  2:2414629/1‑100 (MQ=255)
    aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgtttgc                                                                                              >  1:921322/1‑100 (MQ=255)
    aCGAGAAGTGTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgtttgc                                                                                              <  2:1474687/100‑1 (MQ=255)
      gagaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTg                                                                                            <  2:1629576/100‑1 (MQ=255)
      gagaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTg                                                                                            <  1:3075680/100‑1 (MQ=255)
      gagaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTg                                                                                            <  1:2967479/100‑1 (MQ=255)
       agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTgg                                                                                           <  2:3015248/100‑1 (MQ=255)
       agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTgg                                                                                           <  1:693382/100‑1 (MQ=255)
       agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTgg                                                                                           <  1:2799999/100‑1 (MQ=255)
        gaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGc                                                                                          <  1:1319543/100‑1 (MQ=255)
        gaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGc                                                                                          <  1:2793201/100‑1 (MQ=255)
           gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGt                                                                                       <  2:797167/100‑1 (MQ=255)
           gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGt                                                                                       >  2:3716117/1‑100 (MQ=255)
           gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGt                                                                                       >  1:2709686/1‑100 (MQ=255)
             ttGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCg                                                                                     <  2:1417090/100‑1 (MQ=255)
             ttGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCg                                                                                     >  2:2437878/1‑100 (MQ=255)
                ttAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTg                                                                                  <  2:1675930/100‑1 (MQ=255)
                 tAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTGATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGa                                                                                 >  1:2016713/1‑100 (MQ=255)
                   aaGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAAc                                                                               <  2:2201513/100‑1 (MQ=255)
                    aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACg                                                                              <  2:1282381/100‑1 (MQ=255)
                    aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACg                                                                              <  2:277590/100‑1 (MQ=255)
                    aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACg                                                                              <  2:1706728/100‑1 (MQ=255)
                    aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACg                                                                              <  2:662698/100‑1 (MQ=255)
                       cTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTc                                                                           >  1:4726694/1‑100 (MQ=255)
                       cTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTc                                                                           >  1:1967435/1‑100 (MQ=255)
                       cTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTc                                                                           <  2:2011244/100‑1 (MQ=255)
                         tCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCtt                                                                         <  1:619909/100‑1 (MQ=255)
                          cGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCttt                                                                        <  2:2875022/100‑1 (MQ=255)
                          cGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCttt                                                                        <  2:3785557/100‑1 (MQ=255)
                          cGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCttt                                                                        <  1:2505878/100‑1 (MQ=255)
                           gTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                       >  2:3454696/1‑100 (MQ=255)
                            tGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTAc                                                                      >  2:4108576/1‑100 (MQ=255)
                             ggcggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTAcc                                                                     <  1:2451621/100‑1 (MQ=255)
                             ggcggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTAcc                                                                     <  1:1090230/100‑1 (MQ=255)
                             ggcggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTAcc                                                                     <  1:4533066/100‑1 (MQ=255)
                                 gcATGGGTGAAAGTGATGAGCGTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                 >  1:4395531/1‑100 (MQ=255)
                                 gcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                 <  1:2266124/100‑1 (MQ=255)
                                 gcAGGGGTGAAAGTGATGAGCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                 <  2:3892841/100‑1 (MQ=255)
                                  cATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGACCGTTCTTTACCAGCGt                                                                >  1:5151222/1‑100 (MQ=255)
                                   aTGGGTGAAAGTGATTAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGta                                                               <  1:5002799/100‑1 (MQ=255)
                                    tGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtat                                                              >  1:1115166/1‑100 (MQ=255)
                                    tGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtat                                                              <  2:1551365/100‑1 (MQ=255)
                                    tGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtat                                                              <  2:4134284/100‑1 (MQ=255)
                                       gTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGt                                                           <  1:4564593/100‑1 (MQ=255)
                                         gAAAGTGCTGCACCTCGAACGTTTCGACCTGCTGTAATCAGACCCAGTTCCGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTgg                                                         <  1:2335475/100‑1 (MQ=255)
                                          aacgTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTggg                                                        >  2:430483/4‑100 (MQ=255)
                                             gTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcga                                                     <  2:3254030/100‑1 (MQ=255)
                                               gatgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaac                                                   <  1:62475/100‑1 (MQ=255)
                                                   aaCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaa                                               >  1:3918744/1‑100 (MQ=255)
                                                    tccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaaa                                              <  2:3828221/99‑1 (MQ=255)
                                                    aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaaa                                              >  1:352998/1‑100 (MQ=255)
                                                    aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaaa                                              >  2:1583849/1‑100 (MQ=255)
                                                    aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaaa                                              >  1:3722544/1‑100 (MQ=255)
                                                     ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAAc                                             >  2:4346603/1‑100 (MQ=255)
                                                     ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAAc                                             >  1:3897364/1‑100 (MQ=255)
                                                     ccTCGACCGTTTCGACCTGCTGTAATCGGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTAACAGCGTTTAGTGGTCGATCGAAAAc                                             <  2:1047910/100‑1 (MQ=255)
                                                      cTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                            >  1:868202/1‑100 (MQ=255)
                                                       tCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                           <  2:739559/100‑1 (MQ=255)
                                                       tCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                           >  2:225850/1‑100 (MQ=255)
                                                         gACCGTTTCGACCTGCTGTAATCTTACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTacc                                         <  2:3230403/100‑2 (MQ=255)
                                                         gACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTaca                                         <  1:4271012/100‑1 (MQ=255)
                                                          aCCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTacac                                        <  2:1929381/100‑1 (MQ=255)
                                                          aCCGTTTCGACCTGCTGTAATCTGACCCCCTGCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTacac                                        <  1:2829149/100‑1 (MQ=255)
                                                             gTTTCGACCTGCTGTGATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACt                                     <  2:1991617/100‑1 (MQ=255)
                                                             gTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACt                                     <  1:2610980/100‑1 (MQ=255)
                                                              tttCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTg                                    >  1:1056040/1‑100 (MQ=255)
                                                              tttCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTg                                    <  2:1559893/100‑1 (MQ=255)
                                                               ttCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTgg                                   <  2:3611679/100‑1 (MQ=255)
                                                               ttCGACCTGCGGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTgg                                   >  2:4674671/1‑100 (MQ=255)
                                                                tCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGa                                  <  1:282474/100‑1 (MQ=255)
                                                                tCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGa                                  >  1:1660638/1‑100 (MQ=255)
                                                                 cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATCTTGGGCGAACGCAAACTACACACTGGAt                                 >  2:1213643/1‑100 (MQ=255)
                                                                  gACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtc                                <  2:1875860/100‑1 (MQ=255)
                                                                  gACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtc                                <  2:4406523/100‑1 (MQ=255)
                                                                  gACCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtc                                <  2:313628/100‑1 (MQ=255)
                                                                   aCCTGCTGTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtct                               <  1:1606014/100‑1 (MQ=255)
                                                                     ctgctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctct                             <  1:3686131/100‑1 (MQ=255)
                                                                      tgctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctctc                            <  2:3142359/100‑1 (MQ=255)
                                                                       gctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                           <  2:450222/100‑1 (MQ=255)
                                                                       gctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                           >  2:426055/1‑100 (MQ=255)
                                                                       gctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                           <  1:4954001/100‑1 (MQ=255)
                                                                        ctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGACCGAAAACTACACACTGGATCTCTCAt                          >  1:3605154/1‑100 (MQ=255)
                                                                        ctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCAt                          <  2:2602690/100‑1 (MQ=255)
                                                                        ctgTAATCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATACTGGGCGAACGAAAACTACACACTGGATCCCTCAt                          >  2:92204/1‑100 (MQ=255)
                                                                             aTCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTg                     >  1:1596101/1‑100 (MQ=255)
                                                                              tCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACTAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGc                    <  2:4878836/100‑1 (MQ=255)
                                                                              tCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGc                    <  1:1880536/100‑1 (MQ=255)
                                                                              tCTGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGc                    >  1:1028062/1‑100 (MQ=255)
                                                                                tGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCg                  <  1:3838629/100‑1 (MQ=255)
                                                                                tGACCCCGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCg                  >  1:1136624/1‑100 (MQ=255)
                                                                                   ccccGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg               <  1:4641903/100‑1 (MQ=255)
                                                                                   ccccGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg               >  1:974233/1‑100 (MQ=255)
                                                                                   ccccGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg               >  2:4820760/1‑100 (MQ=255)
                                                                                     ccGTTCAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGa             >  1:2054120/1‑100 (MQ=255)
                                                                                          cAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAagag        >  1:3394067/1‑100 (MQ=255)
                                                                                          cAGCGGCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAagag        <  2:2633371/100‑1 (MQ=255)
                                                                                              ggCTGTTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAc    <  2:3583744/100‑1 (MQ=255)
                                                                                               gctgttTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc   <  2:1942718/100‑1 (MQ=255)
                                                                                                   |                                                                                                
GCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCC  >  NC_000913/4135938‑4136133

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