Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 3,935,825 Δ1 bp coding (1190/1191 nt) rlmI ← large subunit ribosomal RNA m5C1962 methyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,935,8230T.99.3% 1353.7 / 2.1 273intergenic (+35/+1)ywbE/rlmIconserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (2/0);  new base . (136/134);  total (139/134)
Fisher's exact test for biased strand distribution p-value = 4.98e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 285 total aligned reads displayed.

GACGAGAGGGACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAGCCTTT  >  BS168_reference_genome/3935674‑3935967
                                                                                                                                                     |                                                                                                                                                
 aCGAGAGGGACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCatt                                                                                                                                                 <  1:196353/150‑1 (MQ=255)
   gagaGGGACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCGTTTTCattat                                                                                                                                               >  2:2873065/1‑148 (MQ=255)
    agagGGACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGGTGCCCCTGAACAGGGCAGCTTTTTCattatt                                                                                                                                              >  1:2682029/1‑147 (MQ=255)
    agagGGACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCattatt                                                                                                                                              >  1:1018642/1‑147 (MQ=255)
       gggACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAttttt                                                                                                                                          <  1:184667/150‑1 (MQ=255)
          aCAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGATGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCt                                                                                                                                       >  1:2354410/1‑150 (MQ=255)
           cAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAACGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTg                                                                                                                                      <  1:3885654/150‑1 (MQ=255)
               gAAGGATATTTTGACGAAATCAACCTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGa                                                                                                                                  <  1:4303536/150‑1 (MQ=255)
                aaGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCCTTTTCATTATTTTTTCTGCAGaa                                                                                                                                 >  1:1063559/1‑150 (MQ=255)
                 aGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAg                                                                                                                                >  2:4367003/1‑150 (MQ=255)
                   gATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGATGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTGTTCATTATTTTTTCTGCAGAAGGa                                                                                                                              >  1:4404647/1‑150 (MQ=255)
                       ttttGACGAAACCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCAc                                                                                                                          <  2:491562/150‑1 (MQ=255)
                        tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACt                                                                                                                         <  2:2011653/150‑1 (MQ=255)
                        tttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACt                                                                                                                         <  2:87257/150‑1 (MQ=255)
                         ttGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAACTGCCCCCGAACAGGGCAGCGTTTTCATTATTTTTTCTGCAGAAGGTCCACtt                                                                                                                        >  1:3706385/1‑150 (MQ=255)
                            tcGAAATCAACCGTCCACCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCa                                                                                                                     <  2:1353075/149‑1 (MQ=255)
                            aCGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATTAAAAAAAGCTGCCCCTGAACAGGGCAGGCTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCa                                                                                                                     >  1:4151461/1‑150 (MQ=255)
                            aCGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCa                                                                                                                     >  2:1239766/1‑150 (MQ=255)
                             cGAAATCAAACTTTCATCCGCATGGCATTAAAGCCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAg                                                                                                                    <  1:1613919/150‑1 (MQ=255)
                              gAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCTGGAGGACCACATTCAgg                                                                                                                   >  1:3590868/1‑150 (MQ=255)
                                 aTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCTGAAGGACCACTTTCAGGTAg                                                                                                                >  2:3798627/1‑150 (MQ=255)
                                 aTCAAACTTCCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAg                                                                                                                <  2:698536/150‑1 (MQ=255)
                                   cAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGGCCACTTTCAGGTAGtt                                                                                                              >  2:2792596/1‑150 (MQ=255)
                                       cTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCt                                                                                                          >  2:4571939/1‑150 (MQ=255)
                                         ttCATCCGCATGGCATTAAAGTCCGTCTTGAAGATCGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTc                                                                                                        <  2:2626986/150‑1 (MQ=255)
                                          tCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTcc                                                                                                       <  2:539324/150‑1 (MQ=255)
                                            aTCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCgg                                                                                                     <  2:4046597/150‑1 (MQ=255)
                                                  aTGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGc                                                                                               >  2:3499359/1‑150 (MQ=255)
                                                  aTGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGc                                                                                               >  2:4429128/1‑150 (MQ=255)
                                                  aTGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGCTTCTACATAAAAAAAAGCAGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGc                                                                                               <  1:2956648/150‑1 (MQ=255)
                                                      cATTAAGGCCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaa                                                                                           <  2:330970/150‑1 (MQ=255)
                                                      cATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaa                                                                                           >  1:1920253/1‑150 (MQ=255)
                                                       aTTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGaaa                                                                                          >  2:830868/1‑150 (MQ=255)
                                                            aGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGGAATGCCGAAATCGtt                                                                                     >  1:2099303/1‑150 (MQ=255)
                                                             gtccgtccTGAGGATGGCCGCATCGGCCGCGACCAGCAGATCGTGTCTTCATAAAAAAAAGCCGCCCCTGAGCAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGttt                                                                                    <  2:3010539/150‑1 (MQ=255)
                                                              tCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGtttt                                                                                   <  2:173565/150‑1 (MQ=255)
                                                                cGTCTTGAAGATGGCCGCACGGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                                 <  2:3664897/150‑1 (MQ=255)
                                                                   cTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAt                                                                              <  2:1617709/150‑1 (MQ=255)
                                                                    ttGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCCTTTTGAAAtc                                                                             <  1:3476824/150‑1 (MQ=255)
                                                                      gAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcct                                                                           >  1:4292070/1‑150 (MQ=255)
                                                                       aaGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCAACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcctc                                                                          <  1:3920507/150‑1 (MQ=255)
                                                                        aGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTTCATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAtcctcc                                                                         >  1:864328/1‑150 (MQ=255)
                                                                           tGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                                      >  1:2525410/1‑150 (MQ=255)
                                                                            ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCa                                                                     >  2:4683618/1‑150 (MQ=255)
                                                                              ccGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAGGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGc                                                                   <  1:2620015/150‑1 (MQ=255)
                                                                              ccGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGc                                                                   <  2:236257/150‑1 (MQ=255)
                                                                                gCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGGACAGGGCAGCTTTGTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGt                                                                 >  2:1117074/1‑150 (MQ=255)
                                                                                  aTCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAACGGaa                                                               >  2:1397811/1‑150 (MQ=255)
                                                                                   tCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTaaa                                                              >  2:618961/1‑150 (MQ=255)
                                                                                     ggCCGCGCCCAGCAGATCGTTTCTCCATAAAAAACAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAAtt                                                            <  2:1658851/150‑1 (MQ=255)
                                                                                      gCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGATGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                                           >  1:3970198/1‑150 (MQ=255)
                                                                                        cgcgACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGGAAATTCCt                                                         >  2:2751106/1‑150 (MQ=255)
                                                                                         gcgACCAGCAGATCGTTTCTACATCAAAACAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTc                                                        <  2:2655622/150‑1 (MQ=255)
                                                                                          cgACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCa                                                       <  2:1252272/150‑1 (MQ=255)
                                                                                            cccAGCAGATCGTTCCTACATAAAAAGATGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAAt                                                     <  2:228025/149‑1 (MQ=255)
                                                                                               agcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGAt                                                  >  2:2445192/1‑150 (MQ=255)
                                                                                               agcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAAGGCCGAAATCGGTTTGGAGTCCTCCGGCAGCGGACATTCCTCAATGAt                                                  >  1:1829295/1‑150 (MQ=255)
                                                                                               agcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCCTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGCAATGCCGAAATCGTTTTGAAATCCTCCGGCCGCGTAAATTCCTCAATGAt                                                  >  1:569791/1‑150 (MQ=255)
                                                                                                gcagACCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATc                                                 <  2:2606128/150‑1 (MQ=255)
                                                                                                  agATGGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgt                                               <  2:1818304/150‑1 (MQ=255)
                                                                                                    aTCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGGAAATTCCTCAATGATCgtgt                                             >  2:3175827/1‑150 (MQ=255)
                                                                                                     tCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                                            >  1:4515840/1‑150 (MQ=255)
                                                                                                     tCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                                            <  1:785794/150‑1 (MQ=255)
                                                                                                       ttttCTCCATAAAAAAAGGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAATTGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGc                                          <  2:2092098/149‑1 (MQ=255)
                                                                                                        tttCTCCATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCg                                         <  2:3863460/150‑1 (MQ=255)
                                                                                                         ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAGATCCTCCGGCAGCGTAAATTCCTCCATCATCGTGTAGCGt                                        >  1:4839741/1‑150 (MQ=255)
                                                                                                         ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGt                                        <  1:2894465/150‑1 (MQ=255)
                                                                                                         ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGt                                        <  2:1063559/150‑1 (MQ=255)
                                                                                                         ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGt                                        >  1:3071662/1‑150 (MQ=255)
                                                                                                            tACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCTTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGTAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCa                                     >  1:4239524/1‑150 (MQ=255)
                                                                                                                 aaaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                                >  2:564386/1‑150 (MQ=255)
                                                                                                                  aaaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGtt                               <  2:4082754/150‑1 (MQ=255)
                                                                                                                     aaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCt                            >  1:2973270/1‑150 (MQ=255)
                                                                                                                       aaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCttt                          >  1:3059565/1‑150 (MQ=255)
                                                                                                                       aaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCttt                          <  2:2830042/150‑1 (MQ=255)
                                                                                                                       aaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCttc                          >  1:959457/1‑149 (MQ=255)
                                                                                                                           tGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAg                      <  2:792287/150‑1 (MQ=255)
                                                                                                                            gCCCCTGCACAGGGCAGCTTTATCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAgg                     >  1:3772681/1‑150 (MQ=255)
                                                                                                                             ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGGAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGGAGCGTTCATTTGTTTTTGTAAAggc                    >  2:1784293/1‑150 (MQ=255)
                                                                                                                               ccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAGCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCCGCAGCGGAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg                  >  2:4234320/1‑150 (MQ=255)
                                                                                                                               ccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcgg                  <  2:4298679/150‑1 (MQ=255)
                                                                                                                                   aaCAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATc              >  2:2194949/1‑150 (MQ=255)
                                                                                                                                    aCAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGTCCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCa             >  1:1534828/1‑150 (MQ=255)
                                                                                                                                     cAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCaa            >  1:4703955/1‑150 (MQ=255)
                                                                                                                                      aGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAAt           >  1:2282293/1‑150 (MQ=255)
                                                                                                                                      aGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTGTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGAATCaaa           >  2:1495828/1‑149 (MQ=255)
                                                                                                                                        ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAaaaa         >  2:2583030/1‑150 (MQ=255)
                                                                                                                                        ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATaa         <  2:3458756/150‑1 (MQ=255)
                                                                                                                                        ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATaa         <  1:1462101/150‑1 (MQ=255)
                                                                                                                                        ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGAATCAATaa         >  1:3508668/1‑150 (MQ=255)
                                                                                                                                        ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTGTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATaa         >  2:142741/1‑150 (MQ=255)
                                                                                                                                         gCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATaaa        <  1:2797797/150‑1 (MQ=255)
                                                                                                                                         gCAGCTCTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATaaa        <  2:2682029/150‑1 (MQ=255)
                                                                                                                                          cAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGGTCGTGGAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAg       >  1:638153/1‑150 (MQ=255)
                                                                                                                                           aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAg       <  1:1286311/149‑1 (MQ=255)
                                                                                                                                            gCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGGAAATTCCTCAATGATCGTGAAGCGTTCATCTGTCTCTTTAAGGGCGGCATCAATAAAGcc     >  2:3249918/1‑150 (MQ=255)
                                                                                                                                             cTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAGCCt    >  2:3431137/1‑150 (MQ=255)
                                                                                                                                             cTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCCTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAGCCt    >  2:3957620/1‑150 (MQ=255)
                                                                                                                                              tttttCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGATGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAGCCtt   <  1:3344474/150‑1 (MQ=255)
                                                                                                                                               ttttCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAGCCttt  <  2:954745/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GACGAGAGGGACAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCATCAATAAAGCCTTT  >  BS168_reference_genome/3935674‑3935967

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: