Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	309,985	T→A	G14G (GGA→GGT)	ecpB ←	putative fimbrial chaperone EcpB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	309,985	0	T	A	97.6%	149.6 / NA	42	G14G (GGA→GGT)	ecpB	putative fimbrial chaperone EcpB
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base A (17/24);  minor base C (1/0);  total (18/24)
Fisher's exact test for biased strand distribution p-value = 4.29e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01

AGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGT  >  NC_000913/309738‑310219                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           agaagaAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGATCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCg                                                                                                                                                                                                                                         >  1:20136/1‑251 (MQ=255)  gaagaaGCGGAGCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCgg                                                                                                                                                                                                                                        <  1:327310/251‑1 (MQ=255)              cACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCaga                                                                                                                                                                                                                             <  2:84311/250‑1 (MQ=255)                 ttCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGGCGGTACTGTTTTTGATCGTTTTGCTCCGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTTGGCCGGCGATTTACCGGAAAACAGCagaacg                                                                                                                                                                                                                          >  2:175057/1‑247 (MQ=255)                        gCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGGGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATAGCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGGGAGGCAg                                                                                                                                                                                                                  <  2:200791/251‑1 (MQ=255)                                    agcagcAAGCTGGCGGGAGTTAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGGGAAGAGAGCCGTTCGAGTCGGATATTGATAAGGCGACCGCTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAATCGATGATATATCGCCGACATCCAGCGCCTGGGCCGGATTTAGACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCttttt                                                                                                                                                                                                      <  2:129996/251‑1 (MQ=255)                                       agcaAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCAt                                                                                                                                                                                                   >  1:111559/1‑251 (MQ=255)                                        gcaAGCGGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATa                                                                                                                                                                                                  <  2:178127/251‑1 (MQ=255)                                               ggCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGGCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTTCTTTTTCATAATAATTc                                                                                                                                                                                           <  2:45888/251‑1 (MQ=255)                                                     agtgagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGCCGTCAAGCGGTGAAGAGAGCCGTTCGCGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATAGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATCCCGGAAAACAGCAGAGCGAGAGGCAGAAGGAGCTTTTTCATAACAATTCCGTccc                                                                                                                                                                                     >  1:296401/1‑250 (MQ=255)                                                      gtgagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAg                                                                                                                                                                                    <  1:243863/251‑1 (MQ=255)                                                        gagtAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAgg                                                                                                                                                                                   >  1:580/1‑250 (MQ=255)                                                            aGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATaa                                                                                                                                                                              <  2:111559/251‑1 (MQ=255)                                                                     tCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAggggg                                                                                                                                                                     >  1:75682/1‑251 (MQ=255)                                                                            tGTCCATTGAGATAACCTGCCCGTCGTCACGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTAccc                                                                                                                                                              >  1:91066/1‑251 (MQ=255)                                                                                       aTAACCTGCCCGACGTCAAGCGGTGCAGAGAGCCGTTCGAGACGGCTATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGAGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTAc                                                                                                                                                    >  2:153529/1‑250 (MQ=255)                                                                                                ccGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCagag                                                                                                                                             >  1:260642/1‑248 (MQ=255)                                                                                                  gtcgtcAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATACGCCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCAATGTTCATAAACGATGATATATAGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACACCTGCAGGTACATCAGAGAGAtt                                                                                                                                        >  2:130292/1‑251 (MQ=255)                                                                                                  gtcgtcAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGAtt                                                                                                                                        <  2:280106/251‑1 (MQ=255)                                                                                                        aaGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTgg                                                                                                                                  <  1:261952/251‑1 (MQ=255)                                                                                                          gCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTc                                                                                                                                <  2:252318/251‑1 (MQ=255)                                                                                                                        cGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGa                                                                                                                   >  2:221191/1‑250 (MQ=255)                                                                                                                          ttCGAGACGGATATTGCTAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGCGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACt                                                                                                                >  2:246449/1‑251 (MQ=255)                                                                                                                           tCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACt                                                                                                                <  2:94345/250‑1 (MQ=255)                                                                                                                                               gCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTcc                                                                                           <  1:121187/251‑1 (MQ=255)                                                                                                                                                          gggtccTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACAGCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATGCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGAGTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCGCCAGATGCCTTcc                                                                                <  2:301222/246‑1 (MQ=255)                                                                                                                                                           ggTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCg                                                                               >  1:141927/1‑251 (MQ=255)                                                                                                                                                              aCTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGt                                                                             >  2:37562/1‑250 (MQ=255)                                                                                                                                                                      tGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCt                                                                     <  2:88625/250‑1 (MQ=255)                                                                                                                                                                             tttGCTCAGCGTGCTGCTGTCCTTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATGACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTTCTGGTACATCAGAGAGATTAACTGGGCCAGGTCGCGTCGTACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGGGCTGTAATCt                                                              <  2:113609/250‑1 (MQ=255)                                                                                                                                                                                gCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGtt                                                           <  1:64286/250‑1 (MQ=255)                                                                                                                                                                                         gctgctGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCggtgg                                                 >  1:113688/1‑251 (MQ=255)                                                                                                                                                                                                 cACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtg                                         <  2:17987/251‑1 (MQ=255)                                                                                                                                                                                                 cACTGTTCATAAACGATGATATAGCGCCGACAACCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCAtt                                          <  1:161116/250‑1 (MQ=255)                                                                                                                                                                                                  aCTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATtgg                                        <  1:305904/251‑1 (MQ=255)                                                                                                                                                                                                         aTAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTAcc                                 <  2:120407/251‑1 (MQ=255)                                                                                                                                                                                                         aTAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTAcc                                 <  2:75682/251‑1 (MQ=255)                                                                                                                                                                                                                    tatatCGCCGACATCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATgg                      >  1:36730/1‑251 (MQ=255)                                                                                                                                                                                                                                aTCCAGCGCCTGGGCCGGAGATATACCGGAAAACAGCAGCGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAAc          >  2:171187/1‑251 (MQ=255)                                                                                                                                                                                                                                      cgcCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCt    >  1:211328/1‑251 (MQ=255)                                                                                                                                                                                                                                      cgcCTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGGCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACTCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGGTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCt    <  1:295660/251‑1 (MQ=255)                                                                                                                                                                                                                                        ccTGGGCCGGAGATATACCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCtgt  <  2:91066/251‑1 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                           AGAAGAAGCGGATCACTTCGCTGGCTTGGGCGGGTAGCAGCAAGCTGGCGGGAGTGAGTAGCAACTCATCCGGCTTGTCCATTGAGATAACCTGCCCGTCGTCAAGCGGTGAAGAGAGCCGTTCGAGACGGATATTGATAAGGCGACCACTGTCGGTACTGTTTTTGATCGTTTTGCTCAGCGTGCTGCTGTCACTGTTCATAAACGATGATATATCGCCGACATCCAGCGCCTGGGCCGGAGATATTCCGGAAAACAGCAGAGCGAGAGGCAGAAGGTGCTTTTTCATAACAATTCCGTCCAGGAATAAAGCAGGGGGGCTACCCCTGCTGGTACATCAGAGAGATTAACTGGTCCAGGTCGCGTCGAACTGTACGCTAACGTCGCCGCTCCAGATGCCTTCCGGTAGAGTGCTGTAATCTGTTACTGCGGTGGTACCATTGGTGGTACCGCTGATGATGGAGAAGGTGAAACCATCCTGT  >  NC_000913/309738‑310219

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: ‑