Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,816 C→T 100% S256S (AGC→AGT)  stfR → Rac prophage; putative tail fiber protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,8160CT98.3% 7006.1 / 30.9 2462S256S (AGC→AGT) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base C (22/22);  new base T (1301/1108);  total (1326/1136)
Fisher's exact test for biased strand distribution p-value = 6.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.03e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2495 total aligned reads displayed.

GAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGC‑AAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGT  >  NC_000913/1429670‑1429961
                                                                                                                                                   |                                                                                                                                                 
             cAGCGTCACAACAATCAGCCGCCATGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGc                                                                                                                                     <  4:75747/148‑1 (MQ=255)
                cGTCACAACAATCAGCCGCCACGTCTGCCTCCACCCCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTcc                                                                                                                                 >  8:323193/1‑149 (MQ=255)
                  tCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTc                                                                                                                               <  3:350900/149‑1 (MQ=255)
                  tCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTc                                                                                                                               <  4:145410/149‑1 (MQ=255)
                    acaaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg                                                                                                                             >  8:221927/4‑149 (MQ=255)
                    acaaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCg                                                                                                                              <  6:249697/145‑1 (MQ=255)
                    acaaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCACCACGAGATGCGGTGGCTTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg                                                                                                                             >  3:134262/4‑149 (MQ=255)
                     caaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg                                                                                                                             <  2:348354/147‑1 (MQ=255)
                      aaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGc                                                                                                                            >  4:12001/2‑148 (MQ=255)
                       aCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCaa                                                                                                                          <  1:233507/149‑1 (MQ=255)
                       aCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCaa                                                                                                                          >  3:308913/1‑149 (MQ=255)
                         aaTCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcg                                                                                                                        >  4:174103/1‑149 (MQ=255)
                           tCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggc                                                                                                                      <  3:164945/149‑1 (MQ=255)
                           tCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggc                                                                                                                      >  5:358740/1‑149 (MQ=255)
                              gccgcCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggcgg                                                                                                                    >  5:331933/4‑148 (MQ=255)
                                 gcCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAg                                                                                                                 >  8:368116/1‑148 (MQ=255)
                                  ccACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAga                                                                                                                >  7:70919/1‑147 (MQ=255)
                                    acgtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCAGCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGaaaa                                                                                                             <  5:136461/146‑1 (MQ=255)
                                         tGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCg                                                                                                                              <  6:44257/127‑1 (MQ=255)
                                         tGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAAcggcgg                                                                                                                    >  1:201672/1‑137 (MQ=255)
                                           cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCCACGGCGGCAGAAAATTCtgc                                                                                                       >  6:332753/3‑148 (MQ=255)
                                           cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTctg                                                                                                        <  8:64821/145‑3 (MQ=255)
                                           cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTct                                                                                                         <  5:162793/144‑2 (MQ=255)
                                             tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                                                    >  2:211174/1‑146 (MQ=255)
                                             tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                                                                                                    >  5:304409/1‑146 (MQ=255)
                                              ccACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagg                                                                                                   >  8:34353/1‑145 (MQ=255)
                                              ccACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccagg                                                                                                   >  6:191181/1‑145 (MQ=255)
                                               cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAACAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCCACGGCGGCAGAAAATTCtgccaggg                                                                                                  >  3:102721/1‑144 (MQ=255)
                                               cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                                                  >  1:100674/1‑144 (MQ=255)
                                               cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                                                  >  3:857/1‑144 (MQ=255)
                                                  cgcgGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcgg                                                                                               >  5:189896/1‑149 (MQ=255)
                                                   gcgGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCCCGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcggc                                                                                              <  4:330632/149‑1 (MQ=255)
                                                    cggccACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcggc                                                                                              >  1:174552/4‑148 (MQ=255)
                                                    cggccACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCa                                                                                             <  2:79448/146‑1 (MQ=255)
                                                        cACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGCAACGACCGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCaaaaa                                                                                         >  7:152496/1‑149 (MQ=255)
                                                         aCGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGACGTGCAGCTTCCTCGGCAACGGCGGCAGAACATTCTGCCAGGGCGGCCAAAAc                                                                                        >  6:87975/1‑149 (MQ=255)
                                                         aCGAAAGCGTCAGAGGCCGCAACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAAc                                                                                        <  3:196809/149‑1 (MQ=255)
                                                           gAAAGCGTCAGAGGCCGGCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCCTCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGt                                                                                      <  1:283834/149‑1 (MQ=255)
                                                           gAAAGCGTCAGAGGCCGCCACTTCAGCAGGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGt                                                                                      <  8:398907/149‑1 (MQ=255)
                                                           gAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                                                                                                  >  3:283169/1‑132 (MQ=255)
                                                            aaagcgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                                                     >  6:137908/7‑149 (MQ=255)
                                                             aGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGC‑AAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTcc                                                                                    <  2:10261/149‑1 (MQ=255)
                                                              agcgtaagagGCCGCCACTTCAGCACGAGATGCGGTGGACTAAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGGCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTc                                                                                     <  4:270391/141‑1 (MQ=255)
                                                                cgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCgag                                                                                 <  7:34353/147‑1 (MQ=255)
                                                                cgtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCgag                                                                                 <  6:307906/147‑1 (MQ=255)
                                                                 gtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCgaga                                                                                <  4:219441/148‑1 (MQ=255)
                                                                 gtCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCgaga                                                                                <  6:358206/148‑1 (MQ=255)
                                                                  tCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGAc                                                                               <  1:79380/149‑1 (MQ=255)
                                                                    agagGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGa                                                                             <  7:180988/149‑1 (MQ=255)
                                                                                   agcacgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAATACGTCCGAGACGAATGCCAGGtcatct                                                               <  2:278398/141‑1 (MQ=255)
                                                                                   agcacgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGtcatct                                                               >  1:30799/8‑148 (MQ=255)
                                                                                   agcacgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGCtcatct                                                               >  4:240372/8‑148 (MQ=255)
                                                                                    gcacgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGCACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGCAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCCTCTg                                                              >  5:234892/7‑148 (MQ=255)
                                                                                    gcacgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTg                                                              >  6:221662/7‑148 (MQ=255)
                                                                                     cacgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGa                                                             >  3:261258/6‑148 (MQ=255)
                                                                                      acgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg                                                                                                                             <  7:236529/79‑1 (MQ=255)
                                                                                      acgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCTGGTCATCTGaaa                                                           <  1:382251/145‑1 (MQ=255)
                                                                                      acgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                                                           <  1:140318/145‑1 (MQ=255)
                                                                                      acgagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAAGCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                                                           <  6:108705/145‑1 (MQ=255)
                                                                                        gagATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGaaa                                                           >  2:150916/3‑147 (MQ=255)
                                                                                          gATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAcag                                                        >  1:354995/1‑146 (MQ=255)
                                                                                              cggtggCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACagcagc                                                    <  2:302634/148‑1 (MQ=255)
                                                                                                 tggcCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCggaa                                                >  5:85977/2‑147 (MQ=255)
                                                                                                   gcCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCGTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCCGCAAAACCGTCCGAGACGAATGCCAGGTCATCTGAAACAGCCGCGGAAcg                                              >  5:256393/1‑148 (MQ=255)
                                                                                                     cTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                                                                                                       >  5:214929/1‑90 (MQ=255)
                                                                                                     cTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTcc                                                                                    >  7:37982/1‑109 (MQ=255)
                                                                                                     cTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCgga                                                 <  1:262459/144‑2 (MQ=255)
                                                                                                     cTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAAcg                                              >  8:251679/1‑146 (MQ=255)
                                                                                                      tCAAAAGAGGCAGC‑AAAATCATCAGAAACGCACGCATCATCAAGTGCCGGTCGTTCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGCATGCCAGGTCATCTGAAACAGCAGCGGAACggag                                           >  7:317211/1‑149 (MQ=255)
                                                                                                      tCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCa                                                                       <  7:112702/121‑1 (MQ=255)
                                                                                                        aaaaGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAgtgc                                                                                                                                                 >  2:122214/1‑45 (MQ=255)
                                                                                                        aaaaGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGAc                                                                               <  3:326991/111‑1 (MQ=255)
                                                                                                        aaaaGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAAc                                                          >  1:121584/1‑132 (MQ=255)
                                                                                                        aaaaGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcg                                         >  2:44328/1‑149 (MQ=255)
                                                                                                        aaaaGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAgcg                                         >  6:165709/1‑149 (MQ=255)
                                                                                                        aaaaGAGGCAGC‑AAAATAATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGATACGAATCCCAGGTCATCTGAAAAAGCAGAGGAACGGAgcg                                         >  6:280588/1‑149 (MQ=255)
                                                                                                         aaaGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGa                                                                             <  1:327888/112‑1 (MQ=255)
                                                                                                          aagaggCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcg                                                                                                >  1:52587/1‑92 (MQ=255)
                                                                                                             aggcagC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTct                                    >  7:212265/6‑148 (MQ=255)
                                                                                                                     aatatcatcAGTACCGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCgg                             <  3:48807/146‑1 (MQ=255)
                                                                                                                      aaaTCATCAGAAACGCACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGc                            >  8:135155/1‑149 (MQ=255)
                                                                                                                      aaaTCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGc                            <  4:36323/149‑1 (MQ=255)
                                                                                                                          catcaGAAACGAACGCATCATCAAGTGCCGGACGTGCAGCTTCATCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGAAAAAACGTCCGAGACGAATGCCAGGTCATCAGAAACAGCAGCGGAACGGCGCGCCTCTGCCGCGGCAgac                        >  5:42728/1‑147 (MQ=255)
                                                                                                                           atcaGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAgacg                       >  7:44970/1‑146 (MQ=255)
                                                                                                                            tcaGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCgg                             <  7:298260/142‑1 (MQ=255)
                                                                                                                              aGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCGTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAgacgca                     >  8:292040/1‑143 (MQ=255)
                                                                                                                               gAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAgg                                                                     >  5:106395/1‑99 (MQ=255)
                                                                                                                                aaaCGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCaaa                   >  2:372607/1‑148 (MQ=255)
                                                                                                                                 aacgaacgCATCATCAAGTGCCGGACGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCaaaa                  >  3:230664/1‑148 (MQ=255)
                                                                                                                                      acgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCg            >  5:6201/1‑149 (MQ=255)
                                                                                                                                      acgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCg            >  8:191234/1‑149 (MQ=255)
                                                                                                                                       cgCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCgg           <  6:69029/149‑1 (MQ=255)
                                                                                                                                          atcatcAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGAc                                                                               >  7:27396/1‑78 (MQ=255)
                                                                                                                                          atcatcAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCCGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCGGcgg        >  4:304873/1‑147 (MQ=255)
                                                                                                                                           tcatcaAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCgg                                                  <  2:24191/106‑1 (MQ=255)
                                                                                                                                           tcatcaAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAGCggcggc       >  5:280661/1‑149 (MQ=255)
                                                                                                                                               caAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTATGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcgggg   >  8:166649/1‑146 (MQ=255)
                                                                                                                                                aaGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcgggga  <  2:180110/149‑5 (MQ=255)
                                                                                                                                                aaGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcgggga  <  2:32550/149‑5 (MQ=255)
                                                                                                                                                aaGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCAGCCAGGGCGGCAAAAACGTCCGAGACGAATGCCAGGTCATCTGAAACAGCAGCGGAACGGAGCGCCTCTGCCGCGGCAGACGCAAAAACAgcggcggcggg    >  8:104410/1‑145 (MQ=255)
                                                                                                                                                   |                                                                                                                                                 
GAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGC‑AAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCCGAGACGAACGCCAGGTCTTCTGAAACGGCAGCGGGACAGAGCGCCTCGGCTGCGGCAGGCTCAAAAACAGCGGCTGCGTCGT  >  NC_000913/1429670‑1429961

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: ‑