Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1280AG54.4% 23.6 / 81.7 69K27R (AAA→AGA) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base A (15/16);  new base G (24/13);  total (40/29)
Fisher's exact test for biased strand distribution p-value = 2.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.74e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGG  >  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)
               cgttgttGATCTTGATTATGACGGCTCCGGCAGTGGCAACTGGCGTACTGTAGGATTACTCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCACGACGGACCCGGAGGGCTTTTTTGTGGGGGGAATATGGAAGTAAAGATTTCAGGCGACCTGACAGACGGCACAGGAAAACCGGCACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACCGCAACACggtg‑g‑tg                                                                                                                                                                                                                                    >  2:381028/6‑251 (MQ=255)
               cgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACggtg‑g‑tg                                                                                                                                                                                                                                    >  1:187014/6‑251 (MQ=255)
                               tATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTggg                                                                                                                                                                                                                      <  1:334552/249‑2 (MQ=255)
                                aTGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTc                                                                                                                                                                                                                   <  2:234528/251‑1 (MQ=255)
                                   aaGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCTCCACGGTG‑G‑TGGGGAACACGGTGGGCTCAGa                                                                                                                                                                                                                <  1:160798/251‑1 (MQ=255)
                                        tccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTgggctacgagaatc                                                                                                                                                                                                           >  2:177645/2‑239 (MQ=255)
                                          cGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCg                                                                                                                                                                                                         >  1:317236/1‑251 (MQ=255)
                                           ggCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCgg                                                                                                                                                                                                        >  1:346426/1‑251 (MQ=255)
                                              agttgCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAGAAATTCTGATTAGCCAGGTAACCCAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATAGGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATg                                                                                                                                                                                                     <  2:185104/247‑1 (MQ=255)
                                               gTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCAACACGGTG‑GATGG‑GAACACGGTGGGCTCAGAGAATCCGGATgc                                                                                                                                                                                                    >  1:495757/1‑250 (MQ=255)
                                                          cggaCTGACGGATGCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAATCCGGGCGtt                                                                                                                                                                                         <  2:460696/248‑1 (MQ=255)
                                                           gtaCTGACGGATTCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGtt                                                                                                                                                                                         >  1:175382/3‑250 (MQ=255)
                                                                gACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTGCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTTAAAGCAAAACGTAACAGCACCACTGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATTAATCCGGTCGTTACAGc                                                                                                                                                                                    <  2:378134/250‑1 (MQ=255)
                                                                   ggattCATCGTTGGGGTCGGTTATAAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATgg                                                                                                                                                                                >  1:64105/4‑251 (MQ=255)
                                                                           cGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGAGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGa                                                                                                                                                                         <  1:139373/250‑1 (MQ=255)
                                                                                    cGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAg                                                                                                                                                               <  2:171667/251‑1 (MQ=255)
                                                                                        tataAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTAc                                                                                                                                                            <  1:110908/250‑1 (MQ=255)
                                                                                              ttCTGATTAGCCAGGGAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGCCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGCACGGCCGGTACAgtgtcc                                                                                                                                                     >  2:315031/1‑246 (MQ=255)
                                                                                                  gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTACCAGCACCACGGAG‑G‑TGGTGAACACGGTGGGCTAAGAGAATCCGGATGAAGACGGGCGTTACAGCATGGATGGGGAGGACGGTCAGTACAgtggcaccct                                                                                                                                                 >  2:349461/1‑242 (MQ=255)
                                                                                                    ttAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGt                                                                                                                                               >  2:125431/1‑251 (MQ=255)
                                                                                                        ccgggTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCGTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAGCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTTAGTACGGTCAGTACAGCGCTATTCTGTtggt                                                                                                                                           <  2:481989/248‑1 (MQ=255)
                                                                                                         cAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggtt                                                                                                                                          >  2:163313/1‑235 (MQ=255)
                                                                                                          aGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttg                                                                                                                                         <  1:163313/251‑18 (MQ=255)
                                                                                                              aaCACAGTGTGATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGGACAgtgtcatcctgcaggttgacg                                                                                                                                      <  2:62126/250‑21 (MQ=255)
                                                                                                                 acaGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctgcaggttgacggttt                                                                                                                                  >  1:448740/1‑227 (MQ=255)
                                                                                                                                   gccgGTTCAGGCGGGTTTTTTTGGGGGGGGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCAAAGGAAAACCGGGACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGGGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGATTACGGTCATTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCgg                                                                                                                <  2:447809/251‑1 (MQ=255)
                                                                                                                                       gtccAGGCGGGCTTTTTTGGGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCGCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCGGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGAcc                                                                                                            <  1:537475/248‑1 (MQ=255)
                                                                                                                                         gcgggtgggCTTTTGTGTGGGGGGAATAGGTCAGGAAAGATTTCAGTTGTACTTAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCGGCTGAAAGAAAAACGAAACAGCCCCA‑GGTGTG‑GGGTGAACGCGCTGGCTTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTAGGGTCATTTCAGCGGGATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAt                                                                                                          <  1:34094/245‑1 (MQ=255)
                                                                                                                                            ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                                                                                       >  1:320049/1‑251 (MQ=255)
                                                                                                                                               gggCTTTTTTGTGGGGTGAAGATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGAAAAACGTAACAACACCACGGTG‑G‑TGGTTAACACGCTGGCCTCGGAAAATCCGGATGAACCCGGGCGTGACAGCATGGACGTTGAGTAGGGTAAGTACAGCGTTATTCTGTTGGTGGAAGTATTCCCGCCGTCACATGCCGGGACCATTACCgt                                                                                                    <  1:369497/251‑1 (MQ=255)
                                                                                                                                                ggCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGCCCATTACCgtg                                                                                                   >  2:744/1‑251 (MQ=255)
                                                                                                                                                 gCTTTTTTGTGGGGTGAATATGGCAGTAAAGCTTTCAGGTGTACTGAAAGACGGCCCAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑AGGCGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTACCATGCCGGGACCATTCCCgtgt                                                                                                  >  2:284960/1‑251 (MQ=255)
                                                                                                                                                      tttGTGGGGTGAATATGGCAGTAAATATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCTCCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGTGCTCATAGAATCCGGATGAAGCCGGGCGGTACAGCATGGATGTGTTGTACGGTCAGTACAGTGTCATCCTGCGGGGTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGa                                                                                              <  1:43405/250‑1 (MQ=255)
                                                                                                                                                       ttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGaa                                                                                             <  1:366256/250‑1 (MQ=255)
                                                                                                                                                          tgGGGTGAATATGGCAGTAAAGATTTCAGGTGTCCTGAAAGCCGGAACAGGAACACAGGCACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCAAGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACGGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCAAATGCCGGGACCATTACCGTGTATGAAGAtt                                                                                         >  2:561195/1‑251 (MQ=255)
                                                                                                                                                          gggggTGAATAGGGCAGTAAAGATTTCAGGTGGACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAAACCAGATGAAAGCAAAACGTAACAGCACCGCGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAAACCGGTTGAAGCCGGGCGTTACAGCATGGAGGTGGAGTACGTTCAGTACAGCGTTATTTTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGGATGAAGAtg                                                                                         <  2:182601/250‑2 (MQ=255)
                                                                                                                                                            gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACCGGAAAACCGGTAAAGAACTGCACCATTCAGCTGCAAGCCAGACGTAACAGCCCCACGGTG‑A‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGCGTACGGCAAGTACAGTGTCATCCTGCCGGCTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGCTtc                                                                                        >  2:493375/1‑250 (MQ=255)
                                                                                                                                                            gggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCAACATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGCTGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCAACGTGTATGAAGATtcc                                                                                       >  2:502073/1‑250 (MQ=255)
                                                                                                                                                              gTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAACACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGCAGATtctc                                                                                      >  1:26086/1‑250 (MQ=255)
                                                                                                                                                               tGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGACAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATtcaca                                                                                     >  1:502116/1‑247 (MQ=255)
                                                                                                                                                                   aatGGCTGTAAAGATTGCAGGTTAACTGAAAGCCGGCACAGGAAAACCGGTACAGAACTGCAATACCCAGCGGAAAGCAAAACTGAAAAGAACCACGGTG‑G‑GGGTGAACACGCGGGCCTCAGAAAATCCGGAGGAAGCCGGGCGTTACAGCATTGACGTGGAGTCGGGTCAGTCCAGCGTTATTCTTTTGGTGAAATTATTCCCGCCGGCACATGCCGGGACCATTACCGTGTATGAAGATGCTCAACCCg                                                                                <  2:510030/250‑1 (MQ=255)
                                                                                                                                                                        cAGTAAAGATTTCAGGTGTCCTGAAAGACGGCACAGGAAATCCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTG‑G‑TGGTGACCCCGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTTTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACg                                                                           <  1:137197/251‑1 (MQ=255)
                                                                                                                                                                        cAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAAccggggacg                                                                           <  2:198451/251‑8 (MQ=255)
                                                                                                                                                                           tAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTg                                                                        >  2:156236/1‑251 (MQ=255)
                                                                                                                                                                            aaaGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCAGCGCACGCCGGGACCATCACCGGGTATGAAGATTCACAAccgtggacgctga                                                                       <  2:88738/251‑12 (MQ=255)
                                                                                                                                                                               gATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATg                                                                    >  1:180217/1‑251 (MQ=255)
                                                                                                                                                                                   tCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAttt                                                                >  1:517741/1‑251 (MQ=255)
                                                                                                                                                                                    caggagtccTGAAAGACGGCACAGACAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCCGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGAtttt                                                               >  2:258665/9‑251 (MQ=255)
                                                                                                                                                                                           ccTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCtgtg                                                        >  1:138193/2‑251 (MQ=255)
                                                                                                                                                                                              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGc                                                       <  1:175182/249‑1 (MQ=255)
                                                                                                                                                                                              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCa                                                     <  2:446098/251‑1 (MQ=255)
                                                                                                                                                                                              gAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGTTGAACACGCTGGCCTCATAAAATCCGGAGGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGTTCAGTACGGCGTTATTCTGTTGGGGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGcc                                                      <  2:341157/250‑1 (MQ=255)
                                                                                                                                                                                                 aGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATg                                                   >  2:161128/1‑250 (MQ=255)
                                                                                                                                                                                                    cGGCACAGGCAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGCTTTTCTCGGTGCCATGACgg                                               >  1:454630/1‑251 (MQ=255)
                                                                                                                                                                                                    cGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGCCgg                                               >  1:220903/1‑251 (MQ=255)
                                                                                                                                                                                                     ggCAAAGGAAAACCGGTACAGAACTGCACAATCCAGATGAAAGAAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTAAAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACAGGCAGGGACCATTAACGTGTATGAAGATTCTCAACCCGGTAAGCTGAATGATTTTAACGGTGCCATGACgg                                               >  1:279228/1‑250 (MQ=255)
                                                                                                                                                                                                       cacaGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACggagg                                            >  1:19581/1‑251 (MQ=255)
                                                                                                                                                                                                                  cGGTAAAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCCCCACGGTG‑G‑TGGTGAACACTTGGGCCTCAGAAAAGCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCGGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGt                                 <  2:320049/251‑1 (MQ=255)
                                                                                                                                                                                                                   ggTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCggc                                >  2:166983/1‑249 (MQ=255)
                                                                                                                                                                                                                            cggccaaaTCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGTGCTTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCAc                       <  1:200182/245‑1 (MQ=255)
                                                                                                                                                                                                                            cTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTG‑G‑TGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCCCGGCCGGAggtgc                       >  1:56678/1‑248 (MQ=255)
                                                                                                                                                                                                                               cacaATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGCTACGCTGAATGATTTTATCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTgc                    >  2:4505/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGt               >  1:191609/1‑251 (MQ=255)
                                                                                                                                                                                                                                     ccAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGtt              >  1:188442/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaGCAAAACGTAACAGCACCACGGTG‑G‑GGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCTCTGCGCCGTTTTGCACTg      >  2:456390/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aaaaCGTAACAGCACCACGGTG‑G‑TGGTGCACACGCTGGCCTAAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGAtgg  >  2:470524/1‑251 (MQ=255)
                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                          
CCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTATAAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTG‑G‑TGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACTGCGCCGTTTTGAACTGATGG  >  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: ‑