Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,104,614	0	T	C	38.0%	37.0 / 50.2	50	K121R (AAA→AGA)	cpxA	sensor histidine kinase CpxA
Reads supporting (aligned to +/- strand):  ref base T (19/12);  new base C (6/13);  total (25/25)
Fisher's exact test for biased strand distribution p-value = 7.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTATCCAGAAGCTCGGTCATCTGGCGTGAATCGA  >  NC_000913/4104357‑4104887                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  ggCGTAAGTTTCCCTGAGCAACTTCATCGGAAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTACTTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCAATACTTTTTCTTCTGCTGATGATCGGCGTTATCGGCCt                                                                                                                                                                                                                                                     <  1:100196/288‑1 (MQ=255)                           tCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCCGTTCCACGCGGCCATACTTTCTCTTCTGCGga                                                                                                                                                                                                                                                                        >  1:474019/1‑243 (MQ=255)                                                             ccAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCtt                                                                                                                                                                                         <  1:524938/288‑1 (MQ=255)                                                               aGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCg                                                                                                                                                                                       <  1:296714/288‑1 (MQ=255)                                                                nncgcaaacacAAACACAACAGAACAGGCGTACTGACAAACATGGGGACAATCAGCAGTAATCGCGGGCGGTAAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCCGATAAAGGTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCCGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCAcg                                                                                                                                                                                                <  2:88663/268‑1 (MQ=255)                                                                 aCTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGCCGTTATCGGACTGACCAATAAAGTTACGAATGATCTGCATTTCTCTGCGTCCAGCGCCGACCACCCGGCCTTCggt                                                                                                                                                                                     >  1:595408/1‑288 (MQ=255)                                                                  cTCCAGGCCAACCACAACAGCATCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTCTCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGa                                                                                                                                                                                                                                    <  1:277876/240‑1 (MQ=255)                                                                    ccAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGAt                                                                                                                                                                                                                               <  1:343817/243‑1 (MQ=255)                                                                             ccACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCgg                                                                                                                                                                                                                                                                 <  2:182864/200‑1 (MQ=255)                                                                             ccACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCgg                                                                                                                                                                                                                                                                 <  1:518311/200‑1 (MQ=255)                                                                             ccACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCgg                                                                                                                                                                                                                                                                 >  1:182864/1‑200 (MQ=255)                                                                             ccACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCgg                                                                                                                                                                                                                                                                 >  2:518311/1‑200 (MQ=255)                                                                             ccACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCGCGCGGCCTTCGGTGGTCACCAATaa                                                                                                                                                                         >  1:299179/1‑288 (MQ=255)                                                                             ccACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATaa                                                                                                                                                                         >  1:196039/1‑288 (MQ=255)                                                                                      cagcGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGTCATACTTTTTCTTCTTCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATGACGTGTCCTTCTGTGGTACCCAAAAACAAACGcgg                                                                                                                                                                >  1:272155/1‑286 (MQ=255)                                                                                        gcggcgTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCg                                                                                                                                                                                       >  1:88663/1‑263 (MQ=255)                                                                                                      caCATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTTGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGc                                                                                                                                                                  <  1:383032/269‑1 (MQ=255)                                                                                                        cATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGACGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTTGGTGGTCACCAATaa                                                                                                                                                                         >  1:151100/1‑261 (MQ=255)                                                                                                             tGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGa                                                                                                                                                                                                                                    <  2:544816/197‑1 (MQ=255)                                                                                                             tGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGa                                                                                                                                                                                                                                    >  1:544816/1‑197 (MQ=255)                                                                                                                       nggtgagggggagcgtTCAAACAGAAAGTTAAGGAAATCGGATGGGGAACTGCTGGCCGGACGGATCAGATAAAGTTGGTAATTATCGTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTGCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGcc                                                                                                                                    <  2:196039/271‑1 (MQ=255)                                                                                                                                    cggTCAAACAGTAAGTTAATTAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATa                                                                                                                                                                                                                                             <  2:528353/165‑1 (MQ=255)                                                                                                                                    cggTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTa                                                                                                                                                                                                                                                            <  1:262953/150‑1 (MQ=255)                                                                                                                                    cggTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTa                                                                                                                                                                                                                                                            >  2:262953/1‑150 (MQ=255)                                                                                                                                    cggTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATa                                                                                                                                                                                                                                             >  1:528353/1‑165 (MQ=255)                                                                                                                                        cAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACg                                                                                                                                                                   <  1:550792/235‑1 (MQ=255)                                                                                                                                              gTAAGTTAATGAAATCGGATGGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTCTCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATGTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGcc                                                                                                                                                   <  1:294084/245‑1 (MQ=255)                                                                                                                                                                                   aCGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCAATACTTTTTCTTCTGCGGATTATCGGCGTTATCGTCCTGACCAATACAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCGCGCGGCCTTCGGTGGTCACCAATAACAACCGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCCTcaccat                                                                    >  1:587034/1‑282 (MQ=255)                                                                                                                                                                                        tCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGTGGCCTTCGGTGGTCACCAATAACAAACg                                                                                                                                                                   >  1:539439/1‑187 (MQ=255)                                                                                                                                                                                        tCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGTGGCCTTCGGTGGTCACCAATAACAAACg                                                                                                                                                                   <  2:539439/187‑1 (MQ=255)                                                                                                                                                                                        tCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATgctgct                                                              >  1:305495/1‑288 (MQ=255)                                                                                                                                                                                          aGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCAAGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGAACAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCa                                                            <  1:451816/288‑1 (MQ=255)                                                                                                                                                                                                      tAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCAcgcg                                                                                                                                                                                              >  1:129241/1‑146 (MQ=255)                                                                                                                                                                                                      tAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCAcgcg                                                                                                                                                                                              <  2:129241/146‑1 (MQ=255)                                                                                                                                                                                                               tCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGAGCACGTGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTTGGCGGATCGTTCGCCAGCTCCGCTTCGACAt                                                                    >  1:214482/1‑259 (MQ=255)                                                                                                                                                                                                                        cgcacgGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGTGGCCTTCGGTGGTCACCAATaa                                                                                                                                                                         >  1:71061/1‑149 (MQ=255)                                                                                                                                                                                                                        cgcacgGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGTGGCCTTCGGTGGTCACCAATaa                                                                                                                                                                         <  2:71061/149‑1 (MQ=255)                                                                                                                                                                                                                             ggAGAACGGACCGACCAGTTCCACGCGGTCCTACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCGCGCGGCCTTCGGTGGTCACCAATAACACACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGCGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTATCCAGAAg                         <  1:407707/288‑1 (MQ=255)                                                                                                                                                                                                                                         gaccAGTTCCACGCGGCCATACTTTTTCTTCTGCGGCTGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCCATAACAAACGCTGTCCTGGCGGGGCCCACTTATCAATCGCCCGGAACAGACTCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAACCATCAGCCCCTGACGCTGTTCGCTATCCAGAAGCTCGTTCATCTg             >  1:110968/1‑288 (MQ=255)                                                                                                                                                                                                                                         gaccAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGCCGGACCGTTCGCCAGCCCCGCTTCCACATGCTACCCAATAACCAGCCACTCGCGCCGGTCGcagatctc                            >  2:67900/1‑266 (MQ=255)                                                                                                                                                                                                                                         gaccAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGtt                                                                                                                                                                                                              <  2:325884/95‑1 (MQ=255)                                                                                                                                                                                                                                         gaccAGTTCCACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGtt                                                                                                                                                                                                              >  1:325884/1‑95 (MQ=255)                                                                                                                                                                                                                                           ccAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATTATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCa                                                                                                                                                                             <  2:381302/126‑1 (MQ=255)                                                                                                                                                                                                                                           ccAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCa                                                                                                                                                                             >  1:381302/1‑126 (MQ=255)                                                                                                                                                                                                                                           ccAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGATCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTACCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTTGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTCTCCAGAAGCTCGGTCATCTGGc           >  1:580078/1‑288 (MQ=255)                                                                                                                                                                                                                                                 ccACGCGGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGAGCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAAGCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTATCCAGAAGCTCGGTCATCTGGCGTGaa      <  1:415211/287‑1 (MQ=255)                                                                                                                                                                                                                                                    cgcgGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGTGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGcc                                                                                                                                    <  2:557161/158‑1 (MQ=255)                                                                                                                                                                                                                                                    cgcgGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGTGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGCCGGTGCCCACTTATCAATCGcc                                                                                                                                    >  1:557161/1‑158 (MQ=255)                                                                                                                                                                                                                                                    cgcgGCCATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTATCCAGAAGCTCGGTCATCTGGCGTGAATCGa  >  1:445288/1‑288 (MQ=255)                                                                                                                                                                                                                                                          cATACTTTCTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTATCCAGAAGCTCGGTCAt                <  1:264163/268‑1 (MQ=255)                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                  GGCGTAAGTTTCCCTGGGCAACTTCATCGGCAGCGTTTTTCAGCTTACGCGCCGGTTTTGCCAGACTCCAGGCCAACCACAACAGCAGCGGCGTACTGACCAACATGGTGACAATCAGCAGTAATAGCGGGCGGTCAAACAGTAAGTTAATGAAATCGGATTGAGAACTGCTGGCCGGACGAATCAGATAAAGTTGGTAATTATCTTCGCCATCACGCACGGAGAACGGACCGACCAGTTCCACGCGGCCATACTTTTTCTTCTGCGGATGATCGGCGTTATCGGCCTGACCAATAAAGTTACGAATGATCTGCATTTCGCTGCGTTCAGCGCCGATCACGCGGCCTTCGGTGGTCACCAATAACAAACGCTGTCCTGGCGGTGCCCACTTATCAATCGCCCGGAACAGACGCCGCCACCACATTAAATCGTTGGGCGGATCGTTCGCCAGCTCCGCTTCGACATGCTGCTCAATCATCAGCCCCTGACGCTGTTCGCTATCCAGAAGCTCGGTCATCTGGCGTGAATCGA  >  NC_000913/4104357‑4104887

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 10 ≤ ATCG/ATCG < 22 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑