Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,429,732 G→A 100% K228K (AAG→AAA)  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,7320GA98.1% 7367.6 / 39.4 2604K228K (AAG→AAA) stfRRac prophage; putative tail fiber protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (1205/1321);  total (1244/1360)
Fisher's exact test for biased strand distribution p-value = 7.82e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 2605 total aligned reads displayed.

GCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGC‑AAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCC  >  NC_000913/1429586‑1429879
                                                                                                                                                  |                                                                                                                                                    
gCCACTGAAGCGGAAAAAAGTGCCGCAGCCTCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGc                                                                                                                                                    >  7:191489/1‑149 (MQ=255)
  cACTGAAGCGGAAAAAAGTGCCGCAGCTGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGcgt                                                                                                                                                  <  2:33816/149‑3 (MQ=255)
  cACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGcgt                                                                                                                                                  >  4:192758/1‑147 (MQ=255)
  cACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGcgt                                                                                                                                                  <  8:331164/149‑3 (MQ=255)
  cACTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGcgt                                                                                                                                                  <  1:287831/149‑3 (MQ=255)
   aCTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTc                                                                                                                                                 <  1:80474/149‑1 (MQ=255)
    cTGAAGCGGAAAAAAGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCa                                                                                                                                                <  4:262933/149‑1 (MQ=255)
            gaaaaaaGTGCCGCAGCCTCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGc                                                                                                                                        <  5:269230/147‑1 (MQ=255)
            gaaaaaaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGACACGAATGCTGCAGCGTCCCAACCATCAGCCGCCACGTCTGCATCCACCGCGGCCACGAACGCGTCAGAGGCCGc                                                                                                                                        >  1:124438/3‑149 (MQ=255)
              aaaaaGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGcca                                                                                                                                      <  3:274984/149‑3 (MQ=255)
                  aGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCActt                                                                                                                                   >  5:369896/1‑146 (MQ=255)
                  aGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTc                                                                                                                                  >  5:90178/1‑149 (MQ=255)
                  aGTGCCGCAGCCGCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTc                                                                                                                                  <  8:130452/149‑1 (MQ=255)
                          agccgCAGAGGCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACgag                                                                                                                          <  1:376112/145‑3 (MQ=255)
                           gccgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACgaga                                                                                                                         <  7:306765/146‑1 (MQ=255)
                             cgCAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATg                                                                                                                       >  5:343248/2‑149 (MQ=255)
                               cAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCCGCCGCCACGTCTGCCTCCCCCGCGGCCACGACAGCGTCAGAGGCCGCCACTTCAGCACGAGATgcg                                                                                                                     >  3:352397/1‑149 (MQ=255)
                               cAGAGTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATgcg                                                                                                                     >  3:79893/1‑149 (MQ=255)
                                agagTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATgcgg                                                                                                                    <  8:185959/149‑1 (MQ=255)
                                agagTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATgcgg                                                                                                                    >  7:159906/1‑149 (MQ=255)
                                agagTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGGCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATgcgg                                                                                                                    <  5:212264/149‑1 (MQ=255)
                                 gagTCCTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATgcggt                                                                                                                   <  2:9356/149‑2 (MQ=255)
                                      cTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGcc                                                                                                               <  3:331009/148‑1 (MQ=255)
                                      cTCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATAAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCt                                                                                                              >  7:364601/1‑149 (MQ=255)
                                       tCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCt                                                                                                              >  4:102575/1‑148 (MQ=255)
                                       tCAAAAAACGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCt                                                                                                              >  2:126872/1‑148 (MQ=255)
                                         aaaaaaCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCa                                                                                                            >  3:142985/1‑148 (MQ=255)
                                           aaaaCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCCGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCaa                                                                                                           >  3:247690/1‑147 (MQ=255)
                                           aaaaCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCaaaa                                                                                                         <  2:235969/149‑1 (MQ=255)
                                           aaaaCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCaaaa                                                                                                         <  4:351488/149‑1 (MQ=255)
                                           aaaGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCaaaa                                                                                                         <  6:80250/149‑1 (MQ=255)
                                              accccgcggcCAACAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAaagag                                                                                                      <  7:94756/144‑2 (MQ=255)
                                               cgcgGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAaagagg                                                                                                     <  4:295133/149‑3 (MQ=255)
                                               cgcgGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAaagagg                                                                                                     <  1:257252/149‑3 (MQ=255)
                                                  ggcggcCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGcag                                                                                                  >  8:51640/1‑147 (MQ=255)
                                                   gcggcCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGc                                                                                                 <  2:66926/149‑1 (MQ=255)
                                                        caccaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAat                                                                                           <  6:57797/149‑1 (MQ=255)
                                                         accaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGc                                                                                                    <  5:162518/140‑1 (MQ=255)
                                                         accaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatc                                                                                          <  6:166791/149‑1 (MQ=255)
                                                         accaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAat                                                                                           >  7:333021/1‑148 (MQ=255)
                                                          ccaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatca                                                                                         >  8:217022/1‑149 (MQ=255)
                                                          ccaGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatca                                                                                         >  1:100362/1‑149 (MQ=255)
                                                           caGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcat                                                                                        <  7:220005/149‑1 (MQ=255)
                                                           caGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcat                                                                                        <  8:27244/149‑1 (MQ=255)
                                                           caGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCCTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGGCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcat                                                                                        >  4:323261/1‑149 (MQ=255)
                                                            aGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcatc                                                                                       <  3:118269/149‑1 (MQ=255)
                                                            aGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcatc                                                                                       <  5:49121/149‑1 (MQ=255)
                                                            aGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcat                                                                                        >  7:124345/1‑148 (MQ=255)
                                                             gTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAatcat                                                                                        >  6:320811/1‑147 (MQ=255)
                                                             gTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAAtcatca                                                                                      <  1:48360/149‑1 (MQ=255)
                                                                ccGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGa                                                                                    >  6:149092/1‑148 (MQ=255)
                                                                 cGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGaaa                                                                                  <  4:46398/149‑1 (MQ=255)
                                                                 cGGTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGaaa                                                                                  <  7:245121/149‑1 (MQ=255)
                                                                   gTGCGGCGAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAaac                                                                                 >  1:199787/1‑148 (MQ=255)
                                                                         cgAAAACGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAAtcatca                                                                                      >  3:221452/1‑137 (MQ=255)
                                                                               cGTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCAtcatca                                                                    <  8:79001/149‑1 (MQ=255)
                                                                                gTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCatcatc                                                                     >  1:173532/1‑147 (MQ=255)
                                                                                gTCAGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCaa                                                                   <  3:41584/149‑1 (MQ=255)
                                                                                   aGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAgtg                                                                >  5:323303/1‑147 (MQ=255)
                                                                                   aGAAACGAATGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAgtg                                                                >  4:167988/1‑147 (MQ=255)
                                                                                     aaaCGAATGCTGCAGCGTCAAAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGcc                                                              <  6:120457/149‑1 (MQ=255)
                                                                                          aaTGCTGCAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCggtcg                                                         >  4:223664/1‑147 (MQ=255)
                                                                                            tgctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTg                                                       <  3:299202/149‑1 (MQ=255)
                                                                                              ctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGc                                                      <  2:27977/145‑1 (MQ=255)
                                                                                              ctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGc                                                      <  2:310095/145‑1 (MQ=255)
                                                                                              ctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGc                                                      <  5:243462/145‑1 (MQ=255)
                                                                                              ctgcAGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCa                                                     <  1:233191/146‑1 (MQ=255)
                                                                                                 cAGCGTCACAACAATCAGCCGCCATGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGc                                                   <  4:75747/148‑1 (MQ=255)
                                                                                                  aGCGTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCtt                                                 >  2:391252/1‑149 (MQ=255)
                                                                                                     gTCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCt                                              >  5:70405/1‑149 (MQ=255)
                                                                                                      tCACAACAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGCAAAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCt                                              >  1:348354/1‑149 (MQ=255)
                                                                                                        acaaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg                                           >  3:188828/4‑149 (MQ=255)
                                                                                                        acaaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCgg                                           <  3:167988/146‑1 (MQ=255)
                                                                                                         caaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGccg                                                             <  6:265636/128‑2 (MQ=255)
                                                                                                         caaCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCg                                            <  6:346942/145‑1 (MQ=255)
                                                                                                           aCAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAAC‑AAAATCATCAGAaacgaa                                                                              <  7:395977/111‑1 (MQ=255)
                                                                                                            cAATCAGCCGCCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGcatcat                                                                      >  8:375948/1‑118 (MQ=255)
                                                                                                                     gcCACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAg                               >  6:298347/1‑148 (MQ=255)
                                                                                                                      ccACGTCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAga                              >  7:70919/1‑147 (MQ=255)
                                                                                                                        acgtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGaaaa                           <  2:365908/146‑1 (MQ=255)
                                                                                                                        acgtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCAGCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGaaaa                           <  5:136461/146‑1 (MQ=255)
                                                                                                                          gtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAg                                                    >  8:16158/2‑122 (MQ=255)
                                                                                                                          gtCTGCCTCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAAtt                         <  6:19049/148‑1 (MQ=255)
                                                                                                                               cctCCACCGCGGCCACGAAAGCGTCGGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgc                     >  6:382740/3‑148 (MQ=255)
                                                                                                                               cctCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTct                       <  8:300785/144‑2 (MQ=255)
                                                                                                                                 tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                  >  2:211174/1‑146 (MQ=255)
                                                                                                                                 tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                  >  2:188986/1‑146 (MQ=255)
                                                                                                                                 tCCACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccag                  >  6:99113/1‑146 (MQ=255)
                                                                                                                                   cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCatcatc                                                                     >  5:184029/1‑96 (MQ=255)
                                                                                                                                   cACCGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCtgccaggg                <  8:242611/149‑6 (MQ=255)
                                                                                                                                     ccGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcg              >  4:134152/1‑149 (MQ=255)
                                                                                                                                     ccGCGGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcg              <  1:147716/149‑1 (MQ=255)
                                                                                                                                      cgcgGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcgg             >  2:15599/1‑149 (MQ=255)
                                                                                                                                       gcgGCCACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGggcggc            <  6:225008/149‑1 (MQ=255)
                                                                                                                                         ggccACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCaa          <  8:354449/147‑1 (MQ=255)
                                                                                                                                           ccACGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCaaaa        <  1:319391/149‑1 (MQ=255)
                                                                                                                                             aCGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCTAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAAc      <  5:105849/149‑1 (MQ=255)
                                                                                                                                             aCGAAAGCGTCAGAGGCCGCCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAAc      >  5:300111/1‑149 (MQ=255)
                                                                                                                                               gAAAGCGTCAGAGGCCACCACTTCAGCACGAGATGCGGTGGCCTCAAAAGAGGCAGC‑AAAATCATCAGAAACGAACGCATCATCAAGTGCCGGTCGTGCAGCTTCCTCGGCAACGGCGGCAGAAAATTCTGCCAGGGCGGCAAAAACGt    >  7:277067/1‑149 (MQ=255)
                                                                                                                                                 aGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGC‑AAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTcc  <  2:10261/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
GCCACTGAAGCATCAAAAAGTGCTGCCGCTGCAGAGTCCTCAAAAAGCGCGGCGGCCACCAGTGCCGGTGCGGCGAAAACGTCAGAAACGAATGCTTCAGCGTCACTACAATCAGCAGCCACATCTGCATCCACCGCGACCACGAAGGCATCAGAAGCTGCGACCTCGGCCCGGGATGCGGCGGCCTCAAAAGAAGCGGC‑AAAATCATCAGAAACGAACGCATCATCAAGCGCCAGTAGTGCAGCTTCCTCGGCAACGGCGGCAGGAAATTCCGCGAAGGCGGCAAAAACGTCC  >  NC_000913/1429586‑1429879

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