Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,147 T→A 8.1% intergenic (+149/‑42) 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,1470TA 8.1% 181.5 / 6.1 62intergenic (+149/‑42)baeR/yegPDNA‑binding transcriptional activator BaeR/DUF1508 domain‑containing protein YegP
Reads supporting (aligned to +/- strand):  ref base T (20/37);  new base A (3/2);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 3.50e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.20e-01

GTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAA  >  NC_000913/2164998‑2165295
                                                                                                                                                     |                                                                                                                                                    
gTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCttt                                                                                                                                                      >  2:1025328/1‑150 (MQ=255)
gTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCttt                                                                                                                                                      <  2:698392/150‑1 (MQ=255)
 ttttAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTa                                                                                                                                                     <  1:587032/150‑1 (MQ=255)
    tAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGAt                                                                                                                                                  <  1:1115418/150‑1 (MQ=255)
            ttattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCACCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGtttt                                                                                                                                          >  2:791879/1‑150 (MQ=255)
             tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                         <  2:73853/150‑1 (MQ=255)
             tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                         <  2:1290263/150‑1 (MQ=255)
             tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                         >  2:162842/1‑150 (MQ=255)
              gttttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCACCAGAAAATCGCTCATACTTTAATCGGTAAACAGCATCTAAAGGTGCTGtttacc                                                                                                                                        <  1:968791/149‑4 (MQ=255)
              attTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCATCTAAAGGTGCTGtttacc                                                                                                                                        <  1:1252336/150‑4 (MQ=255)
                tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                      <  2:1054393/150‑1 (MQ=255)
                tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                      <  1:1067627/150‑1 (MQ=255)
                 ttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATa                                                                                                                                     <  1:257931/150‑1 (MQ=255)
                  tGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCCTGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATac                                                                                                                                    <  2:934538/150‑1 (MQ=255)
                     tAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTCCGATAcaca                                                                                                                                 >  1:1149137/1‑150 (MQ=255)
                      aGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACaa                                                                                                                                >  1:397981/1‑150 (MQ=255)
                          ggAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtc                                                                                                                            <  1:1150046/150‑1 (MQ=255)
                           gagaTGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtca                                                                                                                           <  2:199557/150‑1 (MQ=255)
                             gaTGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTtcaat                                                                                                                         <  2:529209/150‑1 (MQ=255)
                                   gTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGa                                                                                                                   >  2:56602/1‑150 (MQ=255)
                                     tCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGt                                                                                                                 <  2:800373/150‑1 (MQ=255)
                                        gACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCAt                                                                                                              >  2:1331212/1‑150 (MQ=255)
                                          ccAGATCCGGCAACATTATCCCACGCATGGTCAGCAATCTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTa                                                                                                            >  2:1412354/1‑150 (MQ=255)
                                               tCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCt                                                                                                       <  2:881644/150‑1 (MQ=255)
                                                     aaCATTATCCGACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCtggttgg                                                                                                 <  1:996107/150‑1 (MQ=255)
                                                       cATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTggttggtt                                                                                               >  2:1012345/1‑150 (MQ=255)
                                                       cATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTggttggtt                                                                                               >  2:1012346/1‑150 (MQ=255)
                                                         ttATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTg                                                                                             <  1:682576/150‑1 (MQ=255)
                                                           aTCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGaa                                                                                           <  1:353296/150‑1 (MQ=255)
                                                            tCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAAc                                                                                          <  2:561115/150‑1 (MQ=255)
                                                                   cATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTaa                                                                                   >  1:789028/1‑150 (MQ=255)
                                                                    aTGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAag                                                                                  <  1:290177/150‑1 (MQ=255)
                                                                          aGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGt                                                                            >  1:763946/1‑150 (MQ=255)
                                                                            cAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCCGTGa                                                                          <  1:1012346/150‑1 (MQ=255)
                                                                            cAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGa                                                                          <  1:1012345/150‑1 (MQ=255)
                                                                              aaCTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATa                                                                        >  1:1203194/1‑150 (MQ=255)
                                                                                cTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAAt                                                                      >  2:421880/1‑150 (MQ=255)
                                                                                 tGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATc                                                                     <  2:384302/150‑1 (MQ=255)
                                                                                  gACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCa                                                                    <  1:1049965/150‑1 (MQ=255)
                                                                                  gACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCa                                                                    <  1:162842/150‑1 (MQ=255)
                                                                                       tATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTcc                                                               <  1:47896/150‑1 (MQ=255)
                                                                                       tATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTcc                                                               <  2:396530/150‑1 (MQ=255)
                                                                                         tctctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCgg                                                             <  1:1240729/150‑1 (MQ=255)
                                                                                         tctctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCgg                                                             <  1:1412354/150‑1 (MQ=255)
                                                                                            ctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGttt                                                          >  2:740190/1‑150 (MQ=255)
                                                                                              ccGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTtgt                                                        >  2:533866/1‑150 (MQ=255)
                                                                                               cGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTtgtg                                                       <  2:638046/150‑1 (MQ=255)
                                                                                                       caggcaTCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGc                                               >  2:126025/1‑150 (MQ=255)
                                                                                                           caTCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGc                                           <  1:443230/150‑1 (MQ=255)
                                                                                                              cAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAAt                                        >  2:454656/1‑150 (MQ=255)
                                                                                                                gAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATgg                                      >  2:938734/1‑150 (MQ=255)
                                                                                                                   aaTCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTga                                   <  1:791879/150‑1 (MQ=255)
                                                                                                                    aTCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTgag                                  <  2:397981/150‑1 (MQ=255)
                                                                                                                        cTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTa                              >  1:332407/1‑150 (MQ=255)
                                                                                                                             aCTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTc                         <  1:296566/150‑1 (MQ=255)
                                                                                                                                  aaTCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAg                    <  2:763946/150‑1 (MQ=255)
                                                                                                                                    tCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCg                  <  1:436141/150‑1 (MQ=255)
                                                                                                                                     cGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGGATAGTCTCACCattgc                 <  1:629460/150‑5 (MQ=255)
                                                                                                                                      ggTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAg                >  2:585294/1‑150 (MQ=255)
                                                                                                                                       gaaaaCAGCATCTAAAGGTGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGc               >  2:1035565/3‑150 (MQ=255)
                                                                                                                                       gaaaaCAGAATCTAAAGGTGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGc               >  2:210164/3‑150 (MQ=255)
                                                                                                                                        aaaaCAGCATCTAAAGGTGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCt              >  1:270525/2‑150 (MQ=255)
                                                                                                                                         aaaCAGCACCTTTAGATGCTGTTTCCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCAAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCtt             <  2:1149137/150‑1 (MQ=255)
                                                                                                                                          aaCAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCttt            <  2:419034/150‑1 (MQ=255)
                                                                                                                                                    ttAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCaaa  <  2:1017671/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GTTTTAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAA  >  NC_000913/2164998‑2165295

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: ‑