Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,136,042 T→A 5.6% 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.6% 580.6 / 14.5 176intergenic (+27/‑66)katG/yijEhydroperoxidase I/cystine exporter
Reads supporting (aligned to +/- strand):  ref base T (79/87);  new base A (7/3);  total (86/90)
Fisher's exact test for biased strand distribution p-value = 2.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 186 total aligned reads displayed.

AACTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGATTTTC  >  NC_000913/4135893‑4136191
                                                                                                                                                     |                                                                                                                                                     
aaCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgct                                                                                                                                                       >  2:713130/1‑150 (MQ=255)
aaCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCtgcttgct                                                                                                                                                       >  1:1423110/1‑150 (MQ=255)
 aCTCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTg                                                                                                                                                      <  2:118147/150‑1 (MQ=255)
   tCCGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGc                                                                                                                                                    <  1:709947/150‑1 (MQ=255)
    ccGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCa                                                                                                                                                   >  1:649386/1‑150 (MQ=255)
     cGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAg                                                                                                                                                  <  1:1282914/150‑1 (MQ=255)
     cGTCCTGCGTGCGGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAg                                                                                                                                                  <  1:1288055/150‑1 (MQ=255)
          tgcgtgcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCt                                                                                                                                             >  1:997013/1‑150 (MQ=255)
               gcgGTGGCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACg                                                                                                                                        >  1:872866/1‑150 (MQ=255)
                    ggCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAACGACTTCGTGGCGGCATGGGTGAAAGTGAGGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTGCAGCGACTGCCAGCAAGCAGCCGCTGAACggggtc                                                                                                                                   <  1:1324405/150‑6 (MQ=38)
                     gCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCttt                                                                                                                                  <  2:413526/150‑1 (MQ=255)
                     gCGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCttt                                                                                                                                  >  1:491503/1‑150 (MQ=255)
                      cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                                 >  2:687785/1‑150 (MQ=255)
                      cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                                 <  2:491503/150‑1 (MQ=255)
                      cGGAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTa                                                                                                                                 <  2:347283/150‑1 (MQ=255)
                       ggAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTAc                                                                                                                                >  1:757361/1‑150 (MQ=255)
                       ggAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTAc                                                                                                                                <  2:110410/150‑1 (MQ=255)
                        gAAGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTAcc                                                                                                                               <  2:234332/150‑1 (MQ=255)
                          aGTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAg                                                                                                                             >  1:682922/1‑150 (MQ=255)
                           gTTTACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGc                                                                                                                            >  1:246210/1‑150 (MQ=255)
                             ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCTAACTGCTTTAATCTGACACCGTTCAGAGGCTGATTGCTGGCATTCGCTGAACTTTATTTACCAGctt                                                                                                                          <  2:1194729/150‑3 (MQ=255)
                             ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGaa                                                                                                                                          >  2:181335/1‑134 (MQ=255)
                             ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGaa                                                                                                                                          <  1:181335/134‑1 (MQ=255)
                             ttACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGt                                                                                                                          >  1:590670/1‑150 (MQ=255)
                              tACGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGta                                                                                                                         >  2:385726/1‑150 (MQ=255)
                               aCGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtat                                                                                                                        >  2:221707/1‑150 (MQ=255)
                                cGCCAGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGtata                                                                                                                       >  1:356908/1‑150 (MQ=255)
                                    aGTAGCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTAGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTgg                                                                                                                   >  1:547186/1‑150 (MQ=255)
                                        gCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCt                                                                                                                                             <  2:358266/120‑1 (MQ=255)
                                        gCGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcga                                                                                                               >  1:935680/1‑150 (MQ=255)
                                         cGATGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaa                                                                                                              <  2:1400384/150‑1 (MQ=255)
                                           aTGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaacg                                                                                                            >  1:1235907/1‑150 (MQ=255)
                                            tGCCCACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGcgaacga                                                                                                           >  2:730711/1‑150 (MQ=255)
                                               ccACGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGaaaa                                                                                                        <  1:740157/150‑1 (MQ=255)
                                                 aCGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACt                                                                                                      >  2:1378638/1‑150 (MQ=255)
                                                  cGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGTTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                                                                                     >  1:314545/1‑150 (MQ=255)
                                                  cGAGAAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTa                                                                                                     >  2:300380/1‑150 (MQ=255)
                                                   gagaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTac                                                                                                    <  1:480399/150‑1 (MQ=255)
                                                    agaAGTTTGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTaca                                                                                                   <  2:529619/150‑1 (MQ=255)
                                                           tGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCATTCGCTGAACGTTCTTTACAATAGTATATTGGGCGAATGAAAATTAAAGACTgtt                                                                                            <  2:649386/150‑3 (MQ=255)
                                                           tGTTAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGa                                                                                            <  2:954913/150‑1 (MQ=255)
                                                             ttAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtc                                                                                          <  1:79429/150‑1 (MQ=255)
                                                              tAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtct                                                                                         >  1:426242/1‑150 (MQ=255)
                                                              tAAAGACTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtct                                                                                         >  1:608849/1‑150 (MQ=255)
                                                                    cTTCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTACTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATg                                                                                   <  2:1401952/150‑1 (MQ=255)
                                                                      tCGTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTc                                                                                 <  1:169112/150‑1 (MQ=255)
                                                                        gTGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTg                                                                               <  1:713278/150‑1 (MQ=255)
                                                                        gTGGCGGCATGGGTGAAAGTGATGAACCTCGACCCTTTCGACCTGCTGTAATCTGACCCCGTTAAGCGTCTGCTTGCTGGCAGTCGCTGAAAGTTATTTACCAGCGTAAAGTGCCCGAACCAAAACTACACACTGGATCTCACAAGTCTg                                                                               >  2:311343/1‑150 (MQ=255)
                                                                         tGGCGGCATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGc                                                                              <  1:966611/150‑1 (MQ=255)
                                                                          ggcggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGcc                                                                             >  1:285115/1‑150 (MQ=255)
                                                                             ggcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCa                                                                          >  1:516230/1‑150 (MQ=255)
                                                                              gcATGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAg                                                                         <  2:1423110/150‑1 (MQ=255)
                                                                                aTGGGTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGa                                                                       <  2:66921/150‑1 (MQ=255)
                                                                                    gTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAaga                                                                   <  1:442202/150‑1 (MQ=255)
                                                                                    gTGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAaga                                                                   >  1:203744/1‑150 (MQ=255)
                                                                                     tGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAagag                                                                  <  2:779409/150‑1 (MQ=255)
                                                                                     tGAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTAGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAagag                                                                  >  2:962400/1‑150 (MQ=255)
                                                                                      gAAAGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGGAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGTTCTCTCATGTCTGCCGGTGGAAAGAGc                                                                 <  2:1374764/150‑1 (MQ=255)
                                                                                        aaGTGATGAACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCaa                                                               >  2:904141/1‑150 (MQ=255)
                                                                                            gatgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctctc                                                                                      >  1:106374/1‑123 (MQ=255)
                                                                                            gatgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGAtctctc                                                                                      <  2:106374/123‑1 (MQ=255)
                                                                                            gatgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCa                                                           >  2:1045522/1‑150 (MQ=255)
                                                                                             atgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCAc                                                          >  1:576372/1‑150 (MQ=255)
                                                                                              tgaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACt                                                         <  1:1356530/150‑1 (MQ=255)
                                                                                               gaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTg                                                        >  2:873277/1‑150 (MQ=255)
                                                                                               gaACCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTg                                                        <  1:431592/150‑1 (MQ=255)
                                                                                                 aCCTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGc                                                      <  2:757361/150‑1 (MQ=255)
                                                                                                  ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCa                                                     <  2:516230/150‑1 (MQ=255)
                                                                                                  ccTCGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCa                                                     <  1:899654/150‑1 (MQ=255)
                                                                                                     cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAAta                                                  >  2:159335/1‑149 (MQ=255)
                                                                                                     cGACCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATc                                                  <  1:958347/150‑1 (MQ=255)
                                                                                                       aCCGTTTCGACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAg                                                >  2:629665/1‑150 (MQ=255)
                                                                                                               gACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTgg                                        >  1:826537/1‑150 (MQ=255)
                                                                                                               gACCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTgg                                        <  1:687785/150‑1 (MQ=255)
                                                                                                                aCCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGt                                       <  1:766686/150‑1 (MQ=255)
                                                                                                                aCCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGt                                       >  1:33938/1‑150 (MQ=255)
                                                                                                                aCCTGCTGTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGt                                       <  1:1034453/150‑1 (MQ=255)
                                                                                                                  ctgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtg                                     <  1:1378638/150‑1 (MQ=255)
                                                                                                                   tgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtgt                                    >  2:445888/1‑150 (MQ=255)
                                                                                                                   tgctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtgt                                    >  2:680790/1‑150 (MQ=255)
                                                                                                                    gctgTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTtgtg                                   <  2:203744/150‑1 (MQ=255)
                                                                                                                       gTAATCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTc                                >  1:304365/1‑150 (MQ=255)
                                                                                                                         aaTCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTcac                              >  1:1364128/1‑150 (MQ=255)
                                                                                                                         aaTCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTcac                              >  1:233264/1‑150 (MQ=255)
                                                                                                                           tCTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCacac                            >  1:210592/1‑150 (MQ=255)
                                                                                                                            cTGACCCCGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACt                           >  2:637905/1‑150 (MQ=255)
                                                                                                                                ccccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTAtg                       >  1:950699/1‑149 (MQ=255)
                                                                                                                                ccccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATc                       <  2:256044/150‑1 (MQ=255)
                                                                                                                                 cccGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCt                      <  1:221707/150‑1 (MQ=255)
                                                                                                                                  acGTTCAGCGACTGCCAGCAAGCAGCCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTg                     >  2:1311597/2‑150 (MQ=255)
                                                                                                                                   cGTTCAGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTgg                    <  1:637905/150‑1 (MQ=255)
                                                                                                                                        aGCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTta               >  2:7282/1‑150 (MQ=255)
                                                                                                                                         gCGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTtat              <  2:16988/150‑1 (MQ=255)
                                                                                                                                          cGGCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTtata             >  2:1271616/1‑150 (MQ=255)
                                                                                                                                           ggCTGCTTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTtataa            >  1:628759/1‑149 (MQ=255)
                                                                                                                                            gctgctTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACt                                                         <  1:1045864/104‑1 (MQ=255)
                                                                                                                                             ctgctTGCTGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCt          >  1:374882/1‑150 (MQ=255)
                                                                                                                                              tgcttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCGTGGTTGTGCTCACACTTATCTGGAGTTATAGCTg         <  2:1358402/150‑1 (MQ=255)
                                                                                                                                               gcttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTgg        <  2:374882/150‑1 (MQ=255)
                                                                                                                                                cttgctGGCAGTCGCTGAACGTTCTTTACCAGCGTATAGTGGGCGAACGAAAACTACACACTGGATCTCTCATGTCTGCCGCAGGAAAGAGCAACCCACTGGCAATCAGTGGCCTGGTTGTGCTCACACTTATCTGGAGTTATAGCTGGa       >  1:1138066/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: