Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1260AC55.1% 8.4 / 87.6 69A26A (GCA→GCCstfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base C (25/13);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.43e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGG  >  NC_000913/1428885‑1429376
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            
tccggcagcgtccgtcgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAAc                                                                                                                                                                                                                                                     >  2:110908/21‑251 (MQ=255)
 ccgaccgaATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGCCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACAAAGTGTTATGACAGCCCGACGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCCGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCCAAACGTAAc                                                                                                                                                                                                                                                     >  1:564329/1‑250 (MQ=255)
         aTCCGGGGCCAGTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCGCAATCCAGCTGAAAGCAAAACGTAACAGCACCACg                                                                                                                                                                                                                                            <  2:465789/251‑1 (MQ=255)
             gtcgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATTCAGCTGAAAGCCAGACGTCACAGCACCACggagg                                                                                                                                                                                                                                        >  2:447235/8‑248 (MQ=255)
               cgttgttGATCTTGATTATGACGGCTCCGGCAGTGGCAACTGGCGTACTGTAGGATTACTCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCACGACGGACCCGGAGGGCTTTTTTGTGGGGGGAATATGGAAGTAAAGATTTCAGGCGACCTGACAGACGGCACAGGAAAACCGGCACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACCGCAACACggtggtg                                                                                                                                                                                                                                    >  2:381028/6‑251 (MQ=255)
               cgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACggtggtg                                                                                                                                                                                                                                    >  1:187014/6‑251 (MQ=255)
                               tATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTggg                                                                                                                                                                                                                      <  1:334552/249‑2 (MQ=255)
                                aTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTc                                                                                                                                                                                                                   <  2:234528/251‑1 (MQ=255)
                                   aaGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCTCCACGGTGGTGGGGAACACGGTGGGCTCAGa                                                                                                                                                                                                                <  1:160798/251‑1 (MQ=255)
                                        tccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTgggctacgagaatc                                                                                                                                                                                                           >  2:177645/2‑239 (MQ=255)
                                          cGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCg                                                                                                                                                                                                         >  1:317236/1‑251 (MQ=255)
                                           ggCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCgg                                                                                                                                                                                                        >  1:346426/1‑251 (MQ=255)
                                              agttgCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAGAAATTCTGATTAGCCAGGTAACCCAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATAGGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATg                                                                                                                                                                                                     <  2:185104/247‑1 (MQ=255)
                                               gTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCAACACGGTGGATGGGAACACGGTGGGCTCAGAGAATCCGGATgc                                                                                                                                                                                                    >  1:495757/1‑250 (MQ=255)
                                                          cggaCTGACGGATGCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAATCCGGGCGtt                                                                                                                                                                                         <  2:460696/248‑1 (MQ=255)
                                                           gtaCTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGtt                                                                                                                                                                                         >  1:175382/3‑250 (MQ=255)
                                                                gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTGCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTTAAAGCAAAACGTAACAGCACCACTGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATTAATCCGGTCGTTACAGc                                                                                                                                                                                    <  2:378134/250‑1 (MQ=255)
                                                                   ggattCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATgg                                                                                                                                                                                >  1:64105/4‑251 (MQ=255)
                                                                           cGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGAGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGa                                                                                                                                                                         <  1:139373/250‑1 (MQ=255)
                                                                                    cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAg                                                                                                                                                               <  2:171667/251‑1 (MQ=255)
                                                                                        tataAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAc                                                                                                                                                            <  1:110908/250‑1 (MQ=255)
                                                                                              ttCTGATTAGCCAGGGAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGCCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGCACGGCCGGTACAgtgtcc                                                                                                                                                     >  2:315031/1‑246 (MQ=255)
                                                                                                  gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTACCAGCACCACGGAGGTGGTGAACACGGTGGGCTAAGAGAATCCGGATGAAGACGGGCGTTACAGCATGGATGGGGAGGACGGTCAGTACAgtggcaccct                                                                                                                                                 >  2:349461/1‑242 (MQ=255)
                                                                                                    ttAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGt                                                                                                                                               >  2:125431/1‑251 (MQ=255)
                                                                                                        ccgggTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCGTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAGCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTTAGTACGGTCAGTACAGCGCTATTCTGTtggt                                                                                                                                           <  2:481989/248‑1 (MQ=255)
                                                                                                         cAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggtt                                                                                                                                          >  2:163313/1‑235 (MQ=255)
                                                                                                          aGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttg                                                                                                                                         <  1:163313/251‑18 (MQ=255)
                                                                                                              aaCACAGTGTGATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGGACAgtgtcatcctgcaggttgacg                                                                                                                                      <  2:62126/250‑21 (MQ=255)
                                                                                                                 acaGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggttt                                                                                                                                  >  1:448740/1‑227 (MQ=255)
                                                                                                                                   gccgGTTCAGGCGGGTTTTTTTGGGGGGGGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCAAAGGAAAACCGGGACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGGGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGATTACGGTCATTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCgg                                                                                                                <  2:447809/251‑1 (MQ=255)
                                                                                                                                       gtccAGGCGGGCTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCGGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAcc                                                                                                            <  1:537475/248‑1 (MQ=255)
                                                                                                                                         gcgggtgggCTTTTGTGTGGGGGGAATAGGTCAGGAAAGATTTCAGTTGTACTTAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCGGCTGAAAGAAAAACGAAACAGCCCCAGGTGTGGGGTGAACGCGCTGGCTTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAGGGTCATTTCAGCGGGATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAt                                                                                                          <  1:34094/245‑1 (MQ=255)
                                                                                                                                            ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                                                                                       >  1:320049/1‑251 (MQ=255)
                                                                                                                                               gggCTTTTTTGTGGGGTGAAGATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGAAAAACGTAACAACACCACGGTGGTGGTTAACACGCTGGCCTCGGAAAATCCGGATGAACCCGGGCGTGACAGCATGGACGTTGAGTAGGGTAAGTACAGCGTTATTCTGTTGGTGGAAGTATTCCCGCCGTCACATGCCGGGACCATTACCgt                                                                                                    <  1:369497/251‑1 (MQ=255)
                                                                                                                                                ggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGCCCATTACCgtg                                                                                                   >  2:744/1‑251 (MQ=255)
                                                                                                                                                 gCTTTTTTGTGGGGTGAATATGGCAGTAAAGCTTTCAGGTGTACTGAAAGACGGCCCAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGAGGCGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTACCATGCCGGGACCATTCCCgtgt                                                                                                  >  2:284960/1‑251 (MQ=255)
                                                                                                                                                      tttGTGGGGTGAATATGGCAGTAAATATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCTCCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGTGCTCATAGAATCCGGATGAAGCCGGGCGGTACAGCATGGATGTGTTGTACGGTCAGTACAGTGTCATCCTGCGGGGTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGa                                                                                              <  1:43405/250‑1 (MQ=255)
                                                                                                                                                       ttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGaa                                                                                             <  1:366256/250‑1 (MQ=255)
                                                                                                                                                          tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTCCTGAAAGCCGGAACAGGAACACAGGCACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCAAGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACGGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCAAATGCCGGGACCATTACCGTGTATGAAGAtt                                                                                         >  2:561195/1‑251 (MQ=255)
                                                                                                                                                          gggggTGAATAGGGCAGTAAAGATTTCAGGTGGACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGATGAAAGCAAAACGTAACAGCACCGCGGTGGTGGTGAACACGCTGGCCTCAGAAAAACCGGTTGAAGCCGGGCGTTACAGCATGGAGGTGGAGTACGTTCAGTACAGCGTTATTTTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGGATGAAGAtg                                                                                         <  2:182601/250‑2 (MQ=255)
                                                                                                                                                            gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACCGGAAAACCGGTAAAGAACTGCACCATTCAGCTGCAAGCCAGACGTAACAGCCCCACGGTGATGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGCGTACGGCAAGTACAGTGTCATCCTGCCGGCTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTtc                                                                                        >  2:493375/1‑250 (MQ=255)
                                                                                                                                                            gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAACATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGCTGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCAACGTGTATGAAGATtcc                                                                                       >  2:502073/1‑250 (MQ=255)
                                                                                                                                                              gTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAACACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGCAGATtctc                                                                                      >  1:26086/1‑250 (MQ=255)
                                                                                                                                                               tGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGACAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcaca                                                                                     >  1:502116/1‑247 (MQ=255)
                                                                                                                                                                   aatGGCTGTAAAGATTGCAGGTTAACTGAAAGCCGGCACAGGAAAACCGGTACAGAACTGCAATACCCAGCGGAAAGCAAAACTGAAAAGAACCACGGTGGGGGTGAACACGCGGGCCTCAGAAAATCCGGAGGAAGCCGGGCGTTACAGCATTGACGTGGAGTCGGGTCAGTCCAGCGTTATTCTTTTGGTGAAATTATTCCCGCCGGCACATGCCGGGACCATTACCGTGTATGAAGATGCTCAACCCg                                                                                <  2:510030/250‑1 (MQ=255)
                                                                                                                                                                        cAGTAAAGATTTCAGGTGTCCTGAAAGACGGCACAGGAAATCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTGGTGGTGACCCCGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg                                                                           <  1:137197/251‑1 (MQ=255)
                                                                                                                                                                        cAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggggacg                                                                           <  2:198451/251‑8 (MQ=255)
                                                                                                                                                                           tAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTg                                                                        >  2:156236/1‑251 (MQ=255)
                                                                                                                                                                            aaaGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCAGCGCACGCCGGGACCATCACCGGGTATGAAGATTCACAAccgtggacgctga                                                                       <  2:88738/251‑12 (MQ=255)
                                                                                                                                                                               gATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATg                                                                    >  1:180217/1‑251 (MQ=255)
                                                                                                                                                                                   tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                                                >  1:517741/1‑251 (MQ=255)
                                                                                                                                                                                    caggagtccTGAAAGACGGCACAGACAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCCGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAtttt                                                               >  2:258665/9‑251 (MQ=255)
                                                                                                                                                                                           ccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCtgtg                                                        >  1:138193/2‑251 (MQ=255)
                                                                                                                                                                                              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGc                                                       <  1:175182/249‑1 (MQ=255)
                                                                                                                                                                                              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCa                                                     <  2:446098/251‑1 (MQ=255)
                                                                                                                                                                                              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGTTGAACACGCTGGCCTCATAAAATCCGGAGGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGTTCAGTACGGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGcc                                                      <  2:341157/250‑1 (MQ=255)
                                                                                                                                                                                                 aGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATg                                                   >  2:161128/1‑250 (MQ=255)
                                                                                                                                                                                                    cGGCACAGGCAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGCTTTTCTCGGTGCCATGACgg                                               >  1:454630/1‑251 (MQ=255)
                                                                                                                                                                                                    cGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGCCgg                                               >  1:220903/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCAAAGGAAAACCGGTACAGAACTGCACAATCCAGATGAAAGAAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTAAAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACAGGCAGGGACCATTAACGTGTATGAAGATTCTCAACCCGGTAAGCTGAATGATTTTAACGGTGCCATGACgg                                               >  1:279228/1‑250 (MQ=255)
                                                                                                                                                                                                       cacaGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACggagg                                            >  1:19581/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGGTAAAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTGGTGGTGAACACTTGGGCCTCAGAAAAGCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCGGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGt                                 <  2:320049/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggc                                >  2:166983/1‑249 (MQ=255)
                                                                                                                                                                                                                            cggccaaaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGTGCTTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCAc                       <  1:200182/245‑1 (MQ=255)
                                                                                                                                                                                                                            cTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgc                       >  1:56678/1‑248 (MQ=255)
                                                                                                                                                                                                                               cacaATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGCTACGCTGAATGATTTTATCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgc                    >  2:4505/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGt               >  1:191609/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtt              >  1:188442/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaGCAAAACGTAACAGCACCACGGTGGGGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCTCTGCGCCGTTTTGCACTg      >  2:456390/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaaaCGTAACAGCACCACGGTGGTGGTGCACACGCTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGAtgg  >  2:470524/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                            
CCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGG  >  NC_000913/1428885‑1429376

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