Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,136,043 G→A 5.5% intergenic (+28/‑65) 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,0430GA 5.5% 406.0 / 7.3 126intergenic (+28/‑65)katG/yijEcatalase‑peroxidase HPI, heme b‑containing/EamA‑like transporter family protein
Reads supporting (aligned to +/- strand):  ref base G (61/58);  new base A (4/3);  total (65/61)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.64e-01

Only 100 of 130 total aligned reads displayed.

ACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATTTTCA  >  NC_000913/4135894‑4136192
                                                                                                                                                     |                                                                                                                                                     
aCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCGCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTg                                                                                                                                                       <  2:1449590/150‑1 (MQ=255)
  tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                     <  1:1632526/150‑1 (MQ=255)
     gTCCTGCGTGCAGTGGCGGTAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTCATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGt                                                                                                                                                  <  2:265892/150‑1 (MQ=255)
       ccTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCg                                                                                                                                                <  1:1231573/150‑1 (MQ=255)
        cTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGc                                                                                                                                               >  1:99643/1‑150 (MQ=255)
         tgcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCt                                                                                                                                              >  2:365881/1‑150 (MQ=255)
          gcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACTAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGACTGCCAGCAAGCAGCCGCTg                                                                                                                                             <  2:636391/150‑1 (MQ=255)
          gcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTg                                                                                                                                             <  2:1149772/150‑1 (MQ=255)
          gcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGACTGCCAGCAAGCAGCCGCTg                                                                                                                                             <  1:1119555/150‑1 (MQ=255)
            gtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGaa                                                                                                                                           >  1:721166/1‑150 (MQ=255)
            gtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGaa                                                                                                                                           <  1:1039646/150‑1 (MQ=255)
                  tGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTAGTTAAAGACTTCGTCGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCt                                                                                                                                     >  2:37542/1‑150 (MQ=255)
                        aaGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCa                                                                                                                               <  1:1163712/150‑1 (MQ=255)
                         aGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGGCGCTGAACGTTCTTTACCAg                                                                                                                              <  1:1104732/150‑1 (MQ=255)
                           tttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                                                                            <  2:940013/150‑1 (MQ=255)
                             tACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGta                                                                                                                          >  1:1068069/1‑150 (MQ=255)
                               cGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtata                                                                                                                        >  1:1191900/1‑150 (MQ=255)
                                gCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAg                                                                                                                       >  2:1629718/1‑150 (MQ=255)
                                 ccAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGGGGCGGCATGGGTGAAAGTGATGAACCGCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGt                                                                                                                      <  2:796319/150‑1 (MQ=255)
                                   aGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTgg                                                                                                                    >  2:1093881/1‑150 (MQ=255)
                                     tAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGc                                                                                                                  >  2:1531155/1‑150 (MQ=255)
                                          aTGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaacg                                                                                                             <  1:1231164/150‑1 (MQ=255)
                                            gCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaacgaa                                                                                                           >  1:1489753/1‑150 (MQ=255)
                                              ccACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCAACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaaa                                                                                                         >  2:1705131/1‑150 (MQ=255)
                                               cACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAAc                                                                                                        >  2:966529/1‑150 (MQ=255)
                                                aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                                                                                       >  2:1193048/1‑150 (MQ=255)
                                                aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                                                                                       >  1:163681/1‑150 (MQ=255)
                                                 cGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                                                                                      <  2:555504/150‑1 (MQ=255)
                                                  gagaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTac                                                                                                     <  1:981569/150‑1 (MQ=255)
                                                   agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTaca                                                                                                    <  1:1584997/150‑1 (MQ=255)
                                                     aaGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTacaca                                                                                                  <  1:365881/150‑1 (MQ=255)
                                                      aGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTacacac                                                                                                 >  1:685889/1‑150 (MQ=255)
                                                       gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACt                                                                                                >  2:282867/1‑150 (MQ=255)
                                                          tGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGa                                                                                             <  2:721166/150‑1 (MQ=255)
                                                           gTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTTCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAt                                                                                            >  1:709452/1‑150 (MQ=255)
                                                            ttAACGACTTCGTGGCGGCATGGGTGAAAGTGAAGAACCTCGACCGTTTCGACCTGCTGTACTCTGCCCCCGTTCAGCGGCTGCTTGCTGGCAGTCTCTGAACTTTCTTTACCATCGTATAGTGGTCGAACGAAAACTACACACTGGAtc                                                                                           <  2:571680/150‑1 (MQ=255)
                                                              aaaGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctc                                                                                         >  2:1644580/1‑150 (MQ=255)
                                                               aaGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctct                                                                                        <  2:112542/150‑1 (MQ=255)
                                                                 gACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                                                                                      >  2:28552/1‑150 (MQ=255)
                                                                 gACTTCGTGGCGGCATGGGTGAAAGTGATGAACATCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                                                                                      >  2:1234154/1‑150 (MQ=255)
                                                                    ttCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGt                                                                                   <  1:921282/150‑1 (MQ=255)
                                                                         ggcggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGcc                                                                              <  1:1204042/150‑1 (MQ=255)
                                                                              cATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAgg                                                                         <  1:673746/150‑1 (MQ=255)
                                                                               aTGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGa                                                                        <  2:1392341/150‑1 (MQ=255)
                                                                               aTGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGa                                                                        <  1:1448132/150‑1 (MQ=255)
                                                                               aTGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGa                                                                        <  2:1435673/150‑1 (MQ=255)
                                                                                tGGGTGAAAGTGATGCACCTCGCCCGTTTCGACCTGCTGTATTCTGCCCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACGGGATCTCTCATGTCTGCCGCCGGaa                                                                       >  2:684684/1‑150 (MQ=255)
                                                                                tGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCATCTGCTGCTTGCTGGCAGTCGCTGAACGTTATTTACCAGCGTAAAGTGGGCGAACGAAAACAACACACTCGATCTCTCATGTCAGCCGCAGGaa                                                                       >  2:607802/1‑150 (MQ=255)
                                                                                 gggTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGaaa                                                                      >  2:50785/1‑150 (MQ=255)
                                                                                    tGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAagag                                                                   <  1:1331531/150‑1 (MQ=255)
                                                                                    tGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAagag                                                                   <  1:1028469/150‑1 (MQ=255)
                                                                                     gAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGc                                                                  <  2:163681/150‑1 (MQ=255)
                                                                                       aaGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAATAGCaa                                                                >  1:1318151/1‑150 (MQ=255)
                                                                                         gTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc                                                              <  1:12373/150‑1 (MQ=255)
                                                                                          tgatgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGGATAGTGGGCGAACGAAAACTACACACTTGATCTCTCATGTCTGACGCAGGAAAGAGCCAccc                                                             >  2:203697/1‑150 (MQ=255)
                                                                                              gaACCTCGACCGTTTCGACCTTCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTCTAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTg                                                         <  1:250932/150‑1 (MQ=255)
                                                                                                aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                       >  1:666178/1‑150 (MQ=255)
                                                                                                aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGATGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAACCTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                       <  1:894683/150‑1 (MQ=255)
                                                                                                 ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCa                                                      >  2:899797/1‑150 (MQ=255)
                                                                                                 ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCa                                                      >  1:638289/1‑150 (MQ=255)
                                                                                                 ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCa                                                      >  1:868547/1‑150 (MQ=255)
                                                                                                    cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGCACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCGGCCGCAGGAAAGAGCAACCCACGGGCCATc                                                   >  1:1445185/1‑150 (MQ=255)
                                                                                                    cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATc                                                   >  1:407020/1‑150 (MQ=255)
                                                                                                    cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATc                                                   <  1:966529/150‑1 (MQ=255)
                                                                                                             cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                          >  1:1379528/1‑150 (MQ=255)
                                                                                                             cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                          <  1:459095/150‑1 (MQ=255)
                                                                                                             cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                          >  2:1478451/1‑150 (MQ=255)
                                                                                                              gACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTgg                                         <  1:1005457/150‑1 (MQ=255)
                                                                                                                  tgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATTGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtgt                                     >  2:1670694/1‑150 (MQ=255)
                                                                                                                  tgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATGGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtgt                                     >  2:673550/1‑150 (MQ=255)
                                                                                                                  tgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtgt                                     <  2:868547/150‑1 (MQ=255)
                                                                                                                     tgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCt                                  <  1:8578/150‑1 (MQ=255)
                                                                                                                        aaTCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTcac                               <  1:953895/150‑1 (MQ=255)
                                                                                                                            tGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTAGCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCTTGTCTGCCGCGGGAAAGAGCAACCCACTGGCAATCCGTGGCCTGGTTGTGCTCACACtt                           >  1:296049/1‑150 (MQ=255)
                                                                                                                             gACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTGACACttt                          >  1:509655/1‑149 (MQ=255)
                                                                                                                               ccccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATc                        >  2:676969/1‑150 (MQ=255)
                                                                                                                                cccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCt                       >  2:1404381/1‑150 (MQ=255)
                                                                                                                                 acGTTCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTg                      >  2:453953/2‑150 (MQ=255)
                                                                                                                                  cGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTgg                     <  2:1670454/150‑1 (MQ=255)
                                                                                                                                    ttCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAg                   <  1:282867/150‑1 (MQ=255)
                                                                                                                                     tCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGt                  >  1:908547/1‑150 (MQ=255)
                                                                                                                                     tCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGt                  <  1:769938/150‑1 (MQ=255)
                                                                                                                                     tCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGt                  >  1:449845/1‑150 (MQ=255)
                                                                                                                                     tCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGt                  >  2:1356230/1‑150 (MQ=255)
                                                                                                                                      cAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGtt                 <  1:1629718/150‑1 (MQ=255)
                                                                                                                                      cAGCGGCTGCTTGCTGGCAGTCGCTGAACGGTCAGTACCAGCGTATGGTGGGCGGGCGAAAACTAGCCACTGGCTCTCTCATGTCTGCCGCAGGAAAGACCAACCCACTGGCAAACAGTGGCCTGGTTGTGCTCACACTTATCTGGAggt                 >  1:253362/1‑148 (MQ=255)
                                                                                                                                       aGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTTGTTGTGCTCACACTTATCTGGAGTta                >  1:1114707/1‑150 (MQ=255)
                                                                                                                                       aGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTta                >  1:43023/1‑150 (MQ=255)
                                                                                                                                        gCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTtat               >  1:535794/1‑150 (MQ=255)
                                                                                                                                             tgcttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGTTCTCTCATGTCTGCCGCAGGAAAGAGCAACCCCCTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTg          <  2:722637/150‑1 (MQ=255)
                                                                                                                                             tgcttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTg          >  2:791779/1‑150 (MQ=255)
                                                                                                                                              gctttctGGCAGTCGCTGAAGGTTGTTTACAAGTGTTTAGTGGTCGAACGAAATCTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTgg         <  2:1318151/150‑1 (MQ=255)
                                                                                                                                               cttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGa        <  1:278079/150‑1 (MQ=255)
                                                                                                                                               cttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGCAAAGAGCAACCCACTAGCAATCAGTGGCCTGGTTGTGCTCACACTAATCTGGAGTAAAAGCgggt        >  1:240794/1‑149 (MQ=255)
                                                                                                                                                 tgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTgg         <  2:599901/147‑1 (MQ=255)
                                                                                                                                                 tgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTgg         >  1:599901/1‑147 (MQ=255)
                                                                                                                                                   ctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTAGACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGAtttt    >  2:888508/1‑150 (MQ=255)
                                                                                                                                                   ctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGCCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGAtttt    >  2:503731/1‑150 (MQ=255)
                                                                                                                                                   ctGGCAGTCGCTGAACGTTATTTACCAGCGTATAGTGGGCGAACGAAAACTACACGCTGGCTCTCGCTTGTCTGCCGCAGGAAAGCGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATATGGAGTTATAGATGGAtttt    >  2:97546/1‑150 (MQ=255)
                                                                                                                                                     ggCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATTTTCa  >  1:1144566/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATTTTCA  >  NC_000913/4135894‑4136192

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