Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,136,042 T→A 5.5% intergenic (+27/‑66) katG → / → yijE hydroperoxidase I/cystine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,136,0420TA 5.5% 705.6 / 20.2 218intergenic (+27/‑66)katG/yijEhydroperoxidase I/cystine exporter
Reads supporting (aligned to +/- strand):  ref base T (89/116);  new base A (5/7);  total (94/124)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.40e-01

Only 100 of 230 total aligned reads displayed.

AACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATTTTC  >  NC_000913/4135893‑4136191
                                                                                                                                                     |                                                                                                                                                     
aaCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgct                                                                                                                                                       <  2:756945/150‑1 (MQ=255)
aaCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgct                                                                                                                                                       >  1:1538961/1‑150 (MQ=255)
  caccGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                     <  2:609079/148‑1 (MQ=255)
  cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCTGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                     <  1:499544/150‑1 (MQ=255)
  cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                     <  1:409190/150‑1 (MQ=255)
  cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                     >  2:1819766/1‑150 (MQ=255)
  cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                     >  1:1624010/1‑150 (MQ=255)
   tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                    <  1:490665/150‑1 (MQ=255)
   tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                    <  1:1453251/150‑1 (MQ=255)
          tgcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGACTGCCAGCAAGCAGCCGCt                                                                                                                                             <  1:1037826/150‑1 (MQ=255)
           gcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCAGTTCAGCGACTGCCAGCAAGCAGCCGCTg                                                                                                                                            >  2:549732/1‑150 (MQ=255)
            cgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                    >  2:412899/1‑141 (MQ=255)
            cgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGa                                                                                                                                           <  1:201870/150‑1 (MQ=255)
              tgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAAc                                                                                                                                         >  1:1657896/1‑150 (MQ=255)
               gcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGACTGCCAGCAAGCAGCCGCTGAACg                                                                                                                                        <  2:1685553/150‑1 (MQ=255)
                cgGTGGCGGAAGTTTACGCCAGTACCGATGCCCACGAGAAGTTTGTTAAAGACTTTGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGt                                                                                                                                       >  2:633146/1‑150 (MQ=255)
                   tGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCt                                                                                                                                    <  2:339118/150‑1 (MQ=255)
                   tGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCt                                                                                                                                    <  2:1976922/150‑1 (MQ=255)
                    ggCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGCGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCtt                                                                                                                                   >  2:775797/1‑150 (MQ=255)
                      cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                                 >  1:1314797/1‑150 (MQ=255)
                      cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                                 <  2:1008547/150‑1 (MQ=255)
                      cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                                 >  1:574606/1‑150 (MQ=255)
                       ggAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTAc                                                                                                                                <  2:382473/150‑1 (MQ=255)
                        gAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTAcc                                                                                                                               <  2:234839/150‑1 (MQ=255)
                         aaGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCa                                                                                                                              <  1:423789/150‑1 (MQ=255)
                           gTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc                                                                                                                            >  2:518674/1‑150 (MQ=255)
                            tttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                                                                           <  1:1782214/150‑1 (MQ=255)
                             ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGt                                                                                                                          <  1:1116454/150‑1 (MQ=255)
                             ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGt                                                                                                                          <  1:69598/150‑1 (MQ=255)
                              tACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGta                                                                                                                         <  1:1101207/150‑1 (MQ=255)
                               aCGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtat                                                                                                                        <  2:1486705/150‑1 (MQ=255)
                                cGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGATAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtata                                                                                                                       <  1:1574621/150‑1 (MQ=255)
                                      tAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGc                                                                                                                 <  1:1745537/150‑1 (MQ=255)
                                      tAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGc                                                                                                                 <  1:1335851/150‑1 (MQ=255)
                                           aTGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGGGAAAGTGATGAACGTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCGGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaacg                                                                                                            <  2:656319/150‑1 (MQ=255)
                                              cccACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaa                                                                                                         <  2:1317795/150‑1 (MQ=255)
                                              cccACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaa                                                                                                         <  2:742348/150‑1 (MQ=255)
                                                 aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                                                                                      <  1:788557/150‑1 (MQ=255)
                                                 aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                                                                                      >  2:28650/1‑150 (MQ=255)
                                                  cGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                                                                                     <  1:1731931/150‑1 (MQ=255)
                                                   gagaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTac                                                                                                    <  1:1290959/150‑1 (MQ=255)
                                                    agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTaca                                                                                                   <  1:379088/150‑1 (MQ=255)
                                                     gaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTacac                                                                                                  <  1:559709/150‑1 (MQ=255)
                                                        gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACt                                                                                               >  2:1428374/1‑150 (MQ=255)
                                                        gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTAGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACt                                                                                               >  1:307833/1‑150 (MQ=255)
                                                         tttGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTg                                                                                              >  2:1725528/1‑150 (MQ=255)
                                                          ttGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTgg                                                                                             >  1:52514/1‑150 (MQ=255)
                                                           tGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGa                                                                                            <  1:986112/150‑1 (MQ=255)
                                                            gTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAt                                                                                           >  1:1772103/1‑150 (MQ=255)
                                                             ttAAAGACTTAGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtc                                                                                          <  2:470265/150‑1 (MQ=255)
                                                               aaaGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctc                                                                                        <  2:353650/150‑1 (MQ=255)
                                                               aaaGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctc                                                                                        <  1:927175/150‑1 (MQ=255)
                                                                 aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTATCATTGTATAGTGGTCGAACGAAATCTACACACTGGAtctctc                                                                                      <  2:1541769/150‑1 (MQ=255)
                                                                 aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctctc                                                                                      <  2:911485/150‑1 (MQ=255)
                                                                  gACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                                                                                     <  2:1706257/150‑1 (MQ=255)
                                                                    cTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATg                                                                                   >  1:956201/1‑150 (MQ=255)
                                                                       cGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCt                                                                                >  2:882051/1‑150 (MQ=255)
                                                                            cggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGc                                                                           <  1:1633907/150‑1 (MQ=255)
                                                                             ggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTATCATTGTATAGTGGTCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCa                                                                          <  2:1574144/150‑1 (MQ=255)
                                                                               cATGGGTGAAAGCGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAgg                                                                        <  1:775797/150‑1 (MQ=255)
                                                                                  gggTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGaaa                                                                     >  1:23148/1‑150 (MQ=255)
                                                                                   ggTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGAGCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAag                                                                    <  2:106464/150‑1 (MQ=255)
                                                                                    gTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAaga                                                                   >  2:585027/1‑150 (MQ=255)
                                                                                      gAAAGTGATGAACCTCGACCGGTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGc                                                                 <  1:1811942/150‑1 (MQ=255)
                                                                                         aGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg                                                                         >  1:1866085/1‑139 (MQ=255)
                                                                                          gTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc                                                             <  2:1314797/150‑1 (MQ=255)
                                                                                               gaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTg                                                        <  2:1713340/150‑1 (MQ=255)
                                                                                                aaCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctc                                                                                        >  2:1594547/1‑117 (MQ=255)
                                                                                                 aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                      <  2:1657896/150‑1 (MQ=255)
                                                                                                 aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                      >  2:1407879/1‑150 (MQ=255)
                                                                                                 aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                      >  2:588804/1‑150 (MQ=255)
                                                                                                  ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAATGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCCTGTCTGCCGCAGGAAAGAGCAACCCACTGGCa                                                     >  2:649116/1‑150 (MQ=255)
                                                                                                   cTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCaa                                                    >  1:1471740/1‑150 (MQ=255)
                                                                                                      gACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCa                                                 >  2:1861409/1‑150 (MQ=255)
                                                                                                      gACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCa                                                 >  1:1884031/1‑150 (MQ=255)
                                                                                                        ccGTTTCCACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc                                                             <  1:1424735/136‑1 (MQ=255)
                                                                                                        ccGTTTCAACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc                                                             >  2:1424735/1‑136 (MQ=255)
                                                                                                         cGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTg                                              <  2:1778225/150‑1 (MQ=255)
                                                                                                         cGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTg                                              <  1:645814/150‑1 (MQ=255)
                                                                                                              ggACCGGCTGTAATCTGACCCCGGTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCGCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                         <  2:1927108/149‑1 (MQ=255)
                                                                                                              cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGTGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                         >  1:334075/1‑150 (MQ=255)
                                                                                                              cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACCCACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                         >  1:282083/1‑150 (MQ=255)
                                                                                                                  ctgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtg                                     <  2:1664540/150‑1 (MQ=255)
                                                                                                                     ctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCACCCCACTGGCAATCAGTGGCCTGGTTGTGc                                  >  1:1558091/1‑150 (MQ=255)
                                                                                                                     ctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGc                                  >  1:1666824/1‑150 (MQ=255)
                                                                                                                      tgTAATCTGACCGCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGGCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCt                                 <  2:1867636/150‑1 (MQ=255)
                                                                                                                      tgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCt                                 <  1:1840874/150‑1 (MQ=255)
                                                                                                                         aaTCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTcac                              >  1:1691868/1‑150 (MQ=255)
                                                                                                                         aaTCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTcac                              <  1:1196613/150‑1 (MQ=255)
                                                                                                                           tCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCacac                            >  2:1288474/1‑150 (MQ=255)
                                                                                                                            cTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACt                           <  2:568602/150‑1 (MQ=255)
                                                                                                                             tGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACtt                          >  2:1086652/1‑150 (MQ=255)
                                                                                                                              aagaacGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTa                         <  1:27739/145‑1 (MQ=255)
                                                                                                                                 cccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCt                      >  1:159190/1‑150 (MQ=255)
                                                                                                                                    gTTCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGa                   >  2:957395/1‑150 (MQ=255)
                                                                                                                                      tCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGt                 <  2:1154796/150‑1 (MQ=255)
                                                                                                                                            gctgctTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGc           <  1:588120/150‑1 (MQ=255)
                                                                                                                                             ctgctTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCt                                          >  2:20692/1‑118 (MQ=255)
                                                                                                                                             ctgctTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCt                                          <  1:20692/118‑1 (MQ=255)
                                                                                                                                                cttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACGTATCTGGAGTTATAGCTGGa       >  2:1960057/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
AACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATTTTC  >  NC_000913/4135893‑4136191

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