Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,252,540	T→C	9.8%	intergenic (+27/+216)	ureG → / ← PP_2851	urease accessory protein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,252,540	0	T	C	9.8%	160.6 / 7.8	62	intergenic (+27/+216)	ureG/PP_2851	urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/26);  new base C (3/3);  total (33/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

GCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGG  >  NC_002947/3252392‑3252689                                                                                                                                                      |                                                                                                                                                      gcgACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTc                                                                                                                                                      <  1:1035119/150‑1 (MQ=255)   gACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCg                                                                                                                                                    <  2:362058/150‑1 (MQ=255)     cacCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGca                                                                                                                                                  <  1:284767/150‑1 (MQ=255)          aGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACagag                                                                                                                                             <  2:221632/150‑1 (MQ=255)            cgcATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGAc                                                                                                                                           <  1:305608/150‑1 (MQ=255)                     ccGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGt                                                                                                                                  <  2:934963/150‑1 (MQ=255)                            ggCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGtggctggc                                                                                                                           >  2:936730/1‑150 (MQ=255)                              cccTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCAc                                                                                                                         >  1:619398/1‑150 (MQ=255)                                 tGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCgg                                                                                                                      >  1:948880/1‑150 (MQ=255)                                 tGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCgg                                                                                                                      <  2:948880/150‑1 (MQ=255)                                    aCCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGctc                                                                                                                   <  1:1246656/150‑1 (MQ=255)                                      cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                                 >  2:729574/1‑150 (MQ=255)                                      cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                                 >  1:381524/1‑150 (MQ=255)                                      cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGGGGCTGGCACCGGCTCTg                                                                                                                 >  1:368691/1‑150 (MQ=255)                                       ttCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGc                                                                                                                >  2:497165/1‑150 (MQ=255)                                       tccAGCAAACTGAGTAAGGTGGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGc                                                                                                                <  2:357987/148‑1 (MQ=255)                                             aaCCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGgtgt                                                                                                          >  1:659822/1‑150 (MQ=255)                                                cTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcg                                                                                                       >  1:209530/1‑150 (MQ=255)                                                 tGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTAATAGTAGAGCGCGGAACACTGAGAGCGCAAGAGTGACAGATAG‑GGAACGAGTGTCGAGCATTGTGTCGCTGAGGGCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTcgc                                                                                                      >  2:145233/1‑150 (MQ=255)                                                  gaagaaGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAAc                                                                                                                                       <  1:993469/116‑1 (MQ=255)                                                  gaagaaGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAAc                                                                                                                                       >  2:993469/1‑116 (MQ=255)                                                   aagaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCgg                                                                                                    <  2:224889/150‑1 (MQ=255)                                                   aagaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCgg                                                                                                    <  2:158647/150‑1 (MQ=255)                                                    agaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCg                                                                                                              <  2:509994/139‑1 (MQ=255)                                                    agaagGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCg                                                                                                              >  1:509994/1‑139 (MQ=255)                                                        gggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACaga                                                                                                                                              >  1:844148/1‑103 (MQ=255)                                                        gggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACaga                                                                                                                                              <  2:844148/103‑1 (MQ=255)                                                         ggCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTAc                                                                                              >  1:867814/1‑150 (MQ=255)                                                               gggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGAGCCGCGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGGGa                                                                                        >  1:977953/1‑150 (MQ=255)                                                                ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGaa                                                                                       >  1:1001822/1‑150 (MQ=255)                                                                ggTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTAAACCCGCGaa                                                                                       <  1:729574/150‑1 (MQ=255)                                                                      aGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                          >  1:7132/1‑141 (MQ=255)                                                                      aGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                                          <  2:7132/141‑1 (MQ=255)                                                                         ccGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                              <  2:867814/150‑1 (MQ=255)                                                                          cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACaga                                                                                                                                              <  1:975624/85‑1 (MQ=255)                                                                          cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACaga                                                                                                                                              >  2:975624/1‑85 (MQ=255)                                                                            tGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGc                                                                           >  1:1238764/1‑150 (MQ=255)                                                                             gATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGAACAGAGACCAACGCGGTGGCTGGAACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGcc                                                                          >  1:1095193/1‑150 (MQ=255)                                                                                 gATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑TGGGCATGTCCGCACAGAGACAAACGCGGTGGCTGGCACCGGCTCTGCCGGTGATCGCGGGTTGACCCGCGAATACGGTAGCGGCGGcaa                                                                      >  2:163165/1‑150 (MQ=255)                                                                                     ttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGTTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                                  <  1:92647/150‑1 (MQ=255)                                                                                     ttATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAAt                                                                  <  2:381524/150‑1 (MQ=255)                                                                                      tATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACACTg                                                                 >  2:306103/1‑150 (MQ=255)                                                                                       aTCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGa                                                                >  1:271177/1‑150 (MQ=255)                                                                                               gcgcGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCggg                                                        <  1:936730/150‑1 (MQ=255)                                                                                                 gcgGGATGCTGGGGGTGCTAGGCTGACAGATCCAGGACTCTGTG‑CGAGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTCTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTg                                                      <  1:163165/150‑1 (MQ=255)                                                                                                  cgGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCAGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTgg                                                     >  1:919897/1‑150 (MQ=255)                                                                                                   gggATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGa                                                    <  2:67952/150‑1 (MQ=255)                                                                                                         cTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTgg                                              >  1:14813/1‑150 (MQ=255)                                                                                                          tGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGc                                             <  2:619398/150‑1 (MQ=255)                                                                                                            ggggTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCca                                           >  1:238051/1‑150 (MQ=255)                                                                                                                       cTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGggccggc                                >  2:760455/1‑150 (MQ=255)                                                                                                                        tGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcc                                            <  1:638256/137‑1 (MQ=255)                                                                                                                        tGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcc                                            >  2:638256/1‑137 (MQ=255)                                                                                                                          aCAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATg                                                                 <  2:563965/114‑1 (MQ=255)                                                                                                                          aCAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATg                                                                 >  1:563965/1‑114 (MQ=255)                                                                                                                               tCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAACAGGGCCGGCCTCTTcgc                        >  2:256652/1‑150 (MQ=255)                                                                                                                                      ctctGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTaaa                 >  1:466433/1‑150 (MQ=255)                                                                                                                                       tctGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                <  2:238051/150‑1 (MQ=255)                                                                                                                                       tctGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                <  1:115637/150‑1 (MQ=255)                                                                                                                                       tctGTG‑CGGACATGCCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc                >  1:740039/1‑150 (MQ=255)                                                                                                                                         tgtg‑CGGGCATGTCCGCACAGAGCCCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaa                                                                                       <  1:996962/77‑1 (MQ=18)                                                                                                                                         tgtg‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTATTCGCGGGTAAACCCGCGaa                                                                                       >  2:996962/1‑77 (MQ=25)                                                                                                                                                      tCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACgggg  >  2:926913/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                      GCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACGGGG  >  NC_002947/3252392‑3252689

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