Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	167,295	T→C	20.8%	intergenic (+31/‑189)	mrcB → / → fhuA	fused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	167,295	0	T	C	20.8%	116.5 / 25.3	53	intergenic (+31/‑189)	mrcB/fhuA	fused glycosyl transferase and transpeptidase/ferrichrome outer membrane transporter
Reads supporting (aligned to +/- strand):  ref base T (13/29);  new base C (4/7);  total (17/36)
Fisher's exact test for biased strand distribution p-value = 7.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.14e-02

GCGTGGACTACGACGGCAACTTTGTTTGCAGCGGTGGCATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTG  >  NC_000913/167029‑167522                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    gCGTGGACAAGTACGGGAAATTGGGTTGCAGCGGGGGCATGCGTATCTTGCCGGGCTGGACCAGCGATACGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATcc                                                                                                                                                                                                                              <  4:373190/274‑1 (MQ=255)    tGGACTACGGCGGAAACTGTGTTTGCAGCGGTGGAAGGCGTATCGTGACGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGAAGCAGCAGCCGTAAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTa                                                                                                                                                                                                                          <  4:195899/275‑1 (MQ=255)     ggACTACGACGGCAACTTTGTTTGCAGCGGTGGCATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCACTCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGt                                                                                                                                                                                                                           >  2:91981/1‑273 (MQ=255)           cgacgGCAACTTTGTTTGCAGCGGTGGCATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAAc                                                                                                                                                                                                                   >  2:124878/1‑275 (MQ=255)             acgGCAACTTTGTTTGGAGCGGTGGCATGCGTATCTTGCCGGGCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGAAGTAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACAt                                                                                                                                                                                                                 <  2:475593/275‑1 (MQ=255)              cgGCAACTTTGTTTGCAGCGGTGGCATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACAt                                                                                                                                                                                                                 <  2:528808/274‑1 (MQ=255)               gacAACTTGGTATGCAGGGGTGACATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGGGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATc                                                                                                                                                                                                                <  2:210029/272‑1 (MQ=255)                      ttttctgCAGCGGTGGCAGGCGTATCTTGCCTGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGACAGGCAATCCGTTTGATCAGTCATCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGGCAGCGACGGTGTAGCCGGTTGAATCAAAGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACAGCCCATCTaa                                                                                                                                                                                                        <  1:97349/270‑1 (MQ=255)                      ttgtttgCAGCGGTGGCATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTaa                                                                                                                                                                                                        <  2:153462/275‑1 (MQ=255)                      ttgtttcCAGCGTTGGCATGCGTATCTTGCCGGTCTGGCCCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTATGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTa                                                                                                                                                                                                         <  1:250/274‑1 (MQ=255)                       tgtttgTAGCGGTTGCATGCGTATCTTGCCGGTCTGCACGAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTGTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGGAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTaa                                                                                                                                                                                                        <  2:89456/274‑1 (MQ=255)                                                   ccGGGCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgtt                                                                                                                                                                            <  2:273820/274‑1 (MQ=255)                                                          ggACCAGCGATCCGCAATCGCTGGGCCAGCAGAGCGGGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTttatt                                                                                                                                                                     <  4:91980/274‑1 (MQ=255)                                                                      cGCAATCGCTGTGCCAGCAGAGGGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGAAGCGACGGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGaa                                                                                                                                                        <  4:301763/275‑1 (MQ=255)                                                                      cGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATAAGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGaa                                                                                                                                                        <  2:247354/275‑1 (MQ=255)                                                                      cGCAATCGCTGGGCCAGCAGAGCGAGATGCAGCAGCAGCAGTCAGGCAATCCGGTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGAGCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGaa                                                                                                                                                        <  1:122866/275‑1 (MQ=255)                                                                         aaTCGCGTTGCCAGCCGGGCGCGATGCAGCAGCAGCCGTCAGGAAATCCGTTTGATCAGTCTTCTCAGCGGCAGGAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGGGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACg                                                                                                                                                      <  4:95053/274‑1 (MQ=255)                                                                           tCGCTGTGCCAGTAGAGCGGGATGCAGCACCATCCGTCAGGCCATCCTTTTGATCAGGCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTc                                                                                                                                                   <  4:367582/275‑1 (MQ=255)                                                                           tCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTc                                                                                                                                                   <  2:376131/275‑1 (MQ=255)                                                                              cTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGAAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAg                                                                                                                                                 >  2:532127/1‑274 (MQ=255)                                                                                  gCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTg                                                                                                                                             >  1:124577/1‑274 (MQ=255)                                                                                                             cggcAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGtt                                                                                                                  <  2:29540/271‑1 (MQ=255)                                                                                                                            ttGATCAGTCTTCTCAGCCGAAGCAACAGCCGACACAGCAACCTGCTCAGCAAGAGCAGAAAGACAAGGACGGTGTANCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATAttatcttat                                                                                                  <  3:124577/275‑1 (MQ=255)                                                                                                                                          cAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATcc                                                                                                                                                                                                                              >  6:1017295/1‑137 (MQ=255)                                                                                                                                          cAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATcc                                                                                                                                                                                                                              <  5:1017295/137‑1 (MQ=255)                                                                                                                                           aGCCGCAGCAACAGCCGCAACAGCAACCGGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCAt                                                                                   <  1:70233/275‑1 (MQ=255)                                                                                                                                                    aaCAGCCGCAACAGCAACCTGCTCAGCAAGAGAAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTa                                                                          <  4:532126/275‑1 (MQ=255)                                                                                                                                                              aCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCAt                                                                                                                                                                                                                                >  6:1029634/1‑115 (MQ=255)                                                                                                                                                              aCAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCAt                                                                                                                                                                                                                                <  5:1029634/115‑1 (MQ=255)                                                                                                                                                                                       aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                                                                                                     <  6:148136/85‑1 (MQ=255)                                                                                                                                                                                       aaaGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGtt                                                                                                                                                                                                                                     >  5:148136/1‑85 (MQ=255)                                                                                                                                                                                          gACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATc                                                                                                                                                                                                                               <  6:650999/88‑1 (MQ=255)                                                                                                                                                                                              gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATc                                                                                                                                                                                                                               <  7:243574/84‑1 (MQ=39)                                                                                                                                                                                              gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATc                                                                                                                                                                                                                               <  7:243296/84‑1 (MQ=39)                                                                                                                                                                                              gCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATc                                                                                                                                                                                                                               >  8:243296/1‑84 (MQ=39)                                                                                                                                                                                                  cGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATc                                                                                                                                                                                                                               >  6:87100/1‑80 (MQ=39)                                                                                                                                                                                                  cGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCAACTCGCCATc                                                                                                                                                                                                                               <  5:87100/80‑1 (MQ=39)                                                                                                                                                                                                               tGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATcc                                                                                                                                                                                                                              >  2:191636/1‑68 (MQ=255)                                                                                                                                                                                                                   tCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTa                                                                                                                                                                                                                          <  4:191635/68‑1 (MQ=255)                                                                                                                                                                                                                            tGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCAt                                                                                                                                                                                                            >  4:567584/1‑73 (MQ=255)                                                                                                                                                                                                                            tGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTg  >  2:343202/1‑275 (MQ=255)                                                                                                                                                                                                                               ttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGata                                                                                                                                                                                                    >  8:81460/1‑78 (MQ=255)                                                                                                                                                                                                                               ttGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACACCTCATCTAAGata                                                                                                                                                                                                    <  7:81460/78‑1 (MQ=255)                                                                                                                                                                                                                                tGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTaa                                                                                                                                                                                                        <  2:567585/73‑1 (MQ=255)                                                                                                                                                                                                                                      taATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTtta                                                                                                                                                                       >  2:447467/1‑100 (MQ=255)                                                                                                                                                                                                                                          tAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTaa                                                                                                                                                                   <  4:447466/100‑1 (MQ=255)                                                                                                                                                                                                                                                                ggATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                                                                                                             <  6:616017/84‑1 (MQ=255)                                                                                                                                                                                                                                                                ggATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                                                                                                             >  5:616017/1‑84 (MQ=255)                                                                                                                                                                                                                                                                ggATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGgttgt                                                                                                                                                                          <  6:1083710/71‑1 (MQ=38)                                                                                                                                                                                                                                                                 gATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                                                                                                             <  6:355497/83‑1 (MQ=255)                                                                                                                                                                                                                                                                 gATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTc                                                                                                                                                             >  5:355497/1‑83 (MQ=255)                                                                                                                                                                                                                                                                 gATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                                                                                                                 >  7:716393/1‑127 (MQ=255)                                                                                                                                                                                                                                                                 gATGGCAACTCGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGttt                                                                                                                 <  8:716393/127‑1 (MQ=255)                                                                                                                                                                                                                                                                        aGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAcc                                                                                                                            >  7:1062922/1‑109 (MQ=255)                                                                                                                                                                                                                                                                        aGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAAcc                                                                                                                            <  8:1062922/109‑1 (MQ=255)                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                    GCGTGGACTACGACGGCAACTTTGTTTGCAGCGGTGGCATGCGTATCTTGCCGGTCTGGACCAGCGATCCGCAATCGCTGTGCCAGCAGAGCGAGATGCAGCAGCAGCCGTCAGGCAATCCGTTTGATCAGTCTTCTCAGCCGCAGCAACAGCCGCAACAGCAACCTGCTCAGCAAGAGCAGAAAGACAGCGACGGTGTAGCCGGTTGGATCAAGGATATGTTTGGTAGTAATTAACATCTAAGCGTGAAATACCGGATGGCGAGTTGCCATCCGGTAAAATAACATCCCATCTAAGATATTAACCCTTTCTTTTCATCTGGTTGTTTATTAACCCTTCAGGAACGCTCAGATTGCGTACCGCTTGCGAACCCGCCAGCGTTTCGAATATTATCTTATCTTTATAATAATCATTCTCGTTTACGTTATCATTCACTTTACATCAGAGATATACCAATGGCGCGTTCCAAAACTGCTCAGCCAAAACACTCACTG  >  NC_000913/167029‑167522

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