Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,252,535	G→A	9.4%	intergenic (+22/+221)	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,535	0	G	A	9.4%	164.7 / 7.9	64	intergenic (+22/+221)	ureG/PP_2851	urease accessory protein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/26);  new base A (3/3);  total (35/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 = 2.09e-01

GGAGCGCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAC  >  NC_002947/3252387‑3252675                                                                                                                                                      |                                                                                                                                             ggAGCGCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGc                                                                                                                                             >  2:1035119/1‑150 (MQ=255)    gcgcGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATg                                                                                                                                          >  2:115637/1‑150 (MQ=255)      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                                                                                                   >  1:381524/1‑150 (MQ=255)                                           cTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTg                                                                                                   >  2:729574/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                                                                            <  2:7132/141‑1 (MQ=255)                                                                           aGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCcgc                                                                            >  1:7132/1‑141 (MQ=255)                                                                              ccGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGc                                                                <  2:867814/150‑1 (MQ=255)                                                                               cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACaga                                                                                                                                >  2:975624/1‑85 (MQ=255)                                                                               cGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGACATGCCCGCACaga                                                                                                                                <  1:975624/85‑1 (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                              >  2:638256/1‑137 (MQ=255)                                                                                                                             tGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGcc                              <  1:638256/137‑1 (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  <  1:115637/150‑1 (MQ=255)                                                                                                                                            tctGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAc  <  2:238051/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)                                                                                                                                                      |                                                                                                                                             GGAGCGCGACACCCAGCGCATGCGGCCGCAGCGGCCCTGGACCTTCAGCAACCTGAAGAAGGGCGAGGGGTTGCAGGCCGTGATCGATTTTATCGTCGAGCGCGGGATGCTGGGGGTGCGAGGCTGACAGATCGAGGACTCTGTG‑CGGGCATGTCCGCACAGAGACCAACGCGGTGGCTGGCACCGGCTCTGCCGGTGTTCGCGGGTTTACCCGCGAATACGGTAGCGGCCGCAACAATGAACGGCGGGTGGAAATTGGCCAGCAGGGCCGGCCTCTTCGCGGGTAAAC  >  NC_002947/3252387‑3252675

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