Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,361,982 C→T 5.5% R119R (CGG→CGA ypjB ← sporulation septum‑associated protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,361,9820CT 5.5% 313.4 / 4.1 91R119R (CGG→CGAypjBsporulation septum‑associated protein
Reads supporting (aligned to +/- strand):  ref base C (34/52);  new base T (3/2);  total (37/54)
Fisher's exact test for biased strand distribution p-value = 3.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAAATCAAGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGATTTGACGAATTTGGGCTCCCGTT  >  BS168_reference_genome/2361844‑2362120
                                                                                                                                          |                                                                                                                                          
gTAAATCAAGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCGAAAGCGGGTCAGAc                                                                                                                                            <  7:196308/139‑1 (MQ=255)
 tAAATCAAGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGAcc                                                                                                                                           <  3:110267/139‑1 (MQ=255)
   aaTCAAGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGa                                                                                                                                         >  1:36786/1‑139 (MQ=255)
        aGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGt                                                                                                                                     >  3:175586/1‑138 (MQ=255)
        aGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGt                                                                                                                                     >  7:20331/1‑138 (MQ=255)
        aGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTc                                                                                                                                    <  8:368433/139‑1 (MQ=255)
        aGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTc                                                                                                                                    <  5:160465/139‑1 (MQ=255)
        aGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCCTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGt                                                                                                                                     >  5:51633/1‑138 (MQ=255)
         gTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTc                                                                                                                                    >  2:387948/1‑138 (MQ=255)
               taAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGca                                                                                                                             <  7:44132/139‑1 (MQ=255)
                 caCTGATAAATTCATTCCATTTTTCATGACAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGGCAGACCGATTGTCAACCGcatc                                                                                                                           <  2:259019/138‑1 (MQ=255)
                 aaatGATAAATTGATTCCATTTTTCATGAAAGCTTGAGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGcatc                                                                                                                           <  4:61598/136‑1 (MQ=255)
                 aaCTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGcatc                                                                                                                           <  1:197091/139‑1 (MQ=255)
                    tGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGcatccat                                                                                                                        >  3:309903/1‑139 (MQ=255)
                      atatatTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTAGGTCATCTCGCTTCAGTTGAGTAATAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTGAAGCGgatccat                                                                                                                        <  4:175586/137‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTc                                                                                                                                    >  1:50299/1‑125 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTc                                                                                                                                    <  2:50299/125‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCg                                                                                                                               >  7:290363/1‑130 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCg                                                                                                                               <  8:290363/130‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGcatccat                                                                                                                        <  7:21833/137‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGcatccat                                                                                                                        <  3:296236/137‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGcatccat                                                                                                                        <  3:398002/137‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      <  3:232852/139‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      >  7:134110/1‑139 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      <  5:347947/139‑1 (MQ=255)
                      aTAAATTCATTCCATTTTTCATGAAAGCTCGTACTTTCATTCTTTTGGACATCCTGCTTCAGTTCAGCAAAGGCTTCCATAATCGGTCTTTCCAGTGAGCCCCAAAGCGGATCAGATCGATTGTCAACTGCGTCCATCa                                                                                                                      >  2:176046/1‑139 (MQ=255)
                       tatattcattcCATTTTTCATGAAAGCTTGTGCTTGCATTCTTTTGGACATCTCGCTTGAGTTCAGTAATAGCTTCGATTATCGGTTTTTGCAGTGAGCGCCAAAGCGGGTCAGACCGATTGTCAACCgcatgcatga                                                                                                                      <  8:29956/135‑3 (MQ=255)
                       tAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAg                                                                                                                     <  1:80406/139‑1 (MQ=255)
                           ttcattcCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTGCAGTAAAAGCTTCCATAATCGGGTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCAtt                                                                                                                 <  8:27247/139‑1 (MQ=255)
                                tcCATTTTTCATGAAAGCTCGTACTTTCATTCTTTTGGACATCCTGCTTCAGTTCAGCAAAGGCTTCCATAATCGGTCTTTCCAGTGAGCCCCAAAGCGGATCAGATCGATTGTCAACTGCGTCCATCAGCATTCTGaa                                                                                                            >  4:349901/1‑139 (MQ=255)
                                tcCATTTTTCATGAAAGCTCGTACTTTCATTCTTTTGGACATCCTGCTTCAGTTCAGCAAAGGCTTCCATAATCGGTCTTTCCAGTGAGCCCCAAAGCGGATCAGATCGATTGTCAACTGCGTCCATCAGCATTCTGaa                                                                                                            >  4:121739/1‑139 (MQ=255)
                                      tttCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGGTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTgcg                                                                                                       <  2:11395/138‑1 (MQ=255)
                                               aGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGctct                                                                                             <  1:369716/139‑1 (MQ=255)
                                                    gtgCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGt                                                                                         <  4:309903/138‑1 (MQ=255)
                                                         tCCATTCTTTTGGACATCTCGCTTCAGTACAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCTGGTCAGACCGATTGTCAACCGCAGCCATCAGCATTCTTAAATGCGCGGCTGCTCTGAGTATTtgt                                                                                   >  8:295634/1‑139 (MQ=255)
                                                          tcATTCTTTTGGACATCCTGCTTCAGTTCAGCAAAGGCTTCCATAATCGGTCTTTCCAGTGAGCCCCAAAGCGGATCAGATCGATTGTCAACTGCGTCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGCTTTtg                                                                                    <  1:176046/136‑1 (MQ=255)
                                                          ccATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGgctgct                                                                                               <  8:330682/126‑1 (MQ=255)
                                                          ccATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGgctgct                                                                                               >  7:330682/1‑126 (MQ=255)
                                                          ccATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTtgt                                                                                   <  3:120001/138‑1 (MQ=255)
                                                          ccATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTtg                                                                                    <  8:402610/137‑1 (MQ=255)
                                                               cTTTTGGACATCCTGCTTCAGTTCAGCAAAGGCTTCCATAATCGGTCTTTCCAGTGAGCCCCAAAGCGGATCAGATCGATTGTCAACTGCGTCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGCTTTTGTGTATc                                                                              <  3:121739/138‑1 (MQ=255)
                                                                  ttGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCt                                                                             <  3:142386/136‑1 (MQ=255)
                                                                  ttGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTg                                                                            <  1:387948/137‑1 (MQ=255)
                                                                  ttGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTg                                                                            <  8:134110/137‑1 (MQ=255)
                                                                  ttGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTGTTTCCAGTGAGCCCCAAAGAGGTGCAGACCGATTGTCAACCGCATCCATCAGAATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTg                                                                            <  2:36173/137‑1 (MQ=255)
                                                                          ctcGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTtgtgtg                                                                                <  3:224138/125‑1 (MQ=255)
                                                                          ctcGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTtgtgtg                                                                                >  4:224138/1‑125 (MQ=255)
                                                                          ctcGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  >  4:326322/1‑139 (MQ=255)
                                                                          ctcGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  >  7:413996/1‑139 (MQ=255)
                                                                          ctcGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  >  5:356816/1‑139 (MQ=255)
                                                                          ctcGCTTCAGTTCAGTAAAAGCTTCCATAATCCGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  >  8:257819/1‑139 (MQ=255)
                                                                             gCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgctt                                                                <  7:419223/138‑1 (MQ=255)
                                                                             gCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  <  4:243219/136‑1 (MQ=255)
                                                                             gCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  <  5:414197/136‑1 (MQ=255)
                                                                             gCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  <  7:257819/136‑1 (MQ=255)
                                                                                    ttcagtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGgc                                                                                                   >  6:187919/1‑96 (MQ=255)
                                                                                    ttcagtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGgc                                                                                                   <  5:187919/96‑1 (MQ=255)
                                                                                    ttcagtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTaaa                                                        >  2:28732/1‑139 (MQ=255)
                                                                                       agtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGca                                                                                                                             <  5:116761/67‑1 (MQ=255)
                                                                                       agtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      <  3:322582/74‑1 (MQ=255)
                                                                                       agtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      >  4:322582/1‑74 (MQ=255)
                                                                                       agtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTc                                                     <  8:106112/139‑1 (MQ=255)
                                                                                       agtAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTAAACCGca                                                                                                                             >  6:116761/1‑67 (MQ=255)
                                                                                         tAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      >  7:298422/1‑72 (MQ=255)
                                                                                         tAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      <  8:298422/72‑1 (MQ=255)
                                                                                         tAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCg                                                   <  3:59653/139‑1 (MQ=255)
                                                                                               cTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCat                                             <  1:339973/139‑1 (MQ=255)
                                                                                                 tCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGc                                                                                                                    <  3:288347/66‑1 (MQ=255)
                                                                                                 tCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGc                                                                                                                    >  4:288347/1‑66 (MQ=255)
                                                                                                  ccATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCatat                                           <  6:356816/138‑1 (MQ=255)
                                                                                                  ccATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCata                                            <  5:183794/137‑1 (MQ=255)
                                                                                                  ccATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGGCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCatat                                           <  5:340746/138‑1 (MQ=255)
                                                                                                   cATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCa                                         >  5:28460/1‑139 (MQ=255)
                                                                                                   cATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTCTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCa                                         >  8:109027/1‑139 (MQ=255)
                                                                                                       aTCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGACCTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGGATCAGCTTGCTTTAAACTCCGACCCATATCATtat                                     <  2:99001/139‑1 (MQ=255)
                                                                                                       aTCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATta                                      <  6:51633/138‑1 (MQ=255)
                                                                                                            tttttCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      <  1:262244/53‑1 (MQ=255)
                                                                                                            tttttCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCa                                                                                                                      >  2:262244/1‑53 (MQ=255)
                                                                                                                      gAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGAttt                      >  3:236459/1‑139 (MQ=255)
                                                                                                                      gAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGAttt                      >  3:185345/1‑139 (MQ=255)
                                                                                                                        gCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGATTTGa                    <  2:36786/139‑1 (MQ=255)
                                                                                                                         ccccAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  >  6:268189/1‑92 (MQ=255)
                                                                                                                         ccccAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  <  5:262740/92‑1 (MQ=255)
                                                                                                                         ccccAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  <  5:268189/92‑1 (MQ=255)
                                                                                                                         ccccAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCGCTGAGTTTTTGTGTGTCTGATGTATCAgc                                                                  >  6:262740/1‑92 (MQ=255)
                                                                                                                               aGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGgctgct                                                                                               <  5:336030/57‑1 (MQ=255)
                                                                                                                               aGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGgctgct                                                                                               >  6:336030/1‑57 (MQ=255)
                                                                                                                               aGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGATTTGACGAAttt             <  8:20331/139‑1 (MQ=255)
                                                                                                                                 cGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGATTTGACGAATTTgg           >  2:14008/1‑139 (MQ=255)
                                                                                                                                   ggTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTgtgtgt                                                                               >  4:23667/1‑69 (MQ=255)
                                                                                                                                   ggTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTgtgtgt                                                                               <  3:23667/69‑1 (MQ=255)
                                                                                                                                    gTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTCATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGAGATTTGACGAATTAGGGc         >  1:345686/1‑138 (MQ=255)
                                                                                                                                    gTCAGACCGATTGTCAACCGCATCCATCAGCAATCTGAACTGCGCGGCTGCTCTGAGTATTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGACCAATGTGATTTGACGAATTAGGGCt        >  8:319782/1‑139 (MQ=255)
                                                                                                                                          ccGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGATTTGACGAATTTGGGCTCCCGtt  <  1:14008/139‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                          
GTAAATCAAGTCATATAAACTGATAAATTCATTCCATTTTTCATGAAAGCTTGTGCTTCCATTCTTTTGGACATCTCGCTTCAGTTCAGTAAAAGCTTCCATAATCGGTTTTTCCAGTGAGCCCCAAAGCGGGTCAGACCGATTGTCAACCGCATCCATCAGCATTCTGAACTGCGCGGCTGCTCTGAGTTTTTGTGTGTCTGATGTATCAGCTTGCTTTAAACTCCGAACCATATCATTATAGCCCAATGTGATTTGACGAATTTGGGCTCCCGTT  >  BS168_reference_genome/2361844‑2362120

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: