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.98.9% 396.8 / ‑2.0 90intergenic (+35/+1)ywbE/rlmIconserved protein of unknown function/large subunit ribosomal RNA m5C1962 methyltransferase
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base . (46/42);  total (48/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

CAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCAT  >  BS168_reference_genome/3935685‑3935954
                                                                                                                                           |                                                                                                                                   
cAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGCCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCat                                                                                                                                     >  2:205555/1‑139 (MQ=255)
cAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGCACAGGGCAGCTTTTTGat                                                                                                                                     >  5:324909/1‑139 (MQ=255)
cAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCat                                                                                                                                     >  3:93434/1‑139 (MQ=255)
     aaGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTGTTCattatt                                                                                                                                 >  7:415869/1‑136 (MQ=255)
      agtatatTTTCACCAAATCAAACTTTCCTCCCCATGGCATAAAAGCCCGGCTTGACGATGGCCGCAGGGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAtttt                                                                                                                              <  4:94935/136‑1 (MQ=255)
      agcatatTTTCACCAAATCAAACTTTCAGCCGCAGGGCATAAACGGCCGTCTTTACGATGCCCGCAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAtttt                                                                                                                              <  2:10223/136‑1 (MQ=255)
       ggATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCAGCCGCGACCAGCACATCGTTTCTACATAAAAAAAATCTGCCCCTCAACAGGGCAGCTTTTTAATTAttttt                                                                                                                             >  8:405792/1‑139 (MQ=255)
        gATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGCTCATTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAtttttt                                                                                                                            >  7:97488/1‑139 (MQ=255)
        gATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAtttttt                                                                                                                            >  3:531707/1‑139 (MQ=255)
        gATATTTTGACGAAATCAAACTTTAATCCGCATGGCATTAATGTCGGTCTTGACGAGGGCCGCTTGGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTAtttttt                                                                                                                            <  6:170897/139‑1 (MQ=255)
           atTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTg                                                                                                                         >  1:482493/1‑139 (MQ=255)
                    aatcaaaaaTTTCAGCCGCTTGCCATTAATGTCGGTTTTTAATTTGGCGGCTTCGCCCGCTACCAGCAGATAGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  2:482493/134‑1 (MQ=255)
                    aaagaaaaCTTTCATCCGCAGGGCATAAAAGTCCGTCGGGACGATGCCCCCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGTACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  2:309274/134‑1 (MQ=255)
                    aaaaaaaaTTTTAATCCTAATGCCATTAAAGTCGTTTTTCACTAGGGCCGCTGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  5:545195/134‑1 (MQ=255)
                    aaaaaaaaCTTTCATCCGCTTGGAATAAAAGGCGGTCTTCACGATGGCCGAAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  4:531707/134‑1 (MQ=255)
                    aaaaCAAACTTTCATCCGAATGGCTTTAAAGTCCGTCTTTAAGTTGGCCGCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  3:459737/139‑1 (MQ=255)
                    aaaTCACACTTCCACCCCCATGCAATTAAAGCCCGGCTTGAACATGCCCCCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  4:442396/139‑1 (MQ=255)
                    aaaTCAAACTTTCATCCGCATGGCATTAAAGGCCGTCTTGAAGATGGCCGCATCAGCCGCGACCAGCATATCTTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  7:405792/139‑1 (MQ=255)
                    aaaTCAAACTTTCAACCCCATGCCATAAACGCCCTCCTGCACGAGGCCCCCATCGCCCGCGACCAGCAGATCGTTTCTACAAAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  4:450408/139‑1 (MQ=255)
                    aaaTCAAACTTACATCCGAAGGGCATAAACGGCCGTTTTAAAGAGGCCCGCTGGGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  8:370497/139‑1 (MQ=255)
                    aaaCCAAACTTTCATCCGCTTGCAATTAACGTCCGTCTTGACTTGGGCCGCAGCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  5:130070/139‑1 (MQ=255)
                    aaaCCAAACTTCCATCCCCAGGCCATAAACGCCCGTCTTGACGAGGCCCGCCGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  5:468090/139‑1 (MQ=255)
                     aaTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                >  2:93622/1‑138 (MQ=255)
                      aTCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGCCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCAGCACAGGGCAGCTGTTTCATTATTTTTTCTGCAGAAGGAccc                                                                                                              >  5:34915/1‑138 (MQ=255)
                         aaaaTTTCATCCGCCTGGCATTAATGTCCGTCTTGACGTGGGCCGCAGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACtt                                                                                                           <  8:275655/139‑1 (MQ=255)
                           atttgaaTCCCCTTGGAATTAACGTCGGTTTTGATGTGGGCCGCCGCGCCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTc                                                                                                         <  2:393631/133‑1 (MQ=255)
                           aCTTTCATCCGCATGGCATTAACGTCCGTCTTTATGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTc                                                                                                         <  6:386031/139‑1 (MQ=255)
                                 aTCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCTCGACCAGCACATCGTTTCTACATAAAAAAAAGCTGCCACTGTACACGGGAGCTTTTTCATTATATTTTATGCAGAAGTCCCCCTTTAAGGTAg                                                                                                   >  1:370489/1‑139 (MQ=255)
                                      cATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGcc                                                                                              >  7:431028/1‑139 (MQ=255)
                                                aaGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAAt                                                                                    >  4:181225/1‑139 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATTGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                >  3:178915/1‑106 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATTGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAc                                                                                                                <  4:178915/106‑1 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACTTAAAAAAAAACTGACCCTGAACAGGGCAGCTTTTTGATTATTTTTTCTGCAGAAGGACCCCTATCAGGTAGATGCCTTACGGAAATGACg                                                                                >  7:342971/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                                >  5:138367/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                                >  1:303628/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                                >  4:268408/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                                >  7:30431/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTATTCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCg                                                                                >  6:58746/1‑138 (MQ=255)
                                                     cGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTAATTTTTTTTTATGCAGAAGGACCACTTTCAGGTATTTGCCTTCCGGAAATGCCg                                                                                >  7:215122/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTGATTATTTTTTATGCAGAAGGACCACTTTAAGGTATTTGCCTTCCGGAAATGCCGAAAAGGTTTTGa                                                                    >  2:269210/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                    >  4:224649/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                    >  6:45483/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                    >  7:213430/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                    >  3:161450/1‑138 (MQ=255)
                                                                 ggCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCACGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGa                                                                    >  4:304974/1‑138 (MQ=255)
                                                                   ccGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGTACCAATTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAAt                                                                 >  2:368374/1‑139 (MQ=255)
                                                                          ggCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCgg                                                          <  8:431028/139‑1 (MQ=255)
                                                                          ggCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCgg                                                          <  3:181225/139‑1 (MQ=255)
                                                                           gCCGCGACCAGCAGATCTTTTCTACATAAAAAAAAGCTGCCGCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTTGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATGCTCCGGc                                                         >  6:259638/1‑139 (MQ=255)
                                                                           cccGCGACCAGCATATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGc                                                         <  1:269210/138‑1 (MQ=255)
                                                                            ccGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCa                                                        <  3:3690/139‑1 (MQ=255)
                                                                                 aCCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGt                                                    >  8:170058/1‑138 (MQ=255)
                                                                                   cagaagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTaaa                                                 <  8:502529/139‑1 (MQ=255)
                                                                                     gcagATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGACATCCTCCGGCAGCGTAAAtt                                               >  5:102052/1‑139 (MQ=255)
                                                                                      cagCTCGTTTCTACATAAAAAAACGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTc                                              <  4:56571/139‑1 (MQ=255)
                                                                                         aTCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTc                                           <  8:213430/139‑1 (MQ=255)
                                                                                           cgattCTAAAAAATAAAAACATGACACTGATCAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGAACACTTTCAGGTAGTTGCCTTACGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                         <  2:227319/136‑1 (MQ=255)
                                                                                            gTTTGTAAATAAATAAAAGCTGCCCCTGAAGAGGGCAGCTTTTTGATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGGCTTCCGGTAATGCCGATATCGTTTTGTAATCCTCCGGCAGCGTAAATTCCTCaa                                         <  4:42925/138‑1 (MQ=255)
                                                                                            gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                         <  2:303628/138‑1 (MQ=255)
                                                                                            gTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAAt                                        >  7:521293/1‑139 (MQ=255)
                                                                                            gTTTCTAAATTAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATTCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                         <  6:517800/138‑1 (MQ=255)
                                                                                             tttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAAt                                        >  8:262908/1‑138 (MQ=255)
                                                                                              ttatACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                         <  2:531245/133‑1 (MQ=255)
                                                                                              ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAATTTCCTCAATGa                                      >  2:354848/1‑139 (MQ=255)
                                                                                              ttCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAACTTCCTCAATGa                                      >  8:419258/1‑139 (MQ=255)
                                                                                               tCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGa                                      <  4:557240/138‑1 (MQ=255)
                                                                                                cTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATc                                    >  4:104674/1‑139 (MQ=255)
                                                                                                   cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGACATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                 >  7:30410/1‑139 (MQ=255)
                                                                                                   cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                 <  6:34915/139‑1 (MQ=255)
                                                                                                   cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                 <  8:215122/139‑1 (MQ=255)
                                                                                                   cATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCgtg                                 <  3:268408/139‑1 (MQ=255)
                                                                                                     tAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTa                               <  7:170058/139‑1 (MQ=255)
                                                                                                         aaaaaaGCTGCCCCTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCg                            >  6:525528/1‑139 (MQ=255)
                                                                                                             caGCTGCCCCTGTACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCa                        <  2:370489/138‑1 (MQ=255)
                                                                                                                   ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                   <  4:161450/138‑1 (MQ=255)
                                                                                                                   ccccTGAACAGGGCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                   <  1:499015/138‑1 (MQ=255)
                                                                                                                   ccccTGAACAGGGAAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGt                   <  4:300814/138‑1 (MQ=255)
                                                                                                                             gggCAGCTTTTTCATTATTTTTTATGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTACACAATGAACGTGTAGGGTTCATTTGTATCTTTAAAgg        >  2:107564/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                         <  3:327253/105‑1 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCaa                                         >  4:327253/1‑105 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  1:395989/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  8:283436/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  3:419043/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  3:65320/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  3:531495/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  2:115010/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  6:469167/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  5:426389/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggc       >  7:102217/1‑139 (MQ=255)
                                                                                                                              ggCAGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCAGTTGTTTCTTTAACggc       >  7:170408/1‑139 (MQ=255)
                                                                                                                                 aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc    <  8:102217/139‑1 (MQ=255)
                                                                                                                                 aGCTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAggcggc    <  3:280321/139‑1 (MQ=255)
                                                                                                                                   cTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGGTCATTTGTTTCTTTAAAGGCGGCAt  >  7:484004/1‑139 (MQ=255)
                                                                                                                                   cTTTTTCATTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATAGTGTAGCGTTGATTTGTTTTTTTTAAagcgggaa  >  2:531147/1‑136 (MQ=255)
                                                                                                                                           |                                                                                                                                   
CAGTGAAGGATATTTTGACGAAATCAAACTTTCATCCGCATGGCATTAAAGTCCGTCTTGAAGATGGCCGCATCGGCCGCGACCAGCAGATCGTTTCTACATAAAAAAAAGCTGCCCCTGAACAGGGCAGCTTTTTCATTTATTTTTTCTGCAGAAGGACCACTTTCAGGTAGTTGCCTTCCGGAAATGCCGAAATCGTTTTGAAATCCTCCGGCAGCGTAAATTCCTCAATGATCGTGTAGCGTTCATTTGTTTCTTTAAAGGCGGCAT  >  BS168_reference_genome/3935685‑3935954

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