Predicted mutation
evidence seq id position mutation freq annotation gene description
RA BS168_reference_genome 2,669,936 C→A 5.1% A145A (GCG→GCT yqbT ← putative phage baseplate assembly protein; skin element

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*BS168_reference_genome2,669,9360CA 5.1% 257.1 / 3.9 78A145A (GCG→GCTyqbTputative phage baseplate assembly protein; skin element
Reads supporting (aligned to +/- strand):  ref base C (37/37);  new base A (2/2);  total (39/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 = 1.50e-01

TCTACTTCTTCAGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCTGTGTATTGAAAATAAAGATTATCAACATAA  >  BS168_reference_genome/2669799‑2670071
                                                                                                                                         |                                                                                                                                       
tCTACTTCTTCAGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCg                                                                                                                                        >  8:169737/1‑139 (MQ=255)
tCTACTTCTTAAGCCCATTGTTTATAGTGTGCATTATTGGCACTGACAGCTTCTCGACGCACGCGGGGAAAATATCTAGCTCTTAAACTGTCAACAGCCTCATCTTGAATTCCAGTCAACAAAATCGCTTTCAGAATCg                                                                                                                                        >  4:26418/1‑139 (MQ=255)
      tcttcAGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCGTCTCGACGCACGCGGGTTAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCttttt                                                                                                                                  >  2:262298/1‑139 (MQ=255)
      tcttcAGCCCATTGGTTATAGTGTGCGTTATTGGCACAGACCGCTTCTCGACGCCCGCGGGTAAAATATCAAGCTCTTAAACTGTCATCATACTCTTCTTCACGTCCAGGCATCAGAAACTCTTTCTCAATCGCttttt                                                                                                                                  >  8:306278/1‑139 (MQ=255)
           aGCCCATTGTTTATAGTGTGCTTTTTTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  1:150543/138‑1 (MQ=255)
           aGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  8:200350/138‑1 (MQ=255)
           aGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  2:1508/138‑1 (MQ=255)
           aGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  6:303319/138‑1 (MQ=255)
            gCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  5:126136/139‑1 (MQ=255)
            gCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  7:244190/139‑1 (MQ=255)
             cccATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTTGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGccc                                                                                                                           <  6:116378/139‑1 (MQ=255)
               cATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAg                                                                                                                         >  5:316024/1‑139 (MQ=255)
                 ttGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGccc                                                                                                                           >  8:342062/1‑135 (MQ=255)
                 ttGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGccc                                                                                                                           <  7:342062/135‑1 (MQ=255)
                 ttGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAgg                                                                                                                        >  6:90441/1‑138 (MQ=255)
                  tGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  3:11563/133‑1 (MQ=255)
                  tGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            >  4:11563/1‑133 (MQ=255)
                       ataGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  3:356794/126‑1 (MQ=255)
                       ataGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              >  4:356794/1‑126 (MQ=255)
                          gtgtGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCt                                                                                                              <  7:169737/139‑1 (MQ=255)
                            gtgCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  8:79024/123‑1 (MQ=255)
                            gtgCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            >  7:79024/1‑123 (MQ=255)
                             tgCTTTATTGGCCCTCACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  7:308504/1‑139 (MQ=255)
                             tgCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCt                                                                                                              <  6:268667/136‑1 (MQ=255)
                             tgCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCt                                                                                                              >  5:268667/1‑136 (MQ=255)
                             tgCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  4:5870/1‑139 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  3:345361/119‑1 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              >  4:345361/1‑119 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            >  7:402352/1‑121 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  8:402352/121‑1 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGtata                                                                                                                    >  2:409642/1‑129 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGtata                                                                                                                    <  1:409642/129‑1 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  1:246517/1‑138 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  4:305215/1‑138 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  6:391443/1‑138 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  5:70743/1‑138 (MQ=255)
                              gCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAAt                                                                                                           >  5:380397/1‑138 (MQ=255)
                               cTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACAGTCATCAACCTCTTCTTCACGTCCAGCAATCAGAATCTCTTTCACAATCGCTTTTTGAATCCAAGGAATAGAGTATAATg                                                                                                          >  8:92945/1‑138 (MQ=255)
                                  tATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  8:204006/117‑1 (MQ=255)
                                  tATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            >  7:204006/1‑117 (MQ=255)
                                  tATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGAc                                                                                                        <  1:217672/137‑1 (MQ=255)
                                   aTTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACaa                                                                                                      <  6:316024/138‑1 (MQ=255)
                                    ttGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAAt                                                                                                     <  3:401770/138‑1 (MQ=255)
                                       gCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGccc                                                                                                                           <  5:98862/113‑1 (MQ=255)
                                       gCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGccc                                                                                                                           >  6:98862/1‑113 (MQ=255)
                                       gCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            <  5:378600/112‑1 (MQ=255)
                                       gCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            >  6:378600/1‑112 (MQ=255)
                                         aCTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTc                                                                                               >  3:259130/1‑139 (MQ=255)
                                                     tcGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATAGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATg                                                                                   >  4:181233/1‑139 (MQ=255)
                                                     tcGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATAGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATg                                                                                   >  4:187491/1‑139 (MQ=255)
                                                       gACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCGTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGtt                                                                                 <  2:193994/139‑1 (MQ=255)
                                                           cacgcGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGa                                                                                       <  1:118752/129‑1 (MQ=255)
                                                           cacgcGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGa                                                                                       >  2:118752/1‑129 (MQ=255)
                                                                 ggTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTGCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCt                                                                       >  5:334406/1‑139 (MQ=255)
                                                                          tCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTcccc                                                               <  7:106886/138‑1 (MQ=255)
                                                                          tCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTGCCCTGCTTcccc                                                               <  7:391691/138‑1 (MQ=255)
                                                                            tAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              <  7:216180/73‑1 (MQ=255)
                                                                            tAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAg                                                                                                                              >  8:216180/1‑73 (MQ=255)
                                                                              gCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCg                                                             <  1:262298/136‑1 (MQ=255)
                                                                                  ttAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATg                                                                                                          <  5:354237/87‑1 (MQ=255)
                                                                                  ttAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATg                                                                                                          >  6:354237/1‑87 (MQ=255)
                                                                                    aaaCTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTca                                                     <  5:232983/138‑1 (MQ=255)
                                                                                          tcatcaTCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAgtgt                                                                                                                >  3:164405/1‑73 (MQ=255)
                                                                                          tcatcaTCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAgtgt                                                                                                                <  4:164405/73‑1 (MQ=255)
                                                                                          tcatcaTCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATATTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTcc                                               >  5:325818/1‑138 (MQ=255)
                                                                                               atcCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGt                                          <  6:70743/138‑1 (MQ=255)
                                                                                                tcCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTcc                                        >  6:303904/1‑139 (MQ=255)
                                                                                                    cttcttCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGcc                                                                                                                            >  7:115832/1‑51 (MQ=255)
                                                                                                    cttcttCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACAGGATATGTTTGCGTTCCATGCTTCCCCCGCTGTTTCAAACTCCAACGACCCGt                                     >  5:337688/1‑138 (MQ=255)
                                                                                                     ttcttcACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCAGACTCCAACGTCCCGt                                     <  6:380397/137‑1 (MQ=255)
                                                                                                     ttcttcACGTCCAGGAATCAGAATCTCTTTCACAATAGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGt                                     <  3:181233/137‑1 (MQ=255)
                                                                                                     ttcttcACGTCCAGGAATCAGAATCTCTTTCACAATAGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGt                                     <  3:187491/137‑1 (MQ=255)
                                                                                                       cttcACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAg                                  <  7:306278/138‑1 (MQ=255)
                                                                                                          cACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCt                                <  3:26418/137‑1 (MQ=255)
                                                                                                                cAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCTGTGTAtt                         <  2:246517/138‑1 (MQ=255)
                                                                                                                                        tCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCTGTGTAtt                         >  6:357096/1‑114 (MQ=255)
                                                                                                                                        tCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCTGTGTAtt                         <  5:357096/114‑1 (MQ=255)
                                                                                                                                        tCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCTGTGTATTGAAAATAAAGATTATCAACATaa  >  8:128847/1‑137 (MQ=255)
                                                                                                                                         |                                                                                                                                       
TCTACTTCTTCAGCCCATTGTTTATAGTGTGCTTTATTGGCACTGACAGCTTCTCGACGCACGCGGGTAAAATATCTAGCTCTTAAACTGTCATCATCCTCTTCTTCACGTCCAGGAATCAGAATCTCTTTCACAATCGCTTTTTGAAGCCCAGGTATAGTGTCTAATGACAATAGATTCTGACCGGATATGTTTGCGTTCCCTGCTTCCCCCGCTGTTTCACACTCCAACGTCCCGTCAGCTGTGTATTGAAAATAAAGATTATCAACATAA  >  BS168_reference_genome/2669799‑2670071

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: