Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,134,394 A→G D187G (GAC→GGC)  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,3940AG98.4% 243.1 / NA 64D187G (GAC→GGC) katGcatalase‑peroxidase HPI, heme b‑containing
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base G (42/21);  total (42/22)
Fisher's exact test for biased strand distribution p-value = 3.44e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGG  >  NC_000913/4134246‑4134541
                                                                                                                                                    |                                                                                                                                                   
aGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCg                                                                                                                                                   >  2:397969/1‑151 (MQ=255)
    tCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTg                                                                                                                                               >  2:425517/1‑151 (MQ=255)
                gTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCt                                                                                                                                   >  2:626144/1‑151 (MQ=255)
                  cGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTgg                                                                                                                                 >  1:246834/1‑151 (MQ=255)
                    ccTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGAt                                                                                                                               <  2:897596/151‑1 (MQ=255)
                     cTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATg                                                                                                                              >  1:712004/1‑151 (MQ=255)
                       gttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATg                                                                                                                              <  1:257473/149‑1 (MQ=255)
                       gttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGtt                                                                                                                            <  2:674827/151‑1 (MQ=255)
                       gttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAAACGGATCTGGATGtt                                                                                                                            >  2:789784/1‑151 (MQ=255)
                        ttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGtt                                                                                                                            >  2:187933/1‑150 (MQ=255)
                        ttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGtt                                                                                                                            >  1:985102/1‑150 (MQ=255)
                        ttgtGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGtt                                                                                                                            >  2:530173/1‑150 (MQ=255)
                          gtgGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAAc                                                                                                                         >  2:127155/1‑151 (MQ=255)
                            ggCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTg                                                                                                                       >  1:6141/1‑151 (MQ=255)
                              ccAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTggg                                                                                                                     >  2:409848/1‑151 (MQ=255)
                                    aaaCAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGtgatg                                                                                                               >  2:310095/1‑151 (MQ=255)
                                       cAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGaaa                                                                                                            <  2:448773/151‑1 (MQ=255)
                                        aGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGaaa                                                                                                            >  2:111430/1‑150 (MQ=255)
                                          aaaTATGGTCAGAAAATCTCCTGTGCCGACCTGTTTATCCTCGCGGGTAAGGTGGCGTTAGAAATCTCCGTTTTTATTTCTTACGGCTCTGCTGCCCGTCGTGAAGACGGCGGGGAACCGGATCTTTCCCTTAACTGGGCTCATGcaaatc                                                                                                         <  2:693837/151‑3 (MQ=255)
                                           aaTATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGc                                                                                                        <  1:641217/151‑1 (MQ=255)
                                             tatGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCt                                                                                                      >  1:896669/1‑151 (MQ=255)
                                             tatGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCt                                                                                                      >  1:869216/1‑151 (MQ=255)
                                              atGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTg                                                                                                     >  1:598162/1‑151 (MQ=255)
                                                ggTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTgg                                                                                                    <  2:856128/150‑1 (MQ=255)
                                                ggTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTgg                                                                                                    <  2:691516/150‑1 (MQ=255)
                                                 gTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCt                                                                                                  <  2:524541/151‑1 (MQ=255)
                                                          tctcCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCg                                                                                         >  1:315119/1‑151 (MQ=255)
                                                            tcCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTc                                                                                       >  2:423245/1‑151 (MQ=255)
                                                            tcCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTc                                                                                       >  1:347388/1‑151 (MQ=255)
                                                             ccTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTGCCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCg                                                                                         <  2:316913/148‑1 (MQ=255)
                                                                     gACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGa                                                                                >  1:509106/1‑149 (MQ=255)
                                                                      aCCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAgc                                                                             <  1:71319/151‑1 (MQ=255)
                                                                          gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCCCCGTCATCCGGAAGCGCTg                                                                         <  2:679010/151‑1 (MQ=255)
                                                                          gTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCt                                                                          <  2:930319/150‑1 (MQ=255)
                                                                             tATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCg                                                                      >  1:383088/1‑151 (MQ=255)
                                                                              aTCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCg                                                                      <  2:772372/150‑1 (MQ=255)
                                                                                ccTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGaaa                                                                   <  2:253471/151‑1 (MQ=255)
                                                                                ccTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGaaa                                                                   <  1:856592/151‑1 (MQ=255)
                                                                                ccTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGaa                                                                    <  2:167419/150‑1 (MQ=255)
                                                                                     cgGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCAcc                                                              >  1:499411/1‑151 (MQ=255)
                                                                                           aCGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTggg                                                        >  1:79562/1‑151 (MQ=255)
                                                                                            cGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTggg                                                        >  1:140392/1‑150 (MQ=255)
                                                                                             gTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGt                                                       <  1:173433/150‑1 (MQ=255)
                                                                                                    tAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCg                                                >  2:810420/1‑150 (MQ=255)
                                                                                                              cGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTg                                     >  1:307736/1‑151 (MQ=255)
                                                                                                                gCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAt                                   >  2:364795/1‑151 (MQ=255)
                                                                                                                 cTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAtt                                  >  1:480882/1‑151 (MQ=255)
                                                                                                                   tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAttt                                 >  2:178075/1‑150 (MQ=255)
                                                                                                                   tCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAttt                                 >  1:39824/1‑150 (MQ=255)
                                                                                                                    ccGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTAc                               <  2:428254/151‑1 (MQ=255)
                                                                                                                     cGTACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGAttt                                 <  1:589654/148‑1 (MQ=255)
                                                                                                                       tACCTTCGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGtt                            >  2:71078/1‑151 (MQ=255)
                                                                                                                             cGGTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCg                      >  2:131840/1‑151 (MQ=255)
                                                                                                                              ggTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCgg                     >  2:640373/1‑151 (MQ=255)
                                                                                                                              ggTTTTGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCgg                     >  1:688583/1‑151 (MQ=255)
                                                                                                                                  ttGGTGCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAg                  <  1:573142/150‑1 (MQ=255)
                                                                                                                                       gCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCgg            >  1:239188/1‑151 (MQ=255)
                                                                                                                                       gCCGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCg             <  1:385193/150‑1 (MQ=255)
                                                                                                                                        ccGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCgg            >  1:587172/1‑150 (MQ=255)
                                                                                                                                         cGGTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGAt          >  1:678373/1‑151 (MQ=255)
                                                                                                                                          ggTCGTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATc         >  2:428605/1‑151 (MQ=255)
                                                                                                                                              gTGAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATcac       >  1:204010/1‑149 (MQ=255)
                                                                                                                                                gAAGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcg   >  2:311155/1‑151 (MQ=255)
                                                                                                                                                 aaGGCGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAgcgg  >  1:900331/1‑151 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AGCCTCGATAAAGCGCGTCGCCTGTTGTGGCCAATCAAACAGAAATATGGTCAGAAAATCTCCTGGGCCGACCTGTTTATCCTCGCGGGTAACGTGGCGCTAGAAAACTCCGGCTTCCGTACCTTCGGTTTTGGTGCCGGTCGTGAAGACGTCTGGGAACCGGATCTGGATGTTAACTGGGGTGATGAAAAAGCCTGGCTGACTCACCGTCATCCGGAAGCGCTGGCGAAAGCACCGCTGGGTGCAACCGAGATGGGTCTGATTTACGTTAACCCGGAAGGCCCGGATCACAGCGG  >  NC_000913/4134246‑4134541

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