Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 892,022 A→T N196I (AAT→ATT)  malR → transcriptional activator of the Mal operon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome892,0220AT98.7% 266.7 / ‑4.8 76N196I (AAT→ATT) malRtranscriptional activator of the Mal operon
Reads supporting (aligned to +/- strand):  ref base A (0/1);  new base T (37/38);  total (37/39)
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.78e-01

ATTTAATTTGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGAATTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTC  >  BS168_reference_genome/891888‑892152
                                                                                                                                      |                                                                                                                                  
atttaatttGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGAtttt                                                                                                                                 >  3:137113/1‑138 (MQ=255)
 tttaatttGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTa                                                                                                                                >  4:56926/1‑138 (MQ=255)
            gcgcGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTTCGATTTTaaaaaa                                                                                                                           <  2:201632/132‑1 (MQ=255)
            gcgcGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTaaaaaa                                                                                                                           >  1:201632/1‑132 (MQ=255)
               cgATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTc                                                                                                                 >  7:17990/1‑139 (MQ=255)
                 atgatTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCt                                                                                                                  <  2:129949/136‑1 (MQ=255)
                  tgatTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGcc                                                                                                              <  5:178833/139‑1 (MQ=255)
                  tgatTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGcc                                                                                                              <  5:132009/139‑1 (MQ=255)
                          aattGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCATGCAAGAATGCTTTCTGCAAAGGGTATTCCATTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGAc                                                                                                      <  7:289110/136‑1 (MQ=255)
                          aaGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGAc                                                                                                      <  8:371083/139‑1 (MQ=255)
                               aaCGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCg                                                                                                  >  6:165806/1‑138 (MQ=255)
                               aaCGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCg                                                                                                  >  7:203847/1‑138 (MQ=255)
                                aCGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCa                                                                                                <  4:394264/139‑1 (MQ=255)
                                 cGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCAt                                                                                               <  8:267329/139‑1 (MQ=255)
                                        tttATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATa                                                                                              <  8:196269/133‑1 (MQ=255)
                                              atcatTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCg                                                                                   <  8:117373/138‑1 (MQ=255)
                                                catTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGc                                                                                    >  2:275159/1‑135 (MQ=255)
                                                catTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGAc                                                                                 >  8:57475/1‑138 (MQ=255)
                                                catTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATCCgcc                                                                                 >  7:326348/1‑136 (MQ=255)
                                                  tttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCa                                                                                                >  8:61708/1‑121 (MQ=255)
                                                  tttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCa                                                                                                <  7:61708/121‑1 (MQ=255)
                                                  tttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAg                                                                              >  4:356892/1‑139 (MQ=255)
                                                  tttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGACCAGCTGACTCCGCATAATCTTTATGCGACAAg                                                                              >  3:198626/1‑139 (MQ=255)
                                                   ttCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGAc                                                                                 >  2:337365/1‑135 (MQ=255)
                                                     ctctTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATTCGACAAGCa                                                                            <  5:80475/138‑1 (MQ=255)
                                                     ctctTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCa                                                                            <  8:203847/138‑1 (MQ=255)
                                                     ctctTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCa                                                                            <  6:140313/138‑1 (MQ=255)
                                                        tttCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCaa                                                                           <  4:348695/136‑1 (MQ=255)
                                                          tCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCtg                                                                      >  7:271693/1‑139 (MQ=255)
                                                          tCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCtg                                                                      >  7:55811/1‑139 (MQ=255)
                                                           ccGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAAc                                                                         <  1:337365/135‑1 (MQ=255)
                                                             ggAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGa                                                                   >  8:34179/1‑139 (MQ=255)
                                                                    cacCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAg                                                            >  7:49945/1‑139 (MQ=255)
                                                                       cGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAATCATTGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAg                                                            >  1:170791/1‑136 (MQ=255)
                                                                       cGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGa                                                         <  8:271693/139‑1 (MQ=255)
                                                                       cGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGa                                                         <  3:236851/139‑1 (MQ=255)
                                                                        ggAGCTGAAGCCTCAGGCAAGACTGCTTTCTGCAAAGGGTATTCCCATTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATACTCTTTATGCGACAAGCAAACCTGTGACATTAAGCGAt                                                        <  4:255420/139‑1 (MQ=255)
                                                                         gAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGa                                                         >  8:152724/1‑137 (MQ=255)
                                                                               aTTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACa                                                                 <  4:133665/123‑1 (MQ=255)
                                                                               aTTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACg                                                  >  4:221070/1‑138 (MQ=255)
                                                                               aTTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCAATTACGATTTTAAAAAAAAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACa                                                                 >  3:133665/1‑123 (MQ=255)
                                                                                             aaTGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCtttt                                   >  8:65885/1‑139 (MQ=255)
                                                                                               tGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGc                                 <  3:230011/139‑1 (MQ=255)
                                                                                                 cTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCtttt                                   <  2:144082/135‑1 (MQ=255)
                                                                                                 cTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCa                                <  7:57475/138‑1 (MQ=255)
                                                                                                  tttCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCAcc                              >  5:80618/1‑139 (MQ=255)
                                                                                                  tttCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCAcc                              >  7:57255/1‑139 (MQ=255)
                                                                                                    tCTGCAAAGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGt                            >  8:381255/1‑139 (MQ=255)
                                                                                                          aaGGGTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGtttt                         >  2:15678/1‑136 (MQ=255)
                                                                                                             ggTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCa                                                                                                >  8:396471/1‑62 (MQ=255)
                                                                                                             ggTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCa                                                                                                <  7:396471/62‑1 (MQ=255)
                                                                                                             ggTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACa                                                                 <  3:269154/93‑1 (MQ=255)
                                                                                                             ggTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCg                   <  5:165806/139‑1 (MQ=255)
                                                                                                             ggTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCg                   >  4:273833/1‑139 (MQ=255)
                                                                                                             ggTATTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACACCGCATAATCTTTATGCGACAAGCAAACCTGTGACa                                                                 >  4:269154/1‑93 (MQ=255)
                                                                                                                aTTCCCTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTg                >  7:163943/1‑139 (MQ=255)
                                                                                                                   cccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGGTTTTTCTCGTTg                >  1:215889/1‑136 (MQ=255)
                                                                                                                   cccTTTATTTCCATTACGATTTTAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg             <  3:56926/138‑1 (MQ=255)
                                                                                                                    ccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTg                <  1:80592/135‑1 (MQ=255)
                                                                                                                    ccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg             <  3:273833/138‑1 (MQ=255)
                                                                                                                    ccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg             <  7:381255/138‑1 (MQ=255)
                                                                                                                    ccTTTATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg             <  8:163943/138‑1 (MQ=255)
                                                                                                                       ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGc                                 <  7:407790/115‑1 (MQ=255)
                                                                                                                       ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGc                                 >  8:407790/1‑115 (MQ=255)
                                                                                                                       ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCAc                               >  5:250576/1‑117 (MQ=255)
                                                                                                                       ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCAc                               <  6:250576/117‑1 (MQ=255)
                                                                                                                       ttATTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg          >  4:347973/1‑138 (MQ=255)
                                                                                                                         aTTTCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCa        <  4:137113/138‑1 (MQ=255)
                                                                                                                          tttCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAg          >  2:340766/1‑135 (MQ=255)
                                                                                                                            tCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg             <  3:173115/130‑1 (MQ=255)
                                                                                                                            tCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCg             >  4:173115/1‑130 (MQ=255)
                                                                                                                            tCCATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTg     <  4:198626/138‑1 (MQ=255)
                                                                                                                              ccttACGATTTTAAAAAACAATGGCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTg     <  2:170791/134‑1 (MQ=255)
                                                                                                                              cATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTc  >  6:371561/1‑139 (MQ=255)
                                                                                                                              cATTACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTc  <  8:49945/139‑1 (MQ=255)
                                                                                                                                ttACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGtt   <  2:48947/136‑1 (MQ=255)
                                                                                                                                 tACGATTTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTc  >  2:374366/1‑136 (MQ=255)
                                                                                                                                      |                                                                                                                                  
ATTTAATTTGATGCGCGATGATTTTAAAGTGAACGATCTGTTTATCATCATTTCTCTTTCCGGAGAAACACCGGAGCTGATTCCTCAGGCAAGAATGCTTTCTGCAAAGGGTATTCCCTTTATTTCCATTACGAATTTAAAAAACAATGTCCTCGCCCAGCTGACTCCGCATAATCTTTATGCGACAAGCAAACCTGTGACATTAAGCGATAGAACGGAAATTGTCGCTTTTGCACCGTTTTTTCTCGTTGGCGAAGCATTGTTC  >  BS168_reference_genome/891888‑892152

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: