Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 1,204,963 T→A 22.8% intergenic (‑47/‑202) yjzD ← / → yjaU forespore targeted protein/conserved hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome1,204,9630TA22.8% 113.0 / 9.1 67intergenic (‑47/‑202)yjzD/yjaUforespore targeted protein/conserved hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/18);  new base A (2/14);  total (35/32)
Fisher's exact test for biased strand distribution p-value = 3.64e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGCTGTTCATGCTTGCAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCTTA  >  BS168_reference_genome/1204837‑1205101
                                                                                                                                 |                                                                                                                                           
aCGCGGTTGATGATTGCAACAATTAATCATGCAAAATGTGATAAACTTAACGTCAAAATAAAGGCAATAATATATCCCACGTTTCCACTTCGTCTTTCTTATTTTATTATAATAGTTATCATACCAAAAAAAAAGCGa                                                                                                                                 <  4:154323/138‑1 (MQ=255)
 cGCTGTTCATGCTTGCAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTATATTCTAAAAGTTATCATACCAAACAAAACTCGAAg                                                                                                                               >  6:162000/1‑139 (MQ=255)
      ttCATGCTTGCAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCTTTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAggggg                                                                                                                           >  5:109846/1‑138 (MQ=255)
      ttCATGCTTGCAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAggggg                                                                                                                           >  5:143780/1‑138 (MQ=255)
      ttCATGCTTGCAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAggggg                                                                                                                           >  3:145942/1‑138 (MQ=255)
      ttCAGGCTTCCAACCAAGTAGCTTGCAATATGTGAAAAAAGGAACGTCAAAATAACGGAAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTAAAAAAGTTAACAAACCATACAAAAAGCGAAgggg                                                                                                                            <  6:109846/137‑1 (MQ=255)
               caaaCAAAGTACCTTGCAATATGTTTAAAATTGACCGTCAAATAAATGGAAAAATAATACCCCCAGTTGCCACCTCGTCCTCCTCATTATTTTATAAAATTAATATACCCATACAAAACGCGAAGGGGGAATCCAAtt                                                                                                                  <  6:281922/136‑1 (MQ=255)
                cAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACATAGTCCTACTCAATCTATTCTAAAAGTTATCATACAATACAAAAAGCGAAGGCGGAATCCAAttat                                                                                                                >  5:281922/1‑137 (MQ=255)
                   ccAAGTACCATGCCATATGGTATAACAGTAACGACCAAATAAAGGCAAAAATAAATCGCACGTTTCCACCTCGTCCTTCTCATTATATTAAAAAAGTAATAATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGc                                                                                                             <  2:97924/139‑1 (MQ=255)
                             tGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCTGATAGTGAg                                                                                                   >  6:180250/1‑139 (MQ=255)
                                  tatGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGa                                                                                               >  6:272971/1‑138 (MQ=255)
                                  tatGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGa                                                                                               >  5:100256/1‑138 (MQ=255)
                                     ttgATAAAAGGAACGTCCAAATAAAGGCAATAATAAATCTCACGTTTCCACATCTTCCTTCTCATTCTATTCTAAAAGTTAAAATAACAAACAAAAACATAAGGGGGAATCAAATTTTAGAAGAAAGTGAGGTGACaa                                                                                            >  1:284273/2‑138 (MQ=255)
                                      tgATAAAAGCAACGTCAAAAAAAAGGAAATAATATACCGCACGTTTCCACCTCGTCCTTCACATTCTATTCTAAAAGTTATCAAACCAAACAAAACGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACaaaa                                                                                          <  3:105504/139‑1 (MQ=255)
                                        aTAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACaaaa                                                                                          >  5:82165/1‑137 (MQ=255)
                                        aTAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGAAGAAAGTGAGGTGACaaaa                                                                                          >  6:65314/1‑137 (MQ=255)
                                        aTAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACaaaa                                                                                          >  3:38945/1‑137 (MQ=255)
                                         tAAAAGTAATGTCAAAAAAACGCAAAAAATATCTCCCACGTTTCAACCTGGTCTTTATCTTTCTATTATAAAATTTAACATACAAAAAAAAACGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAaga                                                                                        <  4:38945/138‑1 (MQ=255)
                                         tAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTACCACCTCGGCCTTCTCATTATATTCTAAAAGTTATCAAAACATACAAAAAGCGAAGGGGGAATCAAAATTATGCAGAAAGTGAGGTGACaaaa                                                                                          >  6:174318/1‑136 (MQ=255)
                                               caccGTCCAAATAAAGCCAATAATATATCGCACGTTCCCACCTCGGCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACTAAAGAGTCTGAt                                                                                 <  7:78086/136‑1 (MQ=255)
                                                     ccAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTca                                                                           <  6:143780/139‑1 (MQ=255)
                                                     caaaaTAAAGCAAATAATATATCCCACGTTTCCACCTCGTCCTTCTAATTCTATTATAAAAGTTAACAAACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTca                                                                           <  3:114261/137‑1 (MQ=255)
                                                      cAAATAAAGCCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATAAAAAACGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTcac                                                                          <  8:83598/139‑1 (MQ=255)
                                                              ggCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGTTGACAAAAGAGTCTGATTCGTCACACTTTAAt                                                                  >  5:292027/1‑139 (MQ=255)
                                                              ggCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAAt                                                                  >  7:76306/1‑139 (MQ=255)
                                                                 aataatATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTTTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGt                                                               <  8:211307/139‑1 (MQ=255)
                                                                 aataatATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACCCTTTAATgg                                                                >  4:34475/1‑138 (MQ=255)
                                                                 aataatATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATgg                                                                >  2:114167/1‑138 (MQ=255)
                                                                  ataataTAGCGCACGTTTCCCCCTCTTCATTCTCATTCTATTCTAAATTTTATAAAAACAAACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTa                                                              <  6:100256/139‑1 (MQ=255)
                                                                  aaaataTATCGCCCTATCCCCCCTCGTCCTTCTCATTAAATTCTAAAAGTTATCAAACAAAACAAAACGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTa                                                              <  6:78888/137‑1 (MQ=255)
                                                                   ataaattatCGCACGTATCCACCTGGCCCTTCTCATACTATTCTAAAAGTTACCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTa                                                              <  8:3430/133‑1 (MQ=255)
                                                                       aatttCGCACGTGTCCACTGCGTCCTTCTCATTATATTGAAAAAGTATCAAAACCCATCAAAACGCGCACGCGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTaaaa                                                           <  2:211137/135‑1 (MQ=255)
                                                                         tatCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAata                                                         >  5:148759/1‑138 (MQ=255)
                                                                           tCGCACGTTACCACCTCGTCCTTCTAATTCTATTCAAAAAGTTATCATACAATAAAAAACGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATaa                                                      <  5:73406/139‑1 (MQ=255)
                                                                            cGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATaa                                                      >  7:55220/1‑138 (MQ=255)
                                                                            cGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATATTAGCTGAAAGTGAGGTGACAAAAGAGTCTGCTACTTCACACTTTAATGGGAATATATaa                                                      >  1:211137/1‑138 (MQ=255)
                                                                              cACGTTTCCACCTCGTCCTTCTCATTGTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGACACACTTAAAAGGTAAAATATAAAAc                                                   >  5:49742/1‑139 (MQ=255)
                                                                              cACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAATTGAGGTGACAAAAGAGTCTGATTCGTCGCACTTTAATGGTAAAATATAAAAc                                                   >  3:196455/1‑139 (MQ=255)
                                                                              cACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAATTATAAAAc                                                   >  5:14216/1‑139 (MQ=255)
                                                                              cACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAAc                                                   >  3:202450/1‑139 (MQ=255)
                                                                              cACGTTTCCACCTCGTCCTTCTCAATCTATTCTAAAAGTTATCATACCATACATAAAGCGAACCGCGAATCCAATTTTATCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAAc                                                   >  3:195778/1‑139 (MQ=255)
                                                                                cGTTTCCACCTCGTCCTTCTCATTCTATTCTGAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGt                                                 >  5:86624/1‑139 (MQ=255)
                                                                                  ctcccACCGCGCCCTTCTCATACAATTCAAAAAGTTAACAAACCAAACAAAACGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCa                                               <  2:63828/136‑1 (MQ=255)
                                                                                   ttCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCaa                                              >  3:96618/1‑139 (MQ=255)
                                                                                    tCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCaa                                              >  7:239477/1‑138 (MQ=255)
                                                                                     cccccTCGTCATTCTCATTAAATTCTAAAAGTTATCAACCCAAACAAAAAGCGCAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  8:76306/135‑1 (MQ=255)
                                                                                     ccACCTCTTCCTTTTAATTAAATTATAAAATTTATCATAACAAAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  6:14216/138‑1 (MQ=255)
                                                                                     ccACCTCTTCCTTCTCATTCTATTCTAAAATTTATCATACAAAAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  5:272971/138‑1 (MQ=255)
                                                                                     ccACCTCGTCGTTCTAATTCAATTATAAAAGTTATCATACCAAAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  6:82165/138‑1 (MQ=255)
                                                                                     ccACCTCGTCCTTCTCATTCTATTAAAAAATTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  8:55220/138‑1 (MQ=255)
                                                                                     ccACCTCGTCCTTCTAATTATATACAAAAAGTTATCAAACCAAAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  4:266662/138‑1 (MQ=255)
                                                                                     ccACCTCGTCATTATAATTATATTATAAAATTTTTCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACTAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  6:49742/138‑1 (MQ=255)
                                                                                     ccAACTCGTCCTACTAATTCTATTCAAAAAGTTATAATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  3:272986/138‑1 (MQ=255)
                                                                                     caaCCTCGTCCTTCTCATTCTATTATAAAATTTATCATACAATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  6:292027/136‑1 (MQ=255)
                                                                                     acACCTCGTCCTTCTCATTATATTATAAAATTTAACATACCATAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAAt                                             <  1:141968/137‑1 (MQ=255)
                                                                                           cGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGtat                                       >  3:118333/1‑138 (MQ=255)
                                                                                             tCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACCCTTTAATGGTAAAATATAAAACGTCAATCAGTataa                                     >  3:158355/1‑138 (MQ=255)
                                                                                             tCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTataa                                     >  5:210330/1‑138 (MQ=255)
                                                                                             tCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTataa                                     >  2:290315/1‑138 (MQ=255)
                                                                                             tCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTataa                                     >  3:219982/1‑138 (MQ=255)
                                                                                                           ttctaAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTaa                      <  4:145942/139‑1 (MQ=255)
                                                                                                            tcataaaTTTATCATACCAAAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTaa                      <  6:86624/134‑1 (MQ=255)
                                                                                                                aaaGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGc                 <  4:96618/139‑1 (MQ=255)
                                                                                                                 aattaaTCATACCAAAAAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCt                <  4:202450/134‑1 (MQ=255)
                                                                                                                 aaGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCt                <  8:239477/139‑1 (MQ=255)
                                                                                                                     tATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGa            >  7:245960/1‑139 (MQ=255)
                                                                                                                     tATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGa            >  1:182740/1‑139 (MQ=255)
                                                                                                                               catacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTAcc     >  6:156033/1‑138 (MQ=255)
                                                                                                                               catacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCt    >  8:275234/1‑139 (MQ=255)
                                                                                                                                atacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCt    >  4:180738/1‑138 (MQ=255)
                                                                                                                                atacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCt    >  2:244601/1‑138 (MQ=255)
                                                                                                                                 tacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTa               >  4:38892/1‑126 (MQ=255)
                                                                                                                                 tacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCTTa  >  7:219217/1‑139 (MQ=255)
                                                                                                                                 tacAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCTTa  >  8:200886/1‑139 (MQ=255)
                                                                                                                                 |                                                                                                                                           
ACGCTGTTCATGCTTGCAACCAAGTAGCATGCCATATGTGATAAAAGGAACGTCCAAATAAAGGCAATAATATATCGCACGTTTCCACCTCGTCCTTCTCATTCTATTCTAAAAGTTATCATACCATACAAAAAGCGAAGGGGGAATCCAATTTTAGCAGAAAGTGAGGTGACAAAAGAGTCTGATTCGTCACACTTTAATGGTAAAATATAAAACGTCAATCAGTATAATAGCCAAAAATGTAAAAGGCTAGGACTTTACCTTA  >  BS168_reference_genome/1204837‑1205101

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: