Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,136,039 T→A 5.5% intergenic (+24/‑69) 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,0390TA 5.5% 726.5 / 19.6 218intergenic (+24/‑69)katG/yijEhydroperoxidase I/cystine exporter
Reads supporting (aligned to +/- strand):  ref base T (89/117);  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 = 9.95e-01

Only 100 of 232 total aligned reads displayed.

TCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATT  >  NC_000913/4135890‑4136188
                                                                                                                                                     |                                                                                                                                                     
tCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgctt                                                                                                                                                       >  2:1453251/1‑150 (MQ=255)
  tAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgc                                                                                                                                                     >  1:1005549/1‑150 (MQ=255)
  tAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgc                                                                                                                                                     >  2:1290959/1‑150 (MQ=255)
  tAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGACGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgc                                                                                                                                                     >  1:746542/1‑150 (MQ=255)
   aaCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgct                                                                                                                                                    >  2:419866/1‑150 (MQ=255)
    aCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTg                                                                                                                                                   <  2:39681/150‑1 (MQ=255)
     cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCTGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                  <  1:499544/150‑1 (MQ=255)
     cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                  >  2:1644876/1‑150 (MQ=255)
     cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                  <  1:409190/150‑1 (MQ=255)
     cTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTgg                                                                                                                                                  <  1:154832/150‑1 (MQ=255)
      tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                 >  2:854173/1‑150 (MQ=255)
      tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                 <  1:544117/150‑1 (MQ=255)
        cGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAg                                                                                                                                               <  1:1624465/150‑1 (MQ=255)
            cTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGc                                                                                                                                           <  2:1116590/150‑1 (MQ=255)
              gcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTg                                                                                                                                         <  2:1485054/150‑1 (MQ=255)
               cgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGa                                                                                                                                        >  1:208458/1‑150 (MQ=255)
               cgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGa                                                                                                                                        <  1:201870/150‑1 (MQ=255)
               cgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGACTGCCAGCAAGCAGCCGCTGa                                                                                                                                        <  1:182201/150‑1 (MQ=255)
                 tgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAAc                                                                                                                                      >  1:1657896/1‑150 (MQ=255)
                  gcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACg                                                                                                                                     >  2:379088/1‑150 (MQ=255)
                  gcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGACTGCCAGCAAGCAGCCGCTGAACg                                                                                                                                     <  2:1685553/150‑1 (MQ=255)
                      tGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCt                                                                                                                                 <  2:1976922/150‑1 (MQ=255)
                      tGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCt                                                                                                                                 <  2:339118/150‑1 (MQ=255)
                        gCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCttt                                                                                                                               <  2:1538961/150‑1 (MQ=255)
                         cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                              >  1:574606/1‑150 (MQ=255)
                         cGGAAGTTTACGCCAGTAGCAATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                              <  2:749778/150‑1 (MQ=255)
                            aaGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCa                                                                                                                           <  1:423789/150‑1 (MQ=255)
                              gTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc                                                                                                                         <  2:1825297/150‑1 (MQ=255)
                               tttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                                                                        <  1:1782214/150‑1 (MQ=255)
                               tttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCg                                                                                                                        <  2:218232/150‑1 (MQ=255)
                                ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGt                                                                                                                       >  2:495521/1‑150 (MQ=255)
                                ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGt                                                                                                                       <  1:69598/150‑1 (MQ=255)
                                 tACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGta                                                                                                                      <  2:1005549/150‑1 (MQ=255)
                                 tACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGta                                                                                                                      >  1:1467922/1‑150 (MQ=255)
                                   cGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtata                                                                                                                    >  1:1107234/1‑150 (MQ=255)
                                        gTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc                                                                                                                         <  2:1049454/140‑1 (MQ=255)
                                        gTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc                                                                                                                         >  1:1049454/1‑140 (MQ=255)
                                            cGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaa                                                                                                           >  2:1873887/1‑150 (MQ=255)
                                                 cccACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaa                                                                                                      >  1:1771511/1‑150 (MQ=255)
                                                    aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                                                                                   >  2:28650/1‑150 (MQ=255)
                                                     cGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                                                                                  <  1:1731931/150‑1 (MQ=255)
                                                       agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTaca                                                                                                <  1:1873887/150‑1 (MQ=255)
                                                       agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTaca                                                                                                >  1:1320289/1‑150 (MQ=255)
                                                           gTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTAGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACt                                                                                            >  1:307833/1‑150 (MQ=255)
                                                            tttGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTg                                                                                           >  2:1718552/1‑150 (MQ=255)
                                                             ttGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTgg                                                                                          >  1:52514/1‑150 (MQ=255)
                                                              tGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGTCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGTCGATCGAAAACTACACACTGGa                                                                                         <  2:642318/150‑1 (MQ=255)
                                                               gTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAt                                                                                        >  2:600860/1‑150 (MQ=255)
                                                               gTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAt                                                                                        >  1:1664540/1‑150 (MQ=255)
                                                                 tAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtct                                                                                      >  1:1146514/1‑150 (MQ=255)
                                                                    aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctctc                                                                                   <  1:241469/150‑1 (MQ=255)
                                                                    aGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctctc                                                                                   <  2:556912/150‑1 (MQ=255)
                                                                     gACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTgg                                                                                          >  2:654311/1‑142 (MQ=255)
                                                                     gACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                                                                                  >  1:1296244/1‑150 (MQ=255)
                                                                     gACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCa                                                                                  <  2:1706257/150‑1 (MQ=255)
                                                                       cTTCGTGGAGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATg                                                                                <  1:358477/150‑1 (MQ=255)
                                                                                ggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTATCATTGTATAGTGGTCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCa                                                                       <  2:1574144/150‑1 (MQ=255)
                                                                                  cATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGTATCTCTCATGTCTGCCGCAgg                                                                     <  2:20245/150‑1 (MQ=255)
                                                                                   aTGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGa                                                                    <  2:710714/150‑1 (MQ=255)
                                                                                         gAAAGTGATGAACCTCGTCCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGc                                                              <  1:959126/150‑1 (MQ=255)
                                                                                         gAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGc                                                              >  1:1687693/1‑150 (MQ=255)
                                                                                         gAAAGTGATGAACCTCGACCGGTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGc                                                              <  1:1811942/150‑1 (MQ=255)
                                                                                            aGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg                                                                      <  2:1866085/139‑1 (MQ=255)
                                                                                            aGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg                                                                      >  1:1866085/1‑139 (MQ=255)
                                                                                            aGTGATGAACCTAGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAc                                                           <  2:307833/150‑1 (MQ=255)
                                                                                             gTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc                                                          <  2:1314797/150‑1 (MQ=255)
                                                                                                atgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCAc                                                       >  2:588120/1‑150 (MQ=255)
                                                                                                 tgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACt                                                      >  1:932157/1‑150 (MQ=255)
                                                                                                   aaCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctc                                                                                     >  2:1594547/1‑117 (MQ=255)
                                                                                                    aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                   >  2:1407879/1‑150 (MQ=255)
                                                                                                    aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                   <  1:854173/150‑1 (MQ=255)
                                                                                                        cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATc                                               <  1:1736840/150‑1 (MQ=255)
                                                                                                        cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATc                                               <  1:487610/150‑1 (MQ=255)
                                                                                                          aCCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAg                                             <  2:1321930/150‑1 (MQ=255)
                                                                                                           ccGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGt                                            <  2:532786/150‑1 (MQ=255)
                                                                                                           ccGTTTCAACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAAcc                                                          >  2:1424735/1‑136 (MQ=255)
                                                                                                            cGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTg                                           <  2:1778225/150‑1 (MQ=255)
                                                                                                            cGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTg                                           <  1:645814/150‑1 (MQ=255)
                                                                                                            cGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTg                                           <  1:1428374/150‑1 (MQ=255)
                                                                                                                 cGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGTGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTg                                      >  1:334075/1‑150 (MQ=255)
                                                                                                                     ctgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtg                                  <  2:1664540/150‑1 (MQ=255)
                                                                                                                               cTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACt                        <  2:902987/150‑1 (MQ=255)
                                                                                                                                 aagaacGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTa                      <  1:27739/145‑1 (MQ=255)
                                                                                                                                   ccccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATc                    >  1:1529842/1‑150 (MQ=255)
                                                                                                                                    cccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCt                   >  1:159190/1‑150 (MQ=255)
                                                                                                                                      cGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTgg                 <  2:580512/150‑1 (MQ=255)
                                                                                                                                      cGTTCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTgg                 >  1:1208579/1‑150 (MQ=255)
                                                                                                                                       gTTCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGa                >  2:957395/1‑150 (MQ=255)
                                                                                                                                        ttCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAg               >  2:610045/1‑150 (MQ=255)
                                                                                                                                        ttCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGc                                                                           >  1:778896/1‑90 (MQ=37)
                                                                                                                                         tCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGt              >  1:1921963/1‑150 (MQ=255)
                                                                                                                                          cAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGtt             >  1:1560294/1‑150 (MQ=255)
                                                                                                                                          cAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACAAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGtt             >  1:1616654/1‑150 (MQ=255)
                                                                                                                                           aGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTta            >  2:1741139/1‑150 (MQ=255)
                                                                                                                                             cGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTtata          <  2:1079230/150‑1 (MQ=255)
                                                                                                                                               gctgctTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGc        <  1:588120/150‑1 (MQ=255)
                                                                                                                                                  gcttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTgg     >  1:1406021/1‑150 (MQ=255)
                                                                                                                                                   cttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACGTATCTGGAGTTATAGCTGGa    >  2:1960057/1‑150 (MQ=255)
                                                                                                                                                    ttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCACCCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGAt   >  2:1252521/1‑150 (MQ=255)
                                                                                                                                                     tgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGAtt  >  2:775527/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TCTAACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATT  >  NC_000913/4135890‑4136188

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: