Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,126,886 G→A 9.1% F189F (TTC→TTT ntdC ← biosynthesis of neotrehalosadiamine (3,3'‑diamino‑3,3'‑dideoxy‑alpha,beta‑trehalose), glucose‑6‑P 3‑dehydrogenase; kanosamine biosynthesis (recent HGT island)

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,126,8860GA 9.1% 201.6 / 7.6 66F189F (TTC→TTTntdCbiosynthesis of neotrehalosadiamine (3,3'‑diamino‑3,3'‑dideoxy‑alpha,beta‑trehalose), glucose‑6‑P 3‑dehydrogenase; kanosamine biosynthesis (recent HGT island)
Reads supporting (aligned to +/- strand):  ref base G (25/35);  new base A (3/3);  total (28/38)
Fisher's exact test for biased strand distribution p-value = 6.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATGATATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAGTGTACTTTTATG  >  BS168_reference_genome/1126749‑1127024
                                                                                                                                         |                                                                                                                                          
atgatATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCGTTACTGATTCTGTGATGAAATCACTATCg                                                                                                                                            >  3:191339/1‑138 (MQ=255)
atgatATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCg                                                                                                                                            >  6:119786/1‑138 (MQ=255)
atgatATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCg                                                                                                                                            <  4:267386/138‑1 (MQ=255)
atgatATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCg                                                                                                                                            >  7:250449/1‑138 (MQ=255)
atgatATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCg                                                                                                                                            >  3:148993/1‑138 (MQ=255)
   atatTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAAc                                                                                                                                         <  3:168793/138‑1 (MQ=255)
   atatTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACa                                                                                                                                        <  7:18888/139‑1 (MQ=255)
        tACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAatac                                                                                                                                   >  3:108991/1‑139 (MQ=255)
                ccTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAAtaccatacca                                                                                                                            <  5:164067/138‑1 (MQ=255)
                ccTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAAtaccatacca                                                                                                                            <  5:146636/138‑1 (MQ=255)
                ccTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAAtaccatacca                                                                                                                            <  4:18904/138‑1 (MQ=255)
                  ttGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGACCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATg                                                                                                                          >  3:40687/1‑138 (MQ=255)
                   tGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGAGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGCACAAATACCATACCATGt                                                                                                                         >  4:179118/1‑138 (MQ=255)
                     ttcttcAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAAc                                                                                                                                         >  8:48128/1‑120 (MQ=255)
                     ttcttcAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAAc                                                                                                                                         <  7:48128/120‑1 (MQ=255)
                     ttcttcAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATAACTCTCGAACAAATACCATACCATGTCa                                                                                                                       >  3:66637/1‑138 (MQ=255)
                     ttcttcAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAAt                                                                                                                     <  2:208295/139‑1 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAAtaccatacca                                                                                                                            <  2:54739/128‑1 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAAtaccatacca                                                                                                                            >  1:54739/1‑128 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATaccataccat                                                                                                                           <  6:162030/129‑1 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATaccataccat                                                                                                                           >  5:162030/1‑129 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATaa                                                                                                                   <  2:4855/137‑1 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATaa                                                                                                                   <  1:219252/137‑1 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAgg                                                                                                                 >  5:57959/1‑139 (MQ=255)
                          cAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCATCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCAATCTCGAACAAATACCATACCATGTCAATAAgt                                                                                                                 >  5:168350/1‑138 (MQ=255)
                                gCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATa                                                                                                           <  6:147271/139‑1 (MQ=255)
                                   aTAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATAccc                                                                                                        <  7:235436/139‑1 (MQ=255)
                                    tAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCAAACAAATACCATACCATGTCAATAAGGTGAATACCCa                                                                                                       >  6:10431/1‑139 (MQ=255)
                                    tAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCAAACAAATACCATACCATGTCAATAAGGTGAATACCCa                                                                                                       >  8:34612/1‑139 (MQ=255)
                                      aaTTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGa                                                                                                     <  5:119786/139‑1 (MQ=255)
                                            tGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATcccc                                                                                                >  2:202779/1‑138 (MQ=255)
                                              cATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCAAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCaa                                                                                             <  5:10431/139‑1 (MQ=255)
                                               aTAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCAAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCaa                                                                                             <  7:34612/138‑1 (MQ=255)
                                                   tCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATc                                                                                        <  5:24910/139‑1 (MQ=255)
                                                   tCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATc                                                                                        <  4:61513/139‑1 (MQ=255)
                                                     aTCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAATTCACTCTCGAACAAATACCATACCATGTCAATAATGTTAATACCCAGATCCCCGAATGATCCa                                                                                      <  7:169063/139‑1 (MQ=255)
                                                        tACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCAAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTc                                                                                   >  8:168933/1‑139 (MQ=255)
                                                             gTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATc                                                                                        >  5:82851/1‑129 (MQ=255)
                                                             gTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATc                                                                                        <  6:82851/129‑1 (MQ=255)
                                                             gTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGtctt                                                                              >  2:268447/1‑139 (MQ=255)
                                                             gTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGgct                                                                               >  5:239015/1‑138 (MQ=255)
                                                                ccTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGtcttct                                                                            <  8:191571/138‑1 (MQ=255)
                                                                ccTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCAAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGtcttct                                                                            <  7:168933/138‑1 (MQ=255)
                                                                ccTGTTTATCCTCTTTAGATTTTACATTTGTGTTCCTCTTTGCCCTTACTGAATCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGtcttct                                                                            <  8:267274/138‑1 (MQ=255)
                                                                       aTCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTAt                                                                      >  1:265697/1‑137 (MQ=255)
                                                                        tCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAg                                                                   >  4:289512/1‑139 (MQ=255)
                                                                          ctctTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGc                                                                  <  1:202779/138‑1 (MQ=255)
                                                                               ttGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTc                                                               <  2:265697/136‑1 (MQ=255)
                                                                               ttGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCAt                                                             <  4:66637/138‑1 (MQ=255)
                                                                               ttGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCAt                                                             <  4:162192/138‑1 (MQ=255)
                                                                                    tttACATTTGTGTTCATATTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACACAGATCCCCCAATGATCCACTCGTCTTCTTAATATTAGGGTCAtctaggc                                                       >  2:246029/1‑135 (MQ=255)
                                                                                            tgtgTTCATATTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTaaa                                                <  7:244635/138‑1 (MQ=255)
                                                                                            tgtgATCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGACTACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTaaa                                                <  7:244631/138‑1 (MQ=255)
                                                                                            tgtcTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAAAACCATACCATGTCAATAATGTGAATACCCAGATCCCCCAATGATCCACTCCACATCTTACATTTACCGTCATCACACCaacaaaa                                                >  7:154928/1‑133 (MQ=255)
                                                                                                        tGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCt                                    <  4:105113/138‑1 (MQ=255)
                                                                                                            cTTTCTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTg                               <  4:40687/139‑1 (MQ=255)
                                                                                                            cTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTg                               <  8:250449/139‑1 (MQ=255)
                                                                                                               aCTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCag                            >  1:37760/1‑138 (MQ=255)
                                                                                                               aCTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCAc                            >  1:59619/1‑139 (MQ=255)
                                                                                                                   aTTCTGTGATGAAATCACTCTCGAAGAAATACCATACCATGTCAATAAGGTGAATACCAAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTATACTTCTTACGTCTAAGTGCACTAt                         <  8:154928/138‑1 (MQ=255)
                                                                                                                   aTTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTAt                         <  6:239015/138‑1 (MQ=255)
                                                                                                                     tCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGtcttct                                                                            <  1:73971/85‑1 (MQ=255)
                                                                                                                     tCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGtcttct                                                                            >  2:73971/1‑85 (MQ=255)
                                                                                                                     tCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATATCTCCAAGTAAACTTCTTACGTATAAGTGCACTAtttt                      >  5:199680/1‑139 (MQ=255)
                                                                                                                       tgtgATGAAATCACTCTCGAACAAATACCATACCATTTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCAAATCTCCAAGTAAAAATCTTAAGTCTAACTGCACTAttttta                    >  7:160296/1‑138 (MQ=255)
                                                                                                                          gatgaAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAActtc                                            >  3:288260/1‑112 (MQ=255)
                                                                                                                          gatgaAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAActtc                                            <  4:288260/112‑1 (MQ=255)
                                                                                                                             gaAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAg              >  1:167465/1‑139 (MQ=255)
                                                                                                                                 tCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAGTGTa          <  2:19438/139‑1 (MQ=255)
                                                                                                                                 tCACTCTCGAACAAATACCATACCATGTCAATAAGGGGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAGTGTa          <  1:110262/139‑1 (MQ=255)
                                                                                                                                     tctcGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACATTTATCGTCATCTCTCCAAATAAACTTCTTACGTCTAAGTGCACTATTTTTGTTTAAGTGTACttt      >  5:204432/1‑139 (MQ=255)
                                                                                                                                         gAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAGTGTACTTTTa    <  4:148993/137‑1 (MQ=255)
                                                                                                                                         gAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAGTGTACTTTTATg  <  3:179118/139‑1 (MQ=255)
                                                                                                                                         |                                                                                                                                          
ATGATATTTACAAATACCTTGTTCTTCAGCTGGCCATAAATTTCTGCATAGTCATCTACAAGTACCTGTTTATCCTCTTTTGTTTTTACATTTGTGTTCATCTTTGCCCTTACTGATTCTGTGATGAAATCACTCTCGAACAAATACCATACCATGTCAATAAGGTGAATACCCAGATCCCCCAATGATCCACTCGTCTTCTTACTATTAGCGTCATCTCTCCAAGTAAACTTCTTACGTCTAAGTGCACTATTTTTCTTGAAGTGTACTTTTATG  >  BS168_reference_genome/1126749‑1127024

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: