Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	1	1,070,512	G→A	intergenic (+408/‑413)	unknown → / → unknown	Type II secretory pathway, component PulD/Bacterial non‑heme ferritin (EC 1.16.3.2)

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	1	1,070,512	0	G	A	100.0%	79.0 / NA	31	intergenic (+408/‑413)	unknown/unknown	Type II secretory pathway, component PulD/Bacterial non‑heme ferritin (EC 1.16.3.2)
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/17);  total (14/17)

ATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACTAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACTTATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAAATAATCTCTTAACAATTTAGTTA  >  1/1070280‑1070733                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              aTACGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                                                                                                                              <  1:732398/250‑1 (MQ=255)    cGTGCTTGAATTTATGCTTACATTTTTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGa                                                                                                                                                                                                           >  1:684168/1‑250 (MQ=255)                                                                ttACAAGGCTTTATTTCTTCAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                                                                                                                              >  1:661888/1‑187 (MQ=255)                                                                               tcttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTag                                                                                                                               <  1:352855/251‑1 (MQ=255)                                                                                cttcAGTTTCCATCACGTTTAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTaga                                                                                                                              >  1:313791/1‑251 (MQ=255)                                                                                                  ttAACCAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTTATTAATCCGTAGCGTTTATTAACGAAGTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTtc                                                                                                            <  1:807161/251‑1 (MQ=255)                                                                                                                                           ttaatccgtagcgtttattaacgaagTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCAt                                                                   <  1:385787/226‑1 (MQ=255)                                                                                                                                             aatccgtagcgtttattaacgaagTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGCGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATag                                                                 >  1:425117/24‑251 (MQ=255)                                                                                                                                                 cgtagcgtttattaacgaagTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAggg                                                             <  1:898298/232‑1 (MQ=255)                                                                                                                                                                    gTCACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACAt                                          >  1:360393/1‑251 (MQ=255)                                                                                                                                                                      cACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGtt                                                                                          <  1:397864/201‑1 (MQ=255)                                                                                                                                                                      cACTAAATTCGCATAACATATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATc                                         >  1:296838/1‑250 (MQ=255)                                                                                                                                                                                     acatataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                                                                                                                              <  1:599112/67‑1 (MQ=255)                                                                                                                                                                                     acatataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATa                                                                                                                <  1:702760/161‑1 (MQ=255)                                                                                                                                                                                     acatataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTaa                         >  1:588803/4‑251 (MQ=255)                                                                                                                                                                                     acatataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTaa                         >  1:773672/4‑251 (MQ=255)                                                                                                                                                                                      catataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTaaaaa                                                                                                                                                                                                              >  1:141486/3‑69 (MQ=255)                                                                                                                                                                                      catataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTAc                                            >  1:70602/3‑231 (MQ=255)                                                                                                                                                                                      catataTTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCt                                       >  1:5661/3‑236 (MQ=255)                                                                                                                                                                                         atatTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAgg                                                                                            <  1:526854/180‑1 (MQ=255)                                                                                                                                                                                         atatTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTAc                                            >  1:360519/1‑228 (MQ=255)                                                                                                                                                                                                  taaataaaATAAATGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAAATAATCTCTTAAc            <  1:747323/251‑1 (MQ=255)                                                                                                                                                                                                        aaataaatGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTatgatg                                                                                                                                                      <  1:893031/107‑1 (MQ=255)                                                                                                                                                                                                        aaataaatGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAaat                       <  1:41805/234‑1 (MQ=255)                                                                                                                                                                                                        aaataaatGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAAATAATCTCTTAACAATTTa      >  1:449538/1‑251 (MQ=255)                                                                                                                                                                                                          ataaatGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAAATAATCTCTTAACAAt         >  1:648706/1‑246 (MQ=255)                                                                                                                                                                                                               tGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATAc                                                                                                               <  1:713089/139‑1 (MQ=255)                                                                                                                                                                                                               tGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATaa                                                                                     <  1:878604/165‑1 (MQ=255)                                                                                                                                                                                                               tGGTAAGTCACTTCTGGGAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATc                                         <  1:882549/209‑1 (MQ=255)                                                                                                                                                                                                                                gAATTTGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAAATAATCTCTTAACAATTTAGTTa  <  1:9922/231‑1 (MQ=255)                                                                                                                                                                                                                                   tttGAGATTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAgg                                                                                            <  1:233195/138‑1 (MQ=255)                                                                                                                                                                                                                                         |                                                                                                                                                                                                                              ATACGTGCTTGAATTTATGCTTACACTTCTAAATAAGTCTTGATAAATCAATATTTATTGATGTTACAAGGCTTTATTTCTTCAGTTTTCATAACGTTTAACTAGTAAGAGTGACCTGTGAATAACTCTGTGTTTTCATTAAATCGTAGTGTTCATTTACTTTGTCACTAATTTCGCATAACTTATATTATGTTAAATAAAATAAATGGTAAGTCACTTCTGGGAATTTGAGGTTTGCTCTCTTTAAAAATGACTGGGTAAAAAAATTTATTTCTTGAAGATTTGATTTTTTGAGCGTTTATGATGATGAGATATTCGTGAGTTTTTAGAGGTAGGTGAAAATACTTCTCTAGCTGAATCTAGGTTTATAAAAATCGACTGATATTCATAGAGGGTTTTTATCGTTATGTACATCCTTCCGTTTCCATTAAATAATCTCTTAACAATTTAGTTA  >  1/1070280‑1070733

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