Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,043,935	C→T	N335N (AAC→AAT)	PP_4448 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,043,935	0	C	T	100.0%	125.1 / NA	38	N335N (AAC→AAT)	PP_4448	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (19/19);  total (19/19)

TTAACGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGTAT  >  NC_002947/5043793‑5044076                                                                                                                                               |                                                                                                                                              ttAACGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATggag                                                                                                                                           <  2:167172/147‑1 (MQ=255)     cGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGatat                                                                                                                                      <  1:5745/148‑1 (MQ=255)                  gACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGc                                                                                                                         >  3:57733/1‑148 (MQ=255)                   aCCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCt                                                                                                                        <  4:181465/148‑1 (MQ=255)                     ctactaAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTgg                                                                                                                      >  3:184032/1‑148 (MQ=255)                               aaggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAGGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCta                                                                                                            <  4:37100/148‑1 (MQ=255)                                aggTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGCCAACGAACGGATATTTCGATACTTTGAGAAAAACCTAGACTTCCCCCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTatat                                                                                                           >  7:146550/1‑148 (MQ=255)                                   tAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACt                                                                                                        >  6:54367/1‑148 (MQ=255)                                    aGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTc                                                                                                       <  6:14078/148‑1 (MQ=255)                                           agCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGa                                                                                                <  4:91064/148‑1 (MQ=255)                                           agCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGa                                                                                                <  1:101573/148‑1 (MQ=255)                                                  gACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTaaa                                                                                         >  8:88957/1‑148 (MQ=255)                                                           aGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGATAGCATATCAAATGGAGGATATGCTGATGGCTCGCTTGACAAATTCTATACTCGCTTCGAGTGTAAAGTTTccgat                                                                                <  8:3939/148‑4 (MQ=255)                                                           aGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCa                                                                                 <  2:186580/147‑1 (MQ=255)                                                           aGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCa                                                                                 <  4:188426/147‑1 (MQ=255)                                                             ctctAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGa                                                                              >  2:127358/1‑148 (MQ=255)                                                              tctAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGaa                                                                             >  8:40234/1‑148 (MQ=255)                                                                tAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACg                                                                           >  4:98686/1‑148 (MQ=255)                                                                tAGACGACAACGAACGGATATTTCGATACGTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATCTGCTGATGGCGCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACg                                                                           <  7:178887/148‑1 (MQ=255)                                                                   acgacAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACgg                                                                          >  2:164192/1‑146 (MQ=255)                                                                   acgacAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGt                                                                         <  1:79841/147‑1 (MQ=255)                                                                     gacAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGCTTCCCAGGAACGGTTAg                                                                      <  8:95057/148‑1 (MQ=255)                                                                      acAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAgg                                                                     <  3:98686/148‑1 (MQ=255)                                                                          cgaacgGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGcttc                                                                 <  7:40234/148‑1 (MQ=255)                                                                           gaacgGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGcttct                                                                >  5:109216/1‑148 (MQ=255)                                                                                        gATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCggg                                                   >  8:53912/1‑148 (MQ=255)                                                                                             tttGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGcc                                               >  5:176628/1‑147 (MQ=255)                                                                                                         cTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCacac                                  >  5:194345/1‑148 (MQ=255)                                                                                                         cTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCacac                                  <  4:121766/148‑1 (MQ=255)                                                                                                          tAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACt                                 >  5:153030/1‑148 (MQ=255)                                                                                                            gACTTCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACt                                 >  1:129362/1‑146 (MQ=255)                                                                                                               ttCCACCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACATCACACTCgaga                            >  8:114081/1‑148 (MQ=255)                                                                                                                   aCCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAacac                        >  8:117024/1‑148 (MQ=255)                                                                                                                   aCCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAaca                         <  1:72342/147‑1 (MQ=255)                                                                                                                   aCCCAACCAAGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGAGCCAAAGACACCACACTCGAGAacac                        >  3:141854/1‑148 (MQ=255)                                                                                                                        accaaGCGAGAAAGCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTcc                   <  7:88957/148‑1 (MQ=255)                                                                                                                                     gCATATCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGACCCCTCTCCGAAACCTCATAgg      >  6:42673/1‑148 (MQ=255)                                                                                                                                         atCAAATGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGtat  <  6:153030/148‑1 (MQ=255)                                                                                                                                               |                                                                                                                                              TTAACGATGGTACATACGACCTACTAAAGGAAGGTAGAACAGAGCCTACGACCGGTGTAGCTCTAGACGACAACGAACGGATATTTCGATACTTTGAGAAAAAACTAGACTTCCACCCAACCAAGCGAGAAAGCATATCAAACGGAGGATATGCTGATGGCTCGCTGGACAAATTCTATACTCGCTTCGAGTCTAAAGTTTCCCAGGAACGGTTAGGCTTCTTGTTCGGCGCGGGCGCCAAAGACACCACACTCGAGAACACTCTCCGAAACCTCATAGGGTAT  >  NC_002947/5043793‑5044076

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