Predicted mutation
evidence seq id position mutation annotation gene description
RA BS168_reference_genome 2,271,424 T→C R78R (AGA→AGG uvrX ← lesion bypass phage DNA polymerase; phage SPbeta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,271,4240TC98.9% 287.5 / NA 92R78R (AGA→AGGuvrXlesion bypass phage DNA polymerase; phage SPbeta
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (37/54);  minor base G (0/1);  total (37/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAA  >  BS168_reference_genome/2271291‑2271558
                                                                                                                                     |                                                                                                                                      
ttCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCtt                                                                                                                                     <  8:87566/137‑1 (MQ=255)
 tCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTgg                                                                                                                                   <  4:249985/138‑1 (MQ=255)
   cTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACATTGTGTATCCTTGGATc                                                                                                                                >  8:258685/1‑139 (MQ=255)
   cTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATc                                                                                                                                >  8:117810/1‑139 (MQ=255)
      cATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCtt                                                                                                                              <  8:240792/138‑1 (MQ=255)
          tAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCtt                                                                                                                              >  2:195025/1‑134 (MQ=255)
          tAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCa                                                                                                                            <  2:184323/136‑1 (MQ=255)
          tAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGt                                                                                                                         <  7:323158/139‑1 (MQ=255)
          tAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGt                                                                                                                         <  7:117810/139‑1 (MQ=255)
            aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCa                                                                                                                            <  1:401899/134‑1 (MQ=255)
            aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCa                                                                                                                            <  1:69922/134‑1 (MQ=255)
            aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCa                                                                                                                            <  2:333787/134‑1 (MQ=255)
            aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAg                                                                                                                           <  5:155253/135‑1 (MQ=255)
            aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGt                                                                                                                         <  3:291486/137‑1 (MQ=255)
            aaaaaGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTAt                                                                                                                       <  8:115346/139‑1 (MQ=255)
                         aaTTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAg                                                                                                              <  2:270655/135‑1 (MQ=255)
                              aaTACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATc                                                                                                      <  8:100161/138‑1 (MQ=255)
                               aTACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCt                                                                                                    >  3:164965/1‑139 (MQ=255)
                               aTACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCt                                                                                                    >  3:288780/1‑139 (MQ=255)
                                       ggACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAATAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGa                                                                                             <  5:385863/138‑1 (MQ=255)
                                        gACACAGTTCTCAGGAACAAATATGTAAAACAGCTTTGTAATTACAGATGAAACTCTCATGAAAAGCTTCATTAGTGGATTTACAAAGTCTATCATTGGATCTTCAGGTATCTCAAATAGACACGATCCTGTATTGAtt                                                                                           >  4:180203/1‑139 (MQ=11)
                                           acaTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGAtt                                                                                           >  1:332747/1‑136 (MQ=37)
                                             aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGTATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCaaa                                                                                       >  8:74629/1‑138 (MQ=255)
                                             aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  5:126672/1‑138 (MQ=37)
                                             aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  3:37259/1‑138 (MQ=37)
                                             aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  6:97200/1‑138 (MQ=37)
                                             aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  5:41567/1‑138 (MQ=37)
                                             aTTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  5:283243/1‑138 (MQ=37)
                                               ttCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaaaa                                                                                     >  8:165647/1‑138 (MQ=37)
                                               ttCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAAt                                                                                    >  8:127443/1‑139 (MQ=37)
                                                tctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCtt                                                                                                                              <  2:334465/96‑1 (MQ=255)
                                                tctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCtt                                                                                                                              >  1:334465/1‑96 (MQ=255)
                                                tctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGTTATCTCAAATAGTCTCGATCCTCTTTTGATTCCaa                                                                                       >  2:39809/1‑135 (MQ=255)
                                                tctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  2:299293/1‑135 (MQ=37)
                                                tctcAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCaa                                                                                       >  2:210339/1‑135 (MQ=37)
                                                      gAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCtttttt                                                                             <  5:31945/139‑1 (MQ=21)
                                                       aaCAAATCTGTAAATCATCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGAGTTACAATGTGTATCCTTGTATCTTCAGGTATCTCAAATAGTCTCGCTCCTGTTTTGATTCCAAAATCtttt                                                                               <  2:67686/136‑1 (MQ=11)
                                                       aaCAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGAGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTa                                                                            <  8:188568/139‑1 (MQ=21)
                                                           aaTCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTGTGATTCCAAATTCTTTTTTAAGTg                                                                        >  5:240256/1‑139 (MQ=255)
                                                                 tAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAg                                                                     >  7:37725/1‑136 (MQ=21)
                                                                   aaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGga                                                                   <  2:332747/136‑1 (MQ=21)
                                                                   aaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGga                                                                   <  2:204934/136‑1 (MQ=21)
                                                                   aaaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGAt                                                                <  5:194208/139‑1 (MQ=21)
                                                                    aaCAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGgag                                                                  <  1:138025/136‑1 (MQ=21)
                                                                    aaCAGCTTTGTAATTTCAGTTGAAAATCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGTAGATg                                                               <  6:277666/139‑1 (MQ=255)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGc                                                              <  2:353401/136‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  7:155738/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  7:121503/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  8:379947/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  3:69036/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  4:243156/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  5:285615/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  5:291264/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  4:21836/139‑1 (MQ=21)
                                                                        gCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  5:97200/139‑1 (MQ=21)
                                                                         cTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCa                                                             <  7:325147/136‑1 (MQ=21)
                                                                         cTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCa                                                             >  8:325147/1‑136 (MQ=21)
                                                                         cTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           >  4:69036/1‑138 (MQ=21)
                                                                           ttGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           <  4:186149/136‑1 (MQ=21)
                                                                           ttGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGc                                                           >  3:186149/1‑136 (MQ=21)
                                                                           ttGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTaa                                                        >  5:222038/1‑139 (MQ=21)
                                                                            tGTCATTTCAGGTGAAAATCTGATGAAAAGCGTAATTTGTGGAGTTCCACTGTGTCTCATTGGATCTTCCGGTCTCTCAAATAGTCTCTATCCTGTTTTGCTTCCAACATCTTTTTTAAGTGCAGGAGATGCCGCTaa                                                        <  3:145534/138‑1 (MQ=11)
                                                                            tGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAg                                                                                                              <  2:132662/84‑1 (MQ=255)
                                                                            tGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAg                                                                                                              >  1:132662/1‑84 (MQ=255)
                                                                              tAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTaata                                                      <  7:364299/138‑3 (MQ=21)
                                                                              tAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTaata                                                      <  7:165647/138‑3 (MQ=21)
                                                                              tAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTa                                                         <  1:210339/135‑1 (MQ=21)
                                                                              tAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGAGCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGAGCAGGAGATGCAGCTa                                                         <  2:379304/135‑1 (MQ=255)
                                                                                            aCTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGATAATACTACACTTCCCTGTc                                       <  6:328378/139‑1 (MQ=12)
                                                                                               ctGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCtct                                    >  5:338393/1‑137 (MQ=12)
                                                                                                tgatgaAAAGATTCATTTGTGGATTTACAATGTGTATCCTTGTATCTTCAGGTATCTCAAATAGTCTAGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGCTGCAGCTAATACTACACTTCCCTGTCtc                                     <  7:74629/137‑2 (MQ=11)
                                                                                                 gatgaAAAGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTcc                                  >  4:4003/1‑139 (MQ=12)
                                                                                                     aaaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTCTCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                               <  6:264213/138‑1 (MQ=255)
                                                                                                      aaaGCTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTAtt                             <  8:177174/139‑1 (MQ=12)
                                                                                                          cTTCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATAATACACTTCCCTGTCTCTCCGTATTccc                          >  6:299258/1‑138 (MQ=255)
                                                                                                            tCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                               >  3:270926/1‑131 (MQ=12)
                                                                                                            tCATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTa                               <  4:270926/131‑1 (MQ=12)
                                                                                                             cATTTGTGGATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGGGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTac                       <  5:87221/138‑1 (MQ=255)
                                                                                                                   tgGATTTACAATGTGTATCCTTGTATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTTCAGGAGATGCAGATAATACTACACTTCCCTGTCTCTCCGTATTCCCTacaaca                   <  1:377285/136‑1 (MQ=11)
                                                                                                                     gATTTACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCATGTCTCTCCGTATTCCCTACAACAGc                 >  5:215654/1‑136 (MQ=255)
                                                                                                                       tttACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTa                        >  2:76672/1‑127 (MQ=12)
                                                                                                                       tttACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTa                        <  1:76672/127‑1 (MQ=12)
                                                                                                                         tACAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAg          >  7:26367/1‑139 (MQ=255)
                                                                                                                           cAATGTGTATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGc         <  3:150290/138‑1 (MQ=255)
                                                                                                                                tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGtt    >  6:172062/1‑138 (MQ=255)
                                                                                                                                tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGtt    >  4:208090/1‑138 (MQ=255)
                                                                                                                                tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGtt    >  4:267193/1‑138 (MQ=255)
                                                                                                                                tgtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa   <  3:404912/139‑1 (MQ=255)
                                                                                                                                 gtATGCTTGGATATTCAGTTTTCTCAAATAGTCGCGATCTTGTTTTGATTCCAAATTCTTTTTTAAGTGAAGGTGTTTCAGATAATAGTACACTTACCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATTTTaa  <  6:41567/139‑1 (MQ=11)
                                                                                                                                 gtATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTa   >  5:186246/1‑138 (MQ=255)
                                                                                                                                  tATCCTTGGTTCTTCAGGTATCTCAAATTGTCTCTATCCTGTTTTGATTCCAAATTCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa  <  8:37725/138‑1 (MQ=11)
                                                                                                                                  tATCCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAATACTACACTTCCCTGTCTCTCCGTATTCCCTACAACAGCAAGATAGCATGTTaa  <  5:84256/138‑1 (MQ=255)
                                                                                                                                     |                                                                                                                                      
TTCCTGCATCTAAAAAAGATTCATCAATTGAATACGTATGGACACATTTCTCAGGAACAAATCTGTAAAACAGCTTTGTAATTTCAGTTGAAACTCTGATGAAAAGCTTCATTTGTGGATTTACAATGTGTATTCTTGGATCTTCAGGTATCTCAAATAGTCTCGATCCTGTTTTGATTCCAAAATCTTTTTTAAGTGCAGGAGATGCAGCTAACACTACACTTCCCTGTCTATCCGTATTCCCTACAACAGCAAGATAGCATGTTAA  >  BS168_reference_genome/2271291‑2271558

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: