Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913578,4270GC50.0% ‑4.3 / 106.0 58T107T (ACG→ACC) 
R34P (CGG→CCG) 
rzpD
rzoD
putative prophage endopeptidase RzpD
putative prophage lysis lipoprotein RzoD
Reads supporting (aligned to +/- strand):  ref base G (14/15);  new base C (15/14);  total (29/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.76e-01
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

TCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATT  >  NC_000913/578178‑578671
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
tCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGCACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGCTGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACgg                                                                                                                                                                                                                                                     >  2:97021/1‑251 (MQ=255)
 cATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACgg                                                                                                                                                                                                                                                     <  2:125666/250‑1 (MQ=40)
       cGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCgg                                                                                                                                                                                                                                              >  1:334592/1‑251 (MQ=40)
               cGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGGGGAt                                                                                                                                                                                                                                       >  2:12242/1‑250 (MQ=255)
                                  cAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGAcg                                                                                                                                                                                                                   <  2:331034/251‑2 (MQ=21)
                                           aaaaaTGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGacac                                                                                                                                                                                                          <  1:273786/251‑1 (MQ=34)
                                             aaaggCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGAGGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCTTTGATAATGCAGCCTCCCCTCGACTGGCATACACCg                                                                                                                                                                                                        <  2:253894/247‑1 (MQ=255)
                                              aaTGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGc                                                                                                                                                                                                       >  1:341607/1‑251 (MQ=34)
                                              aaTGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGc                                                                                                                                                                                                       >  1:229576/1‑251 (MQ=34)
                                                               gCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCa                                                                                                                                                                                      <  2:115968/251‑1 (MQ=34)
                                                                      aaCGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCag                                                                                                                                                                               <  2:346381/251‑1 (MQ=33)
                                                                         gcgGCAATTACTGACAAGCAGATGCGTCAGCGTGATTTTGCTGCGCGCGATGTAAAATACACGAAGGAGTTAGCTGAGGCTAACGCTGAAATTGATGCTCTGCGGG‑AGGATGTTGCCGTGGGTCGTCGTCGGTTGCACACCAAAGCAGTCTGTCAGTCAGTGGGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCagaga                                                                                                                                                                            <  2:210515/251‑1 (MQ=11)
                                                                             cAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCACTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAgg                                                                                                                                                                        >  2:74002/1‑251 (MQ=37)
                                                                              aaTTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGc                                                                                                                                                                       >  1:58571/1‑251 (MQ=18)
                                                                                     gACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATc                                                                                                                                                                  <  1:197052/249‑1 (MQ=40)
                                                                                          gCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATg                                                                                                                                                            <  1:106875/250‑1 (MQ=255)
                                                                                            aGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGCCTGGCAGACACCGCTGAACGGGATTATTTCACCCTCCGAGAGAGGCTGATCACTATGCa                                                                                                                                                          >  1:40816/1‑250 (MQ=255)
                                                                                             gATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGAGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCaa                                                                                                                                                         >  1:209823/1‑250 (MQ=255)
                                                                                                 cGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaac                                                                                                                                                     >  1:142598/1‑250 (MQ=255)
                                                                                                   tCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGGCACCGCTGAACGGGGTTATTTCACCCTCAGAGAGAGGCTGCTCACTATCCAAAaacaac                                                                                                                                                  >  1:23287/1‑251 (MQ=255)
                                                                                                    cAGCGTGATGTGGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTTTCAGTCAGTGCTTGAAGCCACCACGGCCTCCGGCGTGGATAATGCGGCCTCCCCCCGACTGGCAGACACCGCTGGACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAaacaac                                                                                                                                                  <  1:266872/250‑1 (MQ=255)
                                                                                                         tGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTgga                                                                                                                                              >  1:332459/1‑249 (MQ=255)
                                                                                                                    cgcTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCAGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGCAGGACCCCAGa                                                                                                                                   >  2:335496/1‑249 (MQ=255)
                                                                                                                    cgcTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCAGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAg                                                                                                                                 <  1:335496/251‑1 (MQ=255)
                                                                                                                          atacAAAATACACGCAGGATTTAGCTGATGCTAAAGCTGAAACTGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGACCGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAACAACAACTGGACGGAACCCAGAAGCATAtt                                                                                                                           >  2:314973/4‑251 (MQ=14)
                                                                                                                           tGCAAACTAGACGAAGGAGTTAGCTGATTAGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCCAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGCGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGCTTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACACAGAAGTATAtt                                                                                                                           >  1:112360/1‑250 (MQ=17)
                                                                                                                            gCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAAGGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACGGGAAGGAACCCAGAAGTATATTa                                                                                                                          <  1:97021/250‑1 (MQ=255)
                                                                                                                                aaTACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACt                                                                                                                                                 <  1:14831/223‑1 (MQ=255)
                                                                                                                                aaTACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACt                                                                                                                                                 >  2:14831/1‑223 (MQ=255)
                                                                                                                                aaTACACGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGACGGAACCCAGAAGTATATTAATg                                                                                                                       >  2:346826/1‑249 (MQ=255)
                                                                                                                                   acactaaggatttagcggagtctaaaGCTGAAAATGATGC‑CTTGGTGGATGATTTTGCCGCTGTGCGTCGTCGGTTGCACATCAAAGCAGTCTTTAAGGCAGTGCGTGGAGCCACCCCCGCCGCCGGCGTGGATAATGCAGCCTCCTCCCGACTGGCAGACACCGCTGAACGGGATTATTTCAGCTTCAGAGATAGTCTGATCATTAGGCAAAAACATCTGGAAGGTATCCAGAAGTATATTAATGAGCAg                                                                                                                  <  2:362248/228‑1 (MQ=9)
                                                                                                                                   acacGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGGG‑ATGATGTTGCCGCTGGTCGTCGGCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGAGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTAGGTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCa                                                                                                                   <  2:209823/250‑1 (MQ=17)
                                                                                                                                      cGAAGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCCAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTg                                                                                                                >  2:320971/1‑250 (MQ=255)
                                                                                                                                         aGGAGTTAGCTGATGCTAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGc                                                                                                               >  2:277479/1‑248 (MQ=255)
                                                                                                                                                aGCTGATGCGGAAGCTGAAAATGCTGCTCTGCGTG‑ATGCAGTTGCCGCTGGTCGACGTCGGATGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTACGGCGTGGATAATGCAGCATCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCCAAAACAACTGGAAGGAACCCAGACGTATATTAATGAGAAGTGCAGAAagaact                                                                                                     >  2:261114/1‑248 (MQ=18)
                                                                                                                                                   tGATGCGCAAGCTGAAAATGATGCTCTGCGTG‑ATGCTGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGAc                                                                                                  >  2:135889/1‑251 (MQ=255)
                                                                                                                                                         taaaGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATc                                                                                            >  1:368898/2‑251 (MQ=255)
                                                                                                                                                          aaaGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCg                                                                                           <  2:248887/251‑1 (MQ=255)
                                                                                                                                                                              gCGTG‑ATGGTGCGGCCGCTGGTCGTCGTCGGTTGCACGTCAAAGCAGGCGGTCAGTCAGTGCGTTAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCGGAACGGGATTAGTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGGATGAACCCAGAAGTATATTAATGAGCAGGGCAGATAGAGCTGATCATATCGTTGGGCTACTCATGCAatt                                                                        <  2:384313/250‑1 (MQ=18)
                                                                                                                                                                              gCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCaa                                                                          <  1:78115/248‑1 (MQ=255)
                                                                                                                                                                              gCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCaa                                                                          <  1:277479/248‑1 (MQ=255)
                                                                                                                                                                                       atgTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGc                                                               >  1:143784/1‑251 (MQ=255)
                                                                                                                                                                                          ttGCCGCTGAGCGTCCTCGGATGCACATCAATGCACTCTGTCAGTCAGCCCGCCAAGCCCACACCGCCTCCAGCGTCGATAATGCAGCCTCCGCCCGACTGGCAGCCAACGCCGACCGGGATTATTTCACCATCACAGAGAGGATGATCACTATGCAAAAACAAATGGAAGGCACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATCTCGCAGGGCAACTCATGCAATTATTTTGAGCAAt                                                            >  2:196246/1‑251 (MQ=255)
                                                                                                                                                                                             ccGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGAATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATaca                                                         <  1:320971/251‑1 (MQ=255)
                                                                                                                                                                                               gCTGGTCGTCGTATGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATGATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGTTGACCATATCGATTGGCAACTCATTCAATTATTTTGAGCAATacaca                                                       <  2:33183/251‑1 (MQ=255)
                                                                                                                                                                                                       tcgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAGTGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTc                                               >  1:245302/1‑251 (MQ=255)
                                                                                                                                                                                                        cgtcgGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAGTGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATCCACACGCGCTTcc                                              >  1:258502/1‑251 (MQ=255)
                                                                                                                                                                                                             gTTGCACATCAAAGCAGTATGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAAAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCgg                                         >  2:282342/1‑251 (MQ=255)
                                                                                                                                                                                                                       aaaGCAGGCTGGCAGCCAGTGCGTGAAGCCACCACCGCCGCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTAGTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATg                               <  2:30470/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGc                              <  1:74002/251‑1 (MQ=255)
                                                                                                                                                                                                                        aaGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGc                              >  1:315289/1‑251 (MQ=255)
                                                                                                                                                                                                                           cAGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTa                           <  1:321626/251‑1 (MQ=255)
                                                                                                                                                                                                                            aGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGAcc                                                                                                 >  2:112788/1‑180 (MQ=255)
                                                                                                                                                                                                                            aGTCTGTCAGTCAGTGCGTGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGAcc                                                                                                 <  1:112788/180‑1 (MQ=255)
                                                                                                                                                                                                                                     gtcaggtcgtaaaGCCCCCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGTGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATaaa                 <  1:394290/240‑1 (MQ=255)
                                                                                                                                                                                                                                        agtGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAAcc              <  1:12242/251‑1 (MQ=255)
                                                                                                                                                                                                                                            cgggAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGc          <  2:368898/248‑1 (MQ=255)
                                                                                                                                                                                                                                              tGAAGCCACCACCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCa         >  1:127358/1‑250 (MQ=255)
                                                                                                                                                                                                                                                aaGCCGCCACGGCCTCAGGCGTGGATAATGCAGCCGCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCGCCGAGAGAGGCTTATCACGATGCAAAAGCAACTGGAAGGAACCCAGAAGTATATTAATGATCAGTTCAGATAGAGCTGACCATATCGATGGGCCAATCATGCAATTATTTTGAGCAATACACACTCGCTTCCAGCGGAATATAAAGGCTTAAAGTAATAAAACCGAGCAATc      <  2:270153/251‑1 (MQ=37)
                                                                                                                                                                                                                                                    caccacCGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCAtt  <  1:346826/251‑1 (MQ=255)
                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                    
TCATTACCGTGATAACGCAATCGCCTACAAAGTCCAGCGCGACAAAAATGCCAGAGAACTGAAGCTAGCGAACGCGGCAATTACTGACATGCAGATGCGTCAGCGTGATGTTGCTGCGCTCGATGCAAAATACACGAAGGAGTTAGCTGATGCGAAAGCTGAAAATGATGCTCTGCGTG‑ATGATGTTGCCGCTGGTCGTCGTCGGTTGCACATCAAAGCAGTCTGTCAGTCAGTGCGTGAAGCCACCACGGCCTCCGGCGTGGATAATGCAGCCTCCCCCCGACTGGCAGACACCGCTGAACGGGATTATTTCACCCTCAGAGAGAGGCTGATCACTATGCAAAAACAACTGGAAGGAACCCAGAAGTATATTAATGAGCAGTGCAGATAGAGCTGACCATATCGATGGGCAACTCATGCAATTATTTTGAGCAATACACACGCGCTTCCAGCGGAGTATAAATGCCTAAAGTAATAAAACCGAGCAATCCATT  >  NC_000913/578178‑578671

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