Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,073,672 A→T 26.4% F16I (TTC→ATC)  scoC ← transcriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,073,6720AT26.4% 513.6 / 209.3 278F16I (TTC→ATC) scoCtranscriptional regulator of extracellular protease production, sporulation and bacilysin production (MarR family)
Reads supporting (aligned to +/- strand):  ref base A (83/121);  new base T (36/37);  total (119/159)
Fisher's exact test for biased strand distribution p-value = 2.17e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

Only 100 of 285 total aligned reads displayed.

TTCAGAAATGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAACCATTGC  >  BS168_reference_genome/1073523‑1073821
                                                                                                                                                     |                                                                                                                                                     
ttCAGAAATGGAACCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGaa                                                                                                                                                       <  2:1502965/150‑1 (MQ=255)
  cAGAAATGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGaaaa                                                                                                                                                     >  1:2609259/1‑150 (MQ=255)
    gAAATGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAAcc                                                                                                                                                   >  1:4379062/1‑150 (MQ=255)
    gAAATGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAAcc                                                                                                                                                   <  1:4267187/150‑1 (MQ=255)
     aaaTGCAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCa                                                                                                                                                  <  1:1379239/150‑1 (MQ=255)
      aaTGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCag                                                                                                                                                 >  1:18418/1‑150 (MQ=255)
       aTGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCaga                                                                                                                                                <  1:5488462/150‑1 (MQ=255)
          gAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGct                                                                                                                                             <  2:3858853/150‑1 (MQ=255)
          gAAGCTCCATTAAATTGATACGCAATCAATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTACAAAGAGACTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGct                                                                                                                                             <  1:3295078/150‑1 (MQ=255)
               tCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCttt                                                                                                                                        <  1:4217647/150‑1 (MQ=255)
                ccATACAATTGATACCCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTAGCAAAGAGCCTTGCTAAGCTGAGTCATTTTCTGGGTGATAACCAGAGCTTCTTTc                                                                                                                                       <  2:162559/150‑1 (MQ=255)
                 cATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGGGTCATTTTCTGGGTGATAACCAGAGCTTCTTTca                                                                                                                                      >  1:2454126/1‑150 (MQ=255)
                     aaaTTCAGACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTGTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCGTCTTTCAcatc                                                                                                                                  <  2:2949557/149‑1 (MQ=255)
                      aaTTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCAcatca                                                                                                                                 >  2:2197031/1‑150 (MQ=255)
                       aTTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCAcatcat                                                                                                                                <  2:3359824/150‑1 (MQ=255)
                        ttGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATa                                                                                                                               >  2:4524907/1‑150 (MQ=255)
                         tGATACGCAATCCATAAAATTTGATGGTCATTAATATTCAGGTGATACGGTTTGTGGTATTGCTGTCAATACTTGTCCATCGATTTCCAAAGAGCCTTTCTAAGCTGAGTCATTTTGTGGGTGAAAACGAGAGCTTCTTTGACATCATAg                                                                                                                              <  2:4252522/150‑1 (MQ=255)
                          gATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGCGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCCTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAgg                                                                                                                             <  2:1233667/150‑1 (MQ=255)
                          gATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGGTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAgg                                                                                                                             >  2:2382832/1‑150 (MQ=255)
                           ntACNCAATCCATAAAATATGATGCNCATTAATATTCAGGTCATACGGTNTGAGCCATTNCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAggg                                                                                                                            >  1:3995394/2‑150 (MQ=255)
                           aTACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAggg                                                                                                                            <  2:2840874/150‑1 (MQ=255)
                           aTACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAggg                                                                                                                            >  2:3555251/1‑150 (MQ=255)
                                cAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGtt                                                                                                                       <  2:4953365/150‑1 (MQ=255)
                                   tCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCa                                                                                                                    >  2:1851037/1‑150 (MQ=255)
                                    ccATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCAc                                                                                                                   <  1:2101263/150‑1 (MQ=255)
                                    ccATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCAc                                                                                                                   >  1:589962/1‑150 (MQ=255)
                                     cATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACt                                                                                                                  >  1:3414920/1‑150 (MQ=255)
                                      aTAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTc                                                                                                                 <  2:60968/150‑1 (MQ=255)
                                      aTAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTc                                                                                                                 >  2:2550998/1‑150 (MQ=255)
                                         aaaTATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGAt                                                                                                              <  2:2088622/150‑1 (MQ=255)
                                          aaTATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGAtt                                                                                                             >  2:3082717/1‑150 (MQ=255)
                                            tatGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAGTCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCa                                                                                                           >  1:666175/1‑150 (MQ=255)
                                            tatGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCa                                                                                                           <  2:4742227/150‑1 (MQ=255)
                                             atgatgCACATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATTCTCCTCGATCGATTTCCAAAGAGTCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCAt                                                                                                          <  2:704460/150‑1 (MQ=255)
                                               gattCTTATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCTTTCTCGATCGATTTCGAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTa                                                                                                        <  1:4779144/150‑1 (MQ=255)
                                                atgCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTAc                                                                                                       >  1:3104273/1‑150 (MQ=255)
                                                   cTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTc                                                                                                    <  1:1005773/150‑1 (MQ=255)
                                                   cTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTc                                                                                                    >  2:3741201/1‑150 (MQ=255)
                                                   cTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTc                                                                                                    >  2:5243102/1‑150 (MQ=255)
                                                          atatTTAGGTCATATGGTTTGAGCCATTGCTGCCAATCCTTCTCCATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTACCACTCGATTCATTACGTCACCTGCt                                                                                             <  1:5270645/150‑1 (MQ=255)
                                                            atTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTtc                                                                                           <  2:1802480/150‑1 (MQ=255)
                                                             ttCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTtct                                                                                          <  2:2228385/150‑1 (MQ=255)
                                                             ttCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTtct                                                                                          >  1:2871042/1‑150 (MQ=255)
                                                             ttCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTtct                                                                                          >  2:3392651/1‑150 (MQ=255)
                                                              tCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTtctc                                                                                         >  1:3078282/1‑150 (MQ=255)
                                                                aGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCtctt                                                                                       >  2:4864715/1‑150 (MQ=255)
                                                                aGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCCCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCtctt                                                                                       <  1:5487810/150‑1 (MQ=255)
                                                                 ggTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCtcttt                                                                                      >  1:3206091/1‑150 (MQ=255)
                                                                   tCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCtcttttt                                                                                    >  1:3705195/1‑150 (MQ=255)
                                                                   tCATACGGTTTGAGCCATTGCTCACAATCCTTCTCGATCGATTTACAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTACACTGGATTCATTACGTCACCTGCTTCtcttttt                                                                                    <  1:2397401/150‑1 (MQ=255)
                                                                   gcATACGGTTTGAGCCTTTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCtcttttt                                                                                    <  1:3261700/149‑1 (MQ=255)
                                                                         ggTTTGACCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTctttttcttttt                                                                              <  1:707484/150‑1 (MQ=255)
                                                                             tGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTctct                                                                          >  2:4051602/1‑150 (MQ=255)
                                                                             tGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTctct                                                                          <  1:991580/150‑1 (MQ=255)
                                                                               aGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCtat                                                                        <  1:825466/150‑1 (MQ=255)
                                                                                gCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCtatt                                                                       <  2:861560/150‑1 (MQ=255)
                                                                                  cATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCtattat                                                                     >  1:913568/1‑150 (MQ=255)
                                                                                   aTTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATg                                                                    <  1:5718106/150‑1 (MQ=255)
                                                                                   aTTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATg                                                                    <  1:4708873/150‑1 (MQ=255)
                                                                                     tgctgcAAATCCTTCTCGATCGATTTCCAAAGAGCCTGGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTATCTATTATGcc                                                                  >  1:179291/1‑150 (MQ=255)
                                                                                        tgcCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCaat                                                               <  1:2820479/150‑1 (MQ=255)
                                                                                           cAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCaataaa                                                            >  1:4982519/1‑150 (MQ=255)
                                                                                            aaTCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCaataaaa                                                           >  1:5097196/1‑150 (MQ=255)
                                                                                            aaTCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCaataaaa                                                           >  2:2907889/1‑150 (MQ=255)
                                                                                             aTCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCaataaaat                                                          >  1:3798656/1‑150 (MQ=255)
                                                                                             aTCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCaataaaat                                                          <  2:5471496/150‑1 (MQ=255)
                                                                                                   ctcGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGAt                                                    <  1:4842393/150‑1 (MQ=255)
                                                                                                   ctcGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGAt                                                    <  2:5108730/150‑1 (MQ=255)
                                                                                                         cgatTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATtaataa                                              <  2:4379062/150‑1 (MQ=255)
                                                                                                          gatTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATaaa                                             >  2:5047541/1‑150 (MQ=255)
                                                                                                           atTTCCTAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAg                                            >  1:3182428/1‑150 (MQ=255)
                                                                                                           atTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAg                                            >  1:1447768/1‑150 (MQ=255)
                                                                                                               ccAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTCCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGACAATAAAATAAAGATTAATAAAGATGa                                        <  2:3328665/150‑1 (MQ=255)
                                                                                                                 aaaGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGCGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAACATGaaa                                      >  1:33730/1‑150 (MQ=255)
                                                                                                                     agCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGt                                  <  2:307353/150‑1 (MQ=255)
                                                                                                                        cTTGCTAAGCTGAGCCATTTTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAg                               >  2:3226158/1‑150 (MQ=255)
                                                                                                                         ttGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGa                              >  2:2628757/1‑150 (MQ=255)
                                                                                                                         ttGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGa                              >  2:5390648/1‑150 (MQ=255)
                                                                                                                          tGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGaa                             >  1:3504689/1‑150 (MQ=255)
                                                                                                                          tGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGaa                             >  1:3807370/1‑150 (MQ=255)
                                                                                                                                gCTGAGCCATATTCTGGGTGATAACCAGAGCTTCTTTCACATCATAGTGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTTTTTTTCTTTTTCTCTATTATGACAATAAAATAAAGATTAATAAAGATGAAAATGTAAAAAtttttt                       >  1:497115/1‑150 (MQ=255)
                                                                                                                                 cTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTta                      <  1:3555251/150‑1 (MQ=255)
                                                                                                                                  tGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTtaa                     >  1:4104802/1‑150 (MQ=255)
                                                                                                                                    aGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGACGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTCCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTtaata                   <  1:3982571/150‑1 (MQ=255)
                                                                                                                                     gtcATTTTCTGGGTGAAAACCTGAGCTTCTTTCACACCATAGGGTGGTTCCACTCGATTCATTACGCCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGTCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTtaataa                  <  2:1030278/148‑1 (MQ=255)
                                                                                                                                     gCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTtaataa                  <  1:2550998/150‑1 (MQ=255)
                                                                                                                                       cATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAAt                <  1:3741201/150‑1 (MQ=255)
                                                                                                                                       cATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAAt                <  2:3995394/150‑1 (MQ=255)
                                                                                                                                         ttttCTTGGTGAAAACCAGAGCTTCTTTCACATCATAGTGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTc              >  1:1549269/1‑150 (MQ=255)
                                                                                                                                          tttCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTca             <  1:1983711/150‑1 (MQ=255)
                                                                                                                                           ttCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTcac            <  1:1286610/150‑1 (MQ=255)
                                                                                                                                             cTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACaa          >  1:1717261/1‑150 (MQ=255)
                                                                                                                                             cTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACaa          >  1:2009462/1‑150 (MQ=255)
                                                                                                                                              tGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACaaa         >  2:5757193/1‑150 (MQ=255)
                                                                                                                                               gggTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAAc        >  2:1931758/1‑150 (MQ=255)
                                                                                                                                               gggTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAAc        >  1:2468895/1‑150 (MQ=255)
                                                                                                                                               gggTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAAc        >  1:1092502/1‑150 (MQ=255)
                                                                                                                                                  tGATAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAACCAt     <  2:2413648/150‑1 (MQ=255)
                                                                                                                                                    aaaaCCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAACCATTg   >  2:5813524/1‑150 (MQ=255)
                                                                                                                                                     aaaCCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAACCATTGc  <  2:1794782/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTCAGAAATGGAAGCTCCATTCAATTGATACGCAATCCATAAAATATGATGCTCATTAATATTCAGGTCATACGGTTTGAGCCATTGCTGCCAATCCTTCTCGATCGATTTCCAAAGAGCCTTGCTAAGCTGAGCCATTTTCTGGGTGAAAACCAGAGCTTCTTTCACATCATAGGGCGGTTCCACTCGATTCATTACGTCACCTGCTTCTCTTTTTCTTTTTCTCTATTATGCCAATAAAATAAAGATTAATAAAGATGAAAATGTAAGAATTTTTTAATAAATTCACAAACCATTGC  >  BS168_reference_genome/1073523‑1073821

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: