Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,429,1140CT50.9% 9.1 / 67.5 53I22I (ATC→ATT) stfRputative prophage side tail fiber protein StfR
Reads supporting (aligned to +/- strand):  ref base C (11/15);  new base T (17/10);  total (28/25)
Fisher's exact test for biased strand distribution p-value = 1.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACT  >  NC_000913/1428870‑1429356
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
tGGCGTGCCGTTTCACCCGACCGAATCCGTGGACATTGATATTGCTTATGAAGTCCCCGGCCGTGGCGACTGGCGCACTGACGGTTTAATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTTTGACAGCCCGCCGGGTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGGAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCt                                                                                                                                                                                                                                               <  2:240301/250‑1 (MQ=255)
             aaCCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGAGGGGTGAATATGGCAGAAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAACCCGGGACAGAACTGCACACTCCAGCTGAAAGCacaaact                                                                                                                                                                                                                                  >  2:84012/1‑244 (MQ=255)
                  gaccgaATACGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCCCTGCCGGTTTCTTCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGAAGAAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAACACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGc                                                                                                                                                                                                                            >  2:372804/1‑251 (MQ=255)
                        aTCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACg                                                                                                                                                                                                                      >  2:158547/1‑251 (MQ=255)
                            gtcgttgttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACggtgg                                                                                                                                                                                                                  >  1:440190/8‑251 (MQ=255)
                            gTGGCCCTTGATATTGCTTATGAAGGCCCTGGCAGGGGCGAGTGGCGCCTTGACGGTTTATTCGTGGGTGTCG‑TTATTAAGTTCTGATTAGCCAGGTAACAGAGTGTTTTGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATCCAGCTTAAAGCAAAACTTAACAGCACCACggtgg                                                                                                                                                                                                                  <  2:218611/251‑1 (MQ=255)
                              cgttgttGATATTGCTTATGAAGGCGCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGGCGTAACAGCACCACggtggt                                                                                                                                                                                                                 <  2:75240/245‑1 (MQ=255)
                                ccATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACggtggtggt                                                                                                                                                                                                              <  1:87927/251‑1 (MQ=255)
                                   ttGATATTGCTTATGAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGa                                                                                                                                                                                                            <  1:279030/250‑1 (MQ=255)
                                                 gAAGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTc                                                                                                                                                                                               >  2:130525/1‑249 (MQ=255)
                                                  aaGGCTCCGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAgg                                                                                                                                                                                            >  2:146274/1‑250 (MQ=255)
                                                   aGGCCCCGGCAGTGGCGCCTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGaa                                                                                                                                                                                           <  1:374511/251‑1 (MQ=255)
                                                    ggctccGGCAGTGGCGACTGGCGTACTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGaga                                                                                                                                                                                          <  1:129823/247‑3 (MQ=255)
                                                                       ggtgcgctgccGGTTTCATTGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCggg                                                                                                                                                                         <  2:386607/239‑1 (MQ=255)
                                                                          gtactaaCGGCTTCCTCGTTGGAGTCGGTTAT‑CAATTCTCATTAGCCAGCACACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTACAGATTTCAGGAGTCCTGACAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGCAAGACAGACGTAACAGCACCACGGAGGTGGTGACCACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGtt                                                                                                                                                                     >  2:192610/7‑250 (MQ=255)
                                                                            aCTGACGGATTCATCGTTGGGGTCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTAc                                                                                                                                                                   >  1:124689/1‑250 (MQ=255)
                                                                                acggattcatcgttgggggaggTTAT‑AAATGCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTGTTGTTTGGGGGGAATATGGCAGTAAAGATTTCAGGAGTCCGGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTATCAGCACCACGGTGGTGGGGAACAGGGTGGGCTCATAGAATCCGTAGGAAGTCGGGCGTTACAGCAt                                                                                                                                                              <  2:406717/231‑1 (MQ=255)
                                                                                   gTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGGGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCGGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGTCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGa                                                                                                                                                           <  1:393765/251‑1 (MQ=255)
                                                                                     ttCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACg                                                                                                                                                         >  1:102782/1‑251 (MQ=255)
                                                                                            ttggggtCGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGt                                                                                                                                                   >  1:424676/6‑250 (MQ=255)
                                                                                             tggggtaggTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACCGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAACGACGGCACCGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCCGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGACGCCGGGCGTTACAGCATggttgtggcgtac                                                                                                                                                 >  2:182395/8‑240 (MQ=255)
                                                                                                   cGGTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCa                                                                                                                                            <  2:6745/250‑1 (MQ=255)
                                                                                                     gTTAT‑AAATTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGt                                                                                                                                          >  1:209510/1‑250 (MQ=255)
                                                                                                           aaGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTGCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGGAACAGCACCACGGGGGGGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGTCGTTACAGCATGGATGTTGAGTACTGTCAGTTCAGCg                                                                                                                                    <  2:210163/251‑1 (MQ=255)
                                                                                                            aattCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGGAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGat                                                                                                                                                          >  1:98805/3‑227 (MQ=255)
                                                                                                            aattCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGGAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGat                                                                                                                                                          <  2:98805/226‑2 (MQ=255)
                                                                                                                 tGATTAGCCAGGTAACATAGTGTTAAGACAGGCCGCCGGTTCAGGGAGGCCTTTTTGTGGGGTGAATATGGCAGCAAAGATTTCAGGTGTACTGAAAGACGGCACCGGAAAACCGGTACAGCACTGCACACTCCAGCTGAAAGCAAACCGTAACAGCACCCAGGTGGGGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCAGGCAGTGCAGAATGGAAGTTGAGTACGGTCAGTACCGCGTTATTc                                                                                                                              >  2:293170/1‑251 (MQ=255)
                                                                                                                  gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTc                                                                                                                              <  2:4102/250‑1 (MQ=255)
                                                                                                                  gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTc                                                                                                                              >  1:4102/1‑250 (MQ=255)
                                                                                                                  gATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgacatcct                                                                                                                             <  1:121974/251‑10 (MQ=255)
                                                                                                                   aTTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAgtgtcatcctg                                                                                                                            >  1:317464/1‑241 (MQ=255)
                                                                                                                        ccAGGTAACACAGTGTTATGACAGCCCGCCGGAACCGGTGGGCTTTTTTGTGGGGGGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATATGGAGTACGGTCAGTACAgtgtcatcctgcaggt                                                                                                                       >  1:341575/1‑236 (MQ=255)
                                                                                                                                                           aGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCAtt                                                                                     >  1:325183/1‑250 (MQ=255)
                                                                                                                                                            ggCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGTGTTACAGCATGGACTTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGTAAGGATTCCCGCCGTCACATGCCGGGACCATTAc                                                                                   <  1:425325/251‑1 (MQ=255)
                                                                                                                                                                ggtttttgtGTGGGGGGAAAATGTCAGTAAAGATTTCAGGTGGACTGAAGGGCGGTACAGGAAAACCGGTACAGAACTGCACAATCCAGCGGAAAGCAAAACGTAACAGCCCCGTGGTGGGGGTAAACTCGCGGGCCTCAGAAAATCCGGATGAAGCCGGGCTTTCCGGCATGGTCGTTGAGTGCGTTCAGTCCAGCGTTATTTTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGTCCATTACCgtg                                                                               <  2:289407/248‑1 (MQ=255)
                                                                                                                                                                    tttttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGGTGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATg                                                                           >  1:86619/1‑251 (MQ=255)
                                                                                                                                                                    ttttgggggggTGAATATGGAAGTAAAGATTTCAGGGGTCCGGAAGGACGGCACAGGAAAACCGGTCCAGAACTGCCCAACCCAGCTGAAAGAAAAACGTAACAGCACCGCGGGGGTGGTGAACCCGCTGGCCTCAGAAAATCCGGATGAAGCGGGGCGTTACGGCATGGACGTTGAGTCCGGTCATTAAAGCGTTTTTCTTTTGGTGGAGGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTAt                                                                            <  2:111399/243‑1 (MQ=255)
                                                                                                                                                                      tttGTGGGGTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGa                                                                          >  1:111536/1‑250 (MQ=255)
                                                                                                                                                                        tgaggggTGAATATGGCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGa                                                                          >  2:133227/4‑248 (MQ=255)
                                                                                                                                                                                       gCAGTAAAGATTTCAGGAGTCCTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTGCCACCATCGCACGCCGGGACCATCACCGTGGATGAAGATTCACAAccggggac                                                        <  1:130525/251‑7 (MQ=255)
                                                                                                                                                                                            aaaGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTg                                                    >  1:434348/1‑250 (MQ=255)
                                                                                                                                                                                                            cTGAAAGCCGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGAACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGc                                   >  1:297843/1‑251 (MQ=255)
                                                                                                                                                                                                                  gACGGCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGAc                             <  1:470896/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggCACAGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACgg                           >  1:151270/1‑250 (MQ=255)
                                                                                                                                                                                                                       cacaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGt                                                                                                                                          >  2:298591/1‑137 (MQ=255)
                                                                                                                                                                                                                       cacaGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGt                                                                                                                                          <  1:298591/137‑1 (MQ=255)
                                                                                                                                                                                                                        acaGGAAAACCGGTACAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACggagga                       <  2:317464/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGAc                                                                                                                                                          >  1:453754/1‑110 (MQ=255)
                                                                                                                                                                                                                                  cGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGGAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGGTACAGCATGGAc                                                                                                                                                          <  2:453754/110‑1 (MQ=255)
                                                                                                                                                                                                                                       cAGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGACCATCACCGTGTATGAAGATTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg         >  1:52600/1‑250 (MQ=255)
                                                                                                                                                                                                                                        aGAACTGCACCATTCAGCTGAAAGCCAGACGTAACAGCACCACGGTGGTGGTGAACACGGTGGGCTCAGAGAATCCGGATGAAGCCGGGCGTTACAGCATGGATGTGGAGTACGGTCAGTACAGTGTCATCCTGCAGGTTGACGGTTTTCCACCATCGCACGCCGGGGCCATCACCGTGTATGAAGCTTCACAACCGGGGACGCTGAATGATTTTCTCTGTGCCATGACGGAGGATGATGCccggccgg         >  2:279879/1‑249 (MQ=255)
                                                                                                                                                                                                                                        aGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGTTTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAg       <  1:158547/251‑1 (MQ=255)
                                                                                                                                                                                                                                             gacacaACCCAGTTGAAAGAAAAACGTAACGGCACCACGGTGGTGGTGAATCCTCTGGCTCCAGAAAAACCGGATGAATCCGGGCGTTACAGCATGATCGTTGAGTACGTTCATCATAGCGTTTTTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACt  <  2:72644/249‑1 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                  
TGGCGTGCAGTTTAACCCGACCGAATCCGTGGCCATTGATATTGCTTATGAAGGCCCCGGCAGTGGCGACTGGCGCACTGACGGTTTCATCGTGGGTGTCGGTTAT‑AAGTTCTGATTAGCCAGGTAACACAGTGTTATGACAGCCCGCCGGTTCAGGCGGGCTTTTTTGTGGGGTGAATATGGCAGTAAAGATTTCAGGTGTACTGAAAGACGGCACAGGAAAACCGGTACAGAACTGCACAATCCAGCTGAAAGCAAAACGTAACAGCACCACGGTGGTGGTGAACACGCTGGCCTCAGAAAATCCGGATGAAGCCGGGCGTTACAGCATGGACGTTGAGTACGGTCAGTACAGCGTTATTCTGTTGGTGGAAGGATTCCCGCCGTCACATGCCGGGACCATTACCGTGTATGAAGATTCTCAACCCGGTACGCTGAATGATTTTCTCGGTGCCATGACGGAGGATGATGCCCGTCCGGAGGCACT  >  NC_000913/1428870‑1429356

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