Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,165,150 A→G 8.1% intergenic (+152/‑39) 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,1500AG 8.1% 175.2 / 5.6 62intergenic (+152/‑39)baeR/yegPDNA‑binding transcriptional activator BaeR/DUF1508 domain‑containing protein YegP
Reads supporting (aligned to +/- strand):  ref base A (20/37);  new base G (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 = 9.56e-01

TAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAA  >  NC_000913/2165002‑2165295
                                                                                                                                                    |                                                                                                                                                 
tAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGAt                                                                                                                                                  <  1:1115418/150‑1 (MQ=255)
        ttattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCACCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGtttt                                                                                                                                          >  2:791879/1‑150 (MQ=255)
         tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                         <  2:73853/150‑1 (MQ=255)
         tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                         >  2:162842/1‑150 (MQ=255)
         tattTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTc                                                                                                                                         <  2:1290263/150‑1 (MQ=255)
          gttttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCACCAGAAAATCGCTCATACTTTAATCGGTAAACAGCATCTAAAGGTGCTGtttacc                                                                                                                                        <  1:968791/149‑4 (MQ=255)
          attTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCATCTAAAGGTGCTGtttacc                                                                                                                                        <  1:1252336/150‑4 (MQ=255)
            tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                      <  1:1067627/150‑1 (MQ=255)
            tttGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGAt                                                                                                                                      <  2:1054393/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:162842/150‑1 (MQ=255)
                                                                              gACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCa                                                                    <  1:1049965/150‑1 (MQ=255)
                                                                                   tATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTcc                                                               <  1:47896/150‑1 (MQ=255)
                                                                                   tATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTcc                                                               <  2:396530/150‑1 (MQ=255)
                                                                                     tctctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCgg                                                             <  1:1412354/150‑1 (MQ=255)
                                                                                     tctctcCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCgg                                                             <  1:1240729/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)
                                                                                                                                                    |                                                                                                                                                 
TAGCGACATTATTTTGTTAGCCGGAGATGGCGTTCAGACCAGATCCGGCAACATTATCCCACGCATGGTCAGCAAACTGACATTATCTCTCCGCCCTGGCAGGCATCAGAAAATCGCTCATACTTTAATCGGTAAACAGCACCTTTAGATGCTGTTTTCGATACACAATTTCAATCAAGGAGTCATTATGGCTGGTTGGTTTGAACTCAGTAAGAGCAGTGATAATCAGTTCCGGTTTGTGCTAAAAGCGGGCAATGGTGAGACTATCCTCACCAGCGAGCTTTATACCTCAAA  >  NC_000913/2165002‑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: