Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 3,253,024 T→G 19.1% K142T (AAG→ACG)  comA ← two‑component response quorum‑sensing regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome3,253,0240TG19.1% 140.5 / 3.1 71K142T (AAG→ACG) comAtwo‑component response quorum‑sensing regulator
Reads supporting (aligned to +/- strand):  ref base T (20/36);  new base G (12/3);  total (32/39)
Fisher's exact test for biased strand distribution p-value = 3.11e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGTATTCAATGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCACCGTG  >  BS168_reference_genome/3252886‑3253148
                                                                                                                                          |                                                                                                                            
cTGTATTCAATGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCt                                                                                                                              <  7:170787/139‑1 (MQ=255)
   tATTCAATGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGttgg                                                                                                                           >  6:55295/1‑137 (MQ=255)
   tATTCAATGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTg                                                                                                                           >  4:126486/1‑139 (MQ=255)
     ttCAATGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGgttgtt                                                                                                                         >  3:296571/1‑139 (MQ=255)
         aTGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGtttttt                                                                                                                     >  8:48312/1‑139 (MQ=255)
          tGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGtttttt                                                                                                                     >  5:402743/1‑138 (MQ=255)
          tGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTg                                                                                                                    >  4:393080/1‑139 (MQ=255)
          tGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTg                                                                                                                    >  2:33982/1‑139 (MQ=255)
            gACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCt                                                                                                                  <  4:41091/139‑1 (MQ=255)
            gACCGCTTGCTTAAATGAAGGGCATCTGAGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCt                                                                                                                  <  4:389506/139‑1 (MQ=255)
             aCCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGGTtttttttttgtg                                                                                                                 >  2:274586/1‑135 (MQ=255)
                gcttgcttAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTTCTggat                                                                                                              >  6:135178/1‑136 (MQ=255)
                 cttgcttAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTTCTGAGt                                                                                                              >  3:7314/1‑138 (MQ=255)
                        aaaTGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAg                                                                                                        <  3:126486/137‑1 (MQ=255)
                        aaaTGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAg                                                                                                        <  7:411144/137‑1 (MQ=255)
                         aaTGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTg                                                                                                     >  4:380061/1‑139 (MQ=255)
                               gggCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTTTTTTTTTATGTGTCATCAGCTGTTTaa                                                                                                >  2:156576/1‑138 (MQ=255)
                                ggCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTTGTGGGTCATCAGCTGTTTAAAg                                                                                              >  4:328540/1‑139 (MQ=255)
                                ggCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTTCTGGGTCATCAGCGGGTTTAAg                                                                                              >  4:129621/1‑139 (MQ=255)
                                   aTCTGCGATTTCTTGGTTAGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACTTCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGATTTTTGCTGAGTCATCACCTGTTTATAGTa                                                                                            >  8:132060/1‑138 (MQ=255)
                                    tCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTTCTGAGTCATAAGCTGTTTAAAGTaat                                                                                          >  4:92776/1‑138 (MQ=255)
                                          aTTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGGTTGTTTTTTTGTGGGTAATAAGCGTTTTAAATTTAGAaaaa                                                                                     >  4:29281/1‑138 (MQ=255)
                                           tttCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATGAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAAt                                                                                    <  1:419763/138‑1 (MQ=255)
                                                          aTCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTcc                                                                    <  3:380061/139‑1 (MQ=255)
                                                          aTCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTcc                                                                    <  4:9218/139‑1 (MQ=255)
                                                                    aaCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTAc                                                           <  1:156576/138‑1 (MQ=255)
                                                                    aaCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACg                                                          <  1:33982/139‑1 (MQ=255)
                                                                    aaCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACg                                                          <  6:422406/139‑1 (MQ=255)
                                                                    aaCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACg                                                          <  7:304605/139‑1 (MQ=255)
                                                                    aaCTTCTTGAAGAATCAGGCATTCTCTAGGTGGGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTGTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACg                                                          <  6:190815/139‑1 (MQ=255)
                                                                     acagcttGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGt                                                         <  5:52298/135‑1 (MQ=255)
                                                                     aatgcttGAAGAATCATTCATTCTCTTGGTGTGAGCTCATCTTGTTATTTTTGAGTGTAAGGAGCAGGCTTTGTTTTTTTCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTAc                                                           <  1:172434/133‑1 (MQ=255)
                                                                     aCTTCTTGAAGAATCAGTCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGt                                                         <  6:402743/139‑1 (MQ=255)
                                                                     aCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGt                                                         <  5:135178/139‑1 (MQ=255)
                                                                     aCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGt                                                         <  4:7314/139‑1 (MQ=255)
                                                                     aCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGt                                                         <  3:328540/139‑1 (MQ=255)
                                                                     aCTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGAAAATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGt                                                         <  6:33900/139‑1 (MQ=255)
                                                                      cttcttGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGGCATGAGCAGTTTAAAGTAATCAAAATAGACTAAAATTTCTACGTTTACTACGTg                                                        >  4:213562/1‑139 (MQ=255)
                                                                            gaagaaTCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGta                                                     <  3:92776/136‑1 (MQ=255)
                                                                            gaagaaTCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGta                                                     <  7:48312/136‑1 (MQ=255)
                                                                            gaagaaTCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGta                                                     <  3:393080/136‑1 (MQ=255)
                                                                            gaagaaTCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGta                                                     <  2:336451/136‑1 (MQ=255)
                                                                               gaaTCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTa                                               <  3:29281/139‑1 (MQ=255)
                                                                               gaaTCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTa                                               <  3:129621/139‑1 (MQ=255)
                                                                                    aGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGGTTGTTTTTTTGTGAGTGATCATCTGTTTTAAGTTAGCAAAATAGACTTAAATTTCTGCGTTGTTTACGGGGTTTATGttagttg                                          >  6:52298/1‑133 (MQ=255)
                                                                                    aGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTAtttg                                           <  7:367250/138‑3 (MQ=255)
                                                                                      gacagctctAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTg                                        <  7:421495/134‑1 (MQ=255)
                                                                                        aTTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGAt                                      <  5:412312/139‑1 (MQ=255)
                                                                                          tctcTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGTTTTTTTTTTTATGTGTCATAATCTGTTTAAAGTTATCAAAATCTACTAAAATTTATCCGTTGAGTACGTGGTATATGTATTTGGTGAtct                                    >  2:384629/1‑139 (MQ=255)
                                                                                           ctctAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGAtctt                                   >  7:402681/1‑139 (MQ=255)
                                                                                                 gtgtGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGAtcttttcttt                             <  1:384629/139‑1 (MQ=255)
                                                                                                    tgAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGa                          <  5:378124/139‑1 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGGTTGTTTTTTTGGGGTTACTGATCGGGTTTAATTTAGCAAAATAGACAAAAATTTATCCGGTGATTTCCTGGTTTTTGTTTTTGGTGATGTtttatttatgtt                        >  2:353564/1‑135 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGGGTTGTTTTTTTCTGGGTCATCAGCTGTTTTAAGTTAGCAAAATAGACTAAAATTTCTCCGTTGAGTACGTGGTATTTGTTTTTGGTGTtcttttctttgtgtg                        >  8:421495/1‑134 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGGTGAGTCATCAGGGGTTTAAAGTAAGCAAAATCGACTAAAATTTCTGCGTTGAGTACGTGGTATATGTATTGGGTGGTGTTTTGTTTTGAtt                        >  6:342364/1‑139 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGAtt                        >  5:164074/1‑139 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACGAAAATTTATCCGGTGCGTACGTGGTATATGTATTTGGTGATCTTTTCTTTATAtt                        >  1:42026/1‑139 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTTTTTTTTTGTGAGTGATGACCAGTTTTAAGTTATCAAAAAACAATAAAAATTATCCGTTTAGTACGTGGTTTATTTTTTTGGTGGtgttttatttttgtt                        >  2:163051/1‑134 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTTCTGAGTCATCAGCTGTTTAAAGTTAGCAAAAAAGACTAAAATTTATACGGTGAGTACGTGGTATTTTTTTTGGGTGATCTTTTATTTAGAtt                        >  3:306981/1‑139 (MQ=255)
                                                                                                      aGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCGTTGTTTTTTGGTTAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACGAAAATTTTTACGTTGAGTTCGTGGTATATGTATTGGGTTATGTTTTGTTTTGAtt                        >  1:221769/1‑139 (MQ=255)
                                                                                                                gTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACg                                                          >  3:168246/1‑95 (MQ=255)
                                                                                                                gTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACg                                                          <  4:168246/95‑1 (MQ=255)
                                                                                                                gTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCt              >  5:151026/1‑139 (MQ=255)
                                                                                                                   cTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGa            <  6:151026/138‑1 (MQ=255)
                                                                                                                      tttGAGAGGAAGGAGCCGGCGTTGTTTTTTGCTGAGTCATCAGCGGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGc        >  1:223534/1‑139 (MQ=255)
                                                                                                                        tGAGATGAAGGAGGCGGCTTTCTTTTATGCTGAGTCATCAGCTGTTTAAAGAAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTATCTTTAGATTCCGTTTTGCTGATGGCa       <  7:132060/138‑1 (MQ=255)
                                                                                                                        tGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCa       <  3:177548/138‑1 (MQ=255)
                                                                                                                        tGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCa       <  2:223534/138‑1 (MQ=255)
                                                                                                                        tGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCa       <  1:353564/138‑1 (MQ=255)
                                                                                                                        tGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCa       <  5:55295/138‑1 (MQ=255)
                                                                                                                            aggaaggaGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCACCGTg  <  3:310654/139‑1 (MQ=255)
                                                                                                                              gaaggaGCCGGCTTTGTTTTTTTCTGAGTCATCACCTGTTTAAAGTAATCAAAATCGACTAAAATTTATGCGGTGAGTACGTGGTATTTTTATTTGGTGATATTTTATTTAGATTACGGTTTTCTGATGGCCCCGTg  >  1:323463/1‑137 (MQ=255)
                                                                                                                                      cGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCa       <  2:232017/124‑1 (MQ=255)
                                                                                                                                      cGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATAGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTTCTGATGGCa       >  1:232017/1‑124 (MQ=255)
                                                                                                                                          |                                                                                                                            
CTGTATTCAATGGACCGCTTGCTTAAATGAAGGGCATCTGCGATTTCTTGGTTTGTAAATCCCTTTTCAACTTCTTGAAGAATCAGGCATTCTCTAGGTGTGAGCACATCTTGTTCTTTTTGAGAGGAAGGAGCCGGCTTTGTTTTTTGCTGAGTCATCAGCTGTTTAAAGTAAGCAAAATCGACTAAAATTTCTCCGTTGAGTACGTGGTATATGTATTGGGTGATCTTTTCTTTAGATTCCGTTTTGCTGATGGCACCGTG  >  BS168_reference_genome/3252886‑3253148

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: