Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,143 C→T 7.3% intergenic (+145/‑46) baeR → / → yegP DNA‑binding transcriptional activator BaeR/DUF1508 domain‑containing protein YegP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,165,1430CT 7.3% 290.1 / 8.3 95intergenic (+145/‑46)baeR/yegPDNA‑binding transcriptional activator BaeR/DUF1508 domain‑containing protein YegP
Reads supporting (aligned to +/- strand):  ref base C (43/45);  new base T (5/2);  total (48/47)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.12e-01

TTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTC  >  NC_000913/2164995‑2165292
                                                                                                                                                    |                                                                                                                                                      
ttAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCAcc                                                                                                                                                       <  2:313444/150‑1 (MQ=255)
   gTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCttt                                                                                                                                                    <  2:1427409/150‑1 (MQ=255)
        aGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATg                                                                                                                                               >  2:544504/1‑150 (MQ=255)
              attattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCATCTTTAGATGCTGttt                                                                                                                                         >  1:1578453/1‑150 (MQ=255)
              attattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCATCTAAAGGTGCTGttt                                                                                                                                         <  2:696311/150‑1 (MQ=255)
               ttattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGtttt                                                                                                                                       <  2:958255/150‑1 (MQ=255)
               ttattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACAATATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGtttt                                                                                                                                       >  2:1049462/1‑150 (MQ=255)
                tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCTCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                      <  1:50736/150‑1 (MQ=255)
                 attTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCg                                                                                                                                     >  2:1561225/1‑150 (MQ=255)
                 attTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCg                                                                                                                                     <  1:537892/150‑1 (MQ=255)
                  ttttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGa                                                                                                                                    <  2:188455/150‑1 (MQ=255)
                   tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACGGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                   >  2:181472/1‑150 (MQ=255)
                   tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                   <  1:302144/150‑1 (MQ=255)
                   tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                   <  1:699333/150‑1 (MQ=255)
                   tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGAAAACAGCATCTAAAGGTGCTGtttaccgat                                                                                                                                   <  1:1634230/150‑7 (MQ=255)
                      gTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATaca                                                                                                                                >  1:146984/1‑150 (MQ=255)
                      gTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATaca                                                                                                                                >  1:786752/1‑150 (MQ=255)
                      gTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATaca                                                                                                                                >  1:121131/1‑150 (MQ=255)
                        tAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATAcaca                                                                                                                              >  2:1356568/1‑150 (MQ=255)
                           ccGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAAtt                                                                                                                           >  2:1290591/1‑150 (MQ=255)
                             ggAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtc                                                                                                                         <  1:628974/150‑1 (MQ=255)
                             ggAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtc                                                                                                                         <  1:623209/150‑1 (MQ=255)
                               agaTGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtcaa                                                                                                                       <  1:548156/150‑1 (MQ=255)
                               agaTGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtcaa                                                                                                                       <  2:973606/150‑1 (MQ=255)
                                   ggCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtcaatcaa                                                                                                                   >  1:1563914/1‑150 (MQ=255)
                                     cGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAgg                                                                                                                 <  2:1643098/150‑1 (MQ=255)
                                     cGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAgg                                                                                                                 >  2:950256/1‑150 (MQ=255)
                                      gTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGa                                                                                                                >  2:1316107/1‑150 (MQ=255)
                                       ttCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAg                                                                                                               <  1:2993/150‑1 (MQ=255)
                                          aGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCa                                                                                                            <  2:148237/150‑1 (MQ=255)
                                           gACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCAt                                                                                                           >  1:90055/1‑150 (MQ=255)
                                             ccAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATAGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTa                                                                                                         >  2:41390/1‑150 (MQ=255)
                                               aGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAATTCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATg                                                                                                       >  2:1567667/1‑150 (MQ=255)
                                                gATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATgg                                                                                                      >  1:431283/1‑150 (MQ=255)
                                                gATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATAAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATAAAGGAATCATTATga                                                                                                      >  2:709682/1‑149 (MQ=255)
                                                    cGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCtgg                                                                                                  <  1:635857/150‑1 (MQ=255)
                                                      gCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGTACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCtggtt                                                                                                <  2:706553/150‑1 (MQ=255)
                                                          cATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTggttggtt                                                                                            <  1:729972/150‑1 (MQ=255)
                                                            ttATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTg                                                                                          <  1:857936/150‑1 (MQ=255)
                                                               tCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAAc                                                                                       <  2:1526804/150‑1 (MQ=255)
                                                               tCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAAc                                                                                       <  1:465984/150‑1 (MQ=255)
                                                                  cACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCa                                                                                    >  2:949133/1‑150 (MQ=255)
                                                                  cACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCa                                                                                    <  1:1595094/150‑1 (MQ=255)
                                                                   aCGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAg                                                                                   <  1:1426337/150‑1 (MQ=255)
                                                                      cATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTaa                                                                                >  1:611441/1‑150 (MQ=255)
                                                                      cATCGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTaa                                                                                >  2:1342543/1‑150 (MQ=255)
                                                                            cAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAg                                                                          >  1:705443/1‑150 (MQ=255)
                                                                             aGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGt                                                                         <  1:994712/150‑1 (MQ=255)
                                                                               cAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGCAAGAGCAGTGa                                                                       <  2:1563914/150‑1 (MQ=255)
                                                                                    tGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATc                                                                  <  1:1376177/150‑1 (MQ=255)
                                                                                      aCATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAg                                                                <  1:1290591/150‑1 (MQ=255)
                                                                                      aCATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAg                                                                <  1:500416/150‑1 (MQ=255)
                                                                                       cATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGt                                                               >  1:505311/1‑150 (MQ=255)
                                                                                       aaTTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATAGGTAAACAGCACCTTTAGATTCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGt                                                               <  1:1102604/149‑1 (MQ=255)
                                                                                        aTTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGtt                                                              >  2:62441/1‑150 (MQ=255)
                                                                                          tATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTct                                                            <  1:709682/150‑2 (MQ=255)
                                                                                          tATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTcc                                                            <  2:181486/150‑1 (MQ=255)
                                                                                            tctctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCgg                                                          <  1:1356568/150‑1 (MQ=255)
                                                                                              tctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCTCCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGt                                                                                  >  1:757532/1‑124 (MQ=255)
                                                                                              tctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCTCCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGt                                                                                  <  2:757532/124‑1 (MQ=255)
                                                                                              tctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGtt                                                        >  1:178973/1‑150 (MQ=255)
                                                                                               ctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGttt                                                       >  2:621500/1‑150 (MQ=255)
                                                                                               ctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATACTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGttt                                                       >  2:125955/1‑150 (MQ=255)
                                                                                                  cGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTtgtg                                                    <  1:41390/150‑1 (MQ=255)
                                                                                                   gCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGc                                                   >  1:138169/1‑150 (MQ=255)
                                                                                                      cTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTaa                                                >  1:331699/1‑150 (MQ=255)
                                                                                                            ggcaTCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCgg                                          >  2:846848/1‑150 (MQ=255)
                                                                                                                 cAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAAt                                     >  2:1233277/1‑150 (MQ=255)
                                                                                                                   taaaaTCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATgg                                   <  1:1316107/149‑1 (MQ=255)
                                                                                                                   gAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATgg                                   >  2:134071/1‑150 (MQ=255)
                                                                                                                     aaaTCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTg                                 <  2:1631469/150‑1 (MQ=255)
                                                                                                                       aTCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTgag                               <  1:950256/150‑1 (MQ=255)
                                                                                                                          gCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACt                            >  2:13002/1‑150 (MQ=255)
                                                                                                                          gCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACt                            >  2:19205/1‑150 (MQ=255)
                                                                                                                              aTACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATcc                        >  1:1482295/1‑150 (MQ=255)
                                                                                                                              aTACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATcc                        <  1:1582614/150‑1 (MQ=255)
                                                                                                                              aTACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATcc                        <  2:178973/150‑1 (MQ=255)
                                                                                                                              aTACATTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATcc                        >  1:1482297/1‑150 (MQ=255)
                                                                                                                                 cTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCtac                     >  2:176284/1‑148 (MQ=255)
                                                                                                                                     gtatcgaaaaCAGCATCTAAAGGTGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAg                 >  2:606147/8‑150 (MQ=255)
                                                                                                                                      tatcgaaaaCAGCATCTAAAGGTGCTGTTTACCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAg                 >  2:441901/7‑150 (MQ=38)
                                                                                                                                       atcgaaaaCAGCATCTAAAGGTGCTGTTTACGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCg               >  2:672589/6‑150 (MQ=255)
                                                                                                                                        cGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGa              >  1:1390485/1‑150 (MQ=255)
                                                                                                                                        cGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAAACAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGa              >  1:1084550/1‑150 (MQ=255)
                                                                                                                                           aaaaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCt           <  2:1486858/149‑1 (MQ=255)
                                                                                                                                            aaaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCTGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCtt          <  1:521925/150‑1 (MQ=255)
                                                                                                                                            aaaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCtt          >  1:283162/1‑150 (MQ=255)
                                                                                                                                            aaaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCtt          <  1:1567667/150‑1 (MQ=255)
                                                                                                                                            aaaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGGTCCGGGTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCCCCAGCGAGCtt          <  2:1070329/150‑1 (MQ=255)
                                                                                                                                             aaCAGCATCTAAAGGTGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCttt         >  2:1281218/1‑150 (MQ=255)
                                                                                                                                             aaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCttt         <  1:1233277/150‑1 (MQ=255)
                                                                                                                                               cAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTtat       >  2:685915/1‑150 (MQ=255)
                                                                                                                                               cAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTtat       >  1:1550702/1‑150 (MQ=255)
                                                                                                                                                 gCACCTTTAGATGCTCTTTTCGATACACAATATCAATCAAGCAGACATTATGGCTGGTAGCTATGAACTCAGAATGAGCAGTCATAAACAGTTCCGGTTTGTGCTAAAAGCGGGCATTCGACAGACTATCCTCACCACCGAGCTTTATAc     >  2:1441286/1‑150 (MQ=255)
                                                                                                                                                 gCACCTTTAGATACTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAc     <  1:125955/150‑1 (MQ=255)
                                                                                                                                                  cACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGGCATTATGGCGGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGATCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATAcc    <  1:767097/150‑1 (MQ=255)
                                                                                                                                                    ccTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACATc  >  2:1447513/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                      
TTAGTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTC  >  NC_000913/2164995‑2165292

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: