Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 4,066,077 A→G 5.8% intergenic (+115/‑133) cotNE → / → yxeD inner spore coat protein/conserved protein of unknown function

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome4,066,0770AG 5.8% 212.4 / 2.9 69intergenic (+115/‑133)cotNE/yxeDinner spore coat protein/conserved protein of unknown function
Reads supporting (aligned to +/- strand):  ref base A (28/37);  new base G (2/2);  total (30/39)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.56e-01

CAGGGGAAGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGATCATTCATGCA  >  BS168_reference_genome/4065943‑4066214
                                                                                                                                      |                                                                                                                                         
cAGGGGAAGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATc                                                                                                                                        <  1:226860/138‑1 (MQ=255)
cAGGGGAAGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCt                                                                                                                                       >  4:120006/1‑139 (MQ=255)
  ggggAAGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGATGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCtt                                                                                                                                      >  2:235228/1‑138 (MQ=255)
   gggAAGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTa                                                                                                                                    >  1:117745/1‑139 (MQ=255)
      aaGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCAAGTATACAATCTTTACCg                                                                                                                                 >  6:96765/1‑139 (MQ=255)
         gCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGttt                                                                                                                              <  5:83568/139‑1 (MQ=255)
             tttCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACGGTTttctt                                                                                                                          >  7:156596/1‑137 (MQ=255)
                cTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTttcattc                                                                                                                        <  4:314945/138‑1 (MQ=255)
                   aGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAATAGGGTTTGTACCCAGTATACAATCTTTACCGTTttcattca                                                                                                                       >  7:9506/1‑136 (MQ=255)
                       aaaCAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATag                                                                                                                <  1:23554/139‑1 (MQ=255)
                             aGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGc                                                                                                          >  5:182021/1‑139 (MQ=255)
                              gAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCt                                                                                                         <  7:300039/139‑1 (MQ=255)
                                  aaaaGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATagag                                                                                                              <  7:305450/130‑1 (MQ=255)
                                  aaaaGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATagag                                                                                                              >  8:305450/1‑130 (MQ=255)
                                  aaaaGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAgg                                                                                                     <  7:134020/139‑1 (MQ=255)
                                    aaGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGaa                                                                                                   <  7:132754/139‑1 (MQ=255)
                                         gAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTttca                                                                                                                           >  4:264481/1‑110 (MQ=255)
                                         gAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTttca                                                                                                                           <  3:264481/110‑1 (MQ=255)
                                            tACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTg                                                                                           <  6:259716/139‑1 (MQ=255)
                                                    tcCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTcc                                                                                    >  5:31220/1‑138 (MQ=255)
                                                    tcCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCct                                                                                   >  7:5544/1‑139 (MQ=255)
                                                    tcCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCct                                                                                   >  5:122734/1‑139 (MQ=255)
                                                    tcCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCct                                                                                   >  8:104027/1‑139 (MQ=255)
                                                         tctcAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACGATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCttt                                                                               >  7:290148/1‑138 (MQ=255)
                                                         tctcAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACACAGTATACGAACATTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCttt                                                                               >  7:295069/1‑138 (MQ=255)
                                                               ggTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGaa                                                                        >  3:297729/1‑139 (MQ=255)
                                                                               tttCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAg                                                        <  6:31220/139‑1 (MQ=255)
                                                                                     atatGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTg                                                     <  3:300593/136‑1 (MQ=255)
                                                                                     atatGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTg                                                     <  3:129828/136‑1 (MQ=255)
                                                                                     atatGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTg                                                     <  3:120006/136‑1 (MQ=255)
                                                                                     atatGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTg                                                     <  2:28935/136‑1 (MQ=255)
                                                                                     atatGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTg                                                     <  2:117745/136‑1 (MQ=255)
                                                                                          ttAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGTGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTATTTTCATTCAAGTTGCATGataa                                             >  8:29982/1‑139 (MQ=255)
                                                                                          ttAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGataa                                             >  2:42133/1‑139 (MQ=255)
                                                                                             aaaCGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATAc                                          >  3:356667/1‑139 (MQ=255)
                                                                                             aaaCGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATAc                                          <  6:182021/139‑1 (MQ=255)
                                                                                                   aaCATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTATCATTAAAAGAAAGAGGGGCTCAGGAAAATAATTCAAGTTACTCTTTGACGTAAGAACTTTTCATTCAAGAAGCATGATAATACCGACt                                     >  8:57738/1‑138 (MQ=255)
                                                                                                           cATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCaa                             <  8:5544/138‑1 (MQ=255)
                                                                                                           cATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCaa                             <  6:122734/138‑1 (MQ=255)
                                                                                                           cATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCaa                             <  4:297729/138‑1 (MQ=255)
                                                                                                           cATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAAt                            <  7:126599/139‑1 (MQ=255)
                                                                                                            ctAAAAAGGGTGTGTACCCAGTATACAATCTTTACAGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATc                           <  7:299159/138‑1 (MQ=255)
                                                                                                            aTAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATc                           <  4:259588/139‑1 (MQ=255)
                                                                                                             tAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATc                           >  2:257243/1‑138 (MQ=255)
                                                                                                               aaaaGGTGTTGTACCCATTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGATGAGGAAAATAATTGTAGTTCCTCTTTGACGTAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGAt                        <  6:130066/139‑1 (MQ=255)
                                                                                                                  aGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTAAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCATCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCTACTCACGTCAATCGATACa                     >  6:213895/1‑139 (MQ=255)
                                                                                                                      tttGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGataata                                           <  4:127512/113‑1 (MQ=255)
                                                                                                                      tttGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGataata                                           >  3:127512/1‑113 (MQ=255)
                                                                                                                        tGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAgg               >  1:82134/1‑139 (MQ=255)
                                                                                                                            cccAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAt            <  1:310553/138‑1 (MQ=255)
                                                                                                                            cccAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAt            <  5:96765/138‑1 (MQ=255)
                                                                                                                            cccAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAt            <  7:155033/138‑1 (MQ=255)
                                                                                                                            cccAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAAAACCGACTCACGTCAATCGATACATGGAGGGAt            <  2:185880/138‑1 (MQ=255)
                                                                                                                            cccAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCGGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAt            <  6:253831/138‑1 (MQ=255)
                                                                                                                              cAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcat         >  1:83603/1‑139 (MQ=255)
                                                                                                                              cAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcat         >  1:99543/1‑139 (MQ=255)
                                                                                                                              cAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcat         >  5:286926/1‑139 (MQ=255)
                                                                                                                               aGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAGTCGATACATGGAGGGAgcat         >  7:165439/1‑138 (MQ=255)
                                                                                                                                gTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattc       <  1:235228/139‑1 (MQ=255)
                                                                                                                                 tataGATTCTTTACGGTTTTGATTAAATGATAGAGGGGCTGAGGTAAATTATTGTATTTGCTCTTTGACGGATTTTCTTTTGTTTCAAGTTGCATGATAATACCGACTCACGTCTATCGTTACATGGAGGGAgcattc       <  7:57738/138‑1 (MQ=255)
                                                                                                                                 tataCAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattca      >  3:274128/1‑139 (MQ=255)
                                                                                                                                 tataCAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattc       <  2:83603/138‑1 (MQ=255)
                                                                                                                                   taCGATCTTTACCGTTTTGATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattcat     <  8:295069/138‑1 (MQ=255)
                                                                                                                                   taCGATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattcat     <  8:290148/138‑1 (MQ=255)
                                                                                                                                   taCAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattcat     <  7:104027/138‑1 (MQ=255)
                                                                                                                                   taCAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGAtcattcat     <  6:286926/138‑1 (MQ=255)
                                                                                                                                    aCAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATAGATGGAGGGATCATTCATg    >  5:53373/1‑138 (MQ=255)
                                                                                                                                     cAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGATCATTCATGCa  <  4:356667/139‑1 (MQ=255)
                                                                                                                                     cAATCTTTACCGTGTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGATCATTCATGCa  <  6:86714/138‑1 (MQ=255)
                                                                                                                                      aaTCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAAttt                                                                     >  3:9566/1‑71 (MQ=255)
                                                                                                                                      aaTCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAAttt                                                                     <  4:9566/71‑1 (MQ=255)
                                                                                                                                      aaTCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGa                         <  2:72857/115‑1 (MQ=255)
                                                                                                                                      aaTCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGa                         >  1:72857/1‑115 (MQ=255)
                                                                                                                                      |                                                                                                                                         
CAGGGGAAGGCGGTTTCTGAGCGAAACAAAGAAAAAAAGCGGAATACGATTCTCCCATCTCAGGGTTTCGTTTTCCGCATTTCTTATATGTTAAAACGGAACATTTACATAAAAAGGGTTTGTACCCAGTATACAATCTTTACCGTTTTCATTCAAAGATAGAGGGGCTGAGGAAAATAATTGTAGTTCCTCTTTGACGGAATTTCTTTTCATTCAAGTTGCATGATAATACCGACTCACGTCAATCGATACATGGAGGGATCATTCATGCA  >  BS168_reference_genome/4065943‑4066214

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: