Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,619,539	(A)8→7	intergenic (+808/‑1452)	PP_5652 → / → PP_5654	hypothetical protein/phosphotriesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,619,532	0	A	.	100.0%	132.2 / NA	31	intergenic (+801/‑1459)	PP_5652/PP_5654	hypothetical protein/phosphotriesterase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  new base . (10/21);  total (10/21)

AAACAAGCAGCCTATACAGAGCCCAACGGCGGCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATATAAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTT  >  NC_002947/4619395‑4619675                                                                                                                                          |                                                                                                                                                aaaCAAGCAGCCTATACAGAGCCCAACGGCGGCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATATaaaaaaattcg                                                                                                                                       <  2:178852/148‑5 (MQ=255)          gcCTATACAGAGCCCAACGGCGGCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGc                                                                                                                             <  6:163537/148‑1 (MQ=255)          gcCTATACAGAGCCCAACGGCGACAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGc                                                                                                                             >  1:164569/1‑148 (MQ=255)                  agagCCCAACGGCGGCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCcaa                                                                                                                      >  7:83095/1‑147 (MQ=255)                  agagCCCAACGGCGGCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCcaa                                                                                                                      >  3:139420/1‑147 (MQ=255)                           cggcggCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGACACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACga                                                                                                            <  3:83438/148‑1 (MQ=255)                                       caacaTAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGaaa                                                                                                <  6:29619/148‑1 (MQ=255)                                        aacaTAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGaaa                                                                                                >  2:130443/1‑147 (MQ=255)                                        aacaTAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGaaa                                                                                                >  5:33032/1‑147 (MQ=255)                                            tAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAata                                                                                           <  2:164569/148‑1 (MQ=255)                                               gAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATAtt                                                                                         <  1:176328/147‑1 (MQ=255)                                               gAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGACACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATAtt                                                                                         <  1:176332/147‑1 (MQ=255)                                                  aTCTAGCGCCACAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTctct                                                                                     <  2:91431/148‑1 (MQ=255)                                                  aTCTAGCGCCAAATGTAAACTTTCCTATCATCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AGAAAAATTCGAATACCCTCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTctct                                                                                     <  4:105879/148‑1 (MQ=18)                                                  aTCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATtctc                                                                                      <  5:132835/147‑1 (MQ=255)                                                          ccAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTACTAGAGTATAT‑AAAAAAATTCGAATAGCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCaa                                                                             <  2:164556/148‑1 (MQ=255)                                                           cAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCATCAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAg                                                                            >  7:92675/1‑148 (MQ=255)                                                           cAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAg                                                                            <  6:90741/148‑1 (MQ=255)                                                                 tAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCg                                                                      >  5:26631/1‑148 (MQ=255)                                                                 tAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCg                                                                      <  6:5623/148‑1 (MQ=255)                                                                  aaaCTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGt                                                                     <  6:33582/148‑1 (MQ=255)                                                                            aTCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATg                                                           >  5:40524/1‑148 (MQ=255)                                                                                gCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTc                                                        <  5:162949/147‑1 (MQ=255)                                                                                      ttattGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGc                                                 <  6:33032/148‑1 (MQ=255)                                                                                       tattGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACAAAGCCCTCGTGAGCAAAATGGTCATTCGGCg                                                >  5:49954/1‑148 (MQ=255)                                                                                            aaTACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGATTATAT‑AAAAAAATTCGAATACCCGATCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCTGCGAttg                                            <  4:156068/147‑3 (MQ=17)                                                                                             aTACTGCATATTATCCAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATgga                                           <  2:76632/147‑1 (MQ=255)                                                                                              tACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATggag                                          <  1:34893/147‑1 (MQ=255)                                                                                                 tGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGaaa                                      >  5:27653/1‑148 (MQ=255)                                                                                                             aGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAgg                          <  4:155324/148‑1 (MQ=255)                                                                                                             aGTGAAAAAAGTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAg                           <  4:76897/147‑1 (MQ=255)                                                                                                               tGAAAAAACTGAATTATTAGAGTATAT‑AAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGcc                        <  8:159674/148‑1 (MQ=255)                                                                                                                                     tatat‑aAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGtt  <  5:90985/148‑1 (MQ=255)                                                                                                                                          |                                                                                                                                                AAACAAGCAGCCTATACAGAGCCCAACGGCGGCAGGGACAACATAAGAAATCTAGCGCCAAAGGTAAACTTTCCTATCAGCCACATTATTGAATACTGCATATTATCAAGTGAAAAAAGTGAATTATTAGAGTATATAAAAAAAATTCGAATACCCGCTCCCCAACAATACACGAGAGAAATTGAAAAAATATTCTCTGAAACCAAGCCCTCGTGAGCAAAATGGTCATTCGGCGATGGAGGAAACGCTGTACGAGGCCACGCGGATCCACGCGCCGAGTT  >  NC_002947/4619395‑4619675

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