Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,246,242	A→G	8.6%	intergenic (+31/+60)	ibpA → / ← PP_1983	small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,246,242	0	A	G	8.6%	158.6 / 5.3	58	intergenic (+31/+60)	ibpA/PP_1983	small heat shock protein IbpA/PAS/PAC/GAF sensor‑containing diguanylate cyclase/phosphodiesterase
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base G (2/3);  total (30/28)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

AGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGG  >  NC_002947/2246099‑2246389                                                                                                                                                |                                                                                                                                                        aGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGt                                                                                                                                                    <  4:237933/149‑1 (MQ=255)    ccgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                                 >  5:85436/1‑149 (MQ=255)    ccgccgGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                                 <  1:96799/149‑1 (MQ=255)         ggCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑tttt                                                                                                                                          <  6:85436/149‑1 (MQ=255)               gCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCAAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑tttttttcgt                                                                                                                                    <  6:234585/149‑6 (MQ=255)                 cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                        >  5:244250/1‑143 (MQ=255)                 cAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTc                                                                                                                                        <  6:244250/143‑1 (MQ=255)                              aGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTccc                                                                                                                     <  6:109055/149‑1 (MQ=255)                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                      >  6:294426/1‑111 (MQ=255)                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCagag                                                                                                                                                      <  5:294426/111‑1 (MQ=255)                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg‑g                                                                                                             >  1:71406/1‑149 (MQ=255)                                     aCCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGATGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg‑g                                                                                                             >  7:266436/1‑149 (MQ=255)                                       ctgctgCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑ttt                                                                                                                                           >  8:12097/1‑118 (MQ=255)                                       ctgctgCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCTGGAGGTGCCC‑‑ttt                                                                                                                                           <  7:12097/118‑1 (MQ=255)                                       ctgctgCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTc                                                                                                                       >  2:199064/1‑138 (MQ=255)                                       ctgctgCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTc                                                                                                                       <  1:199064/138‑1 (MQ=255)                                              gcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCcgcg                                                                                                                  >  5:350576/1‑136 (MQ=255)                                              gcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCcgcg                                                                                                                  <  6:350576/136‑1 (MQ=255)                                              gcATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGcc                                                                                                   <  6:84355/149‑1 (MQ=255)                                               cATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                                >  7:173284/1‑137 (MQ=255)                                               cATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTa                                                                                                                <  8:173284/137‑1 (MQ=255)                                                    tGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATg                                                                                             >  1:177803/1‑149 (MQ=255)                                                       cTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGa                                                                                          >  6:133430/1‑149 (MQ=255)                                                           agaagCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCa                                                                                       <  2:177803/148‑1 (MQ=255)                                                               gccaaggccaatCGAATTCCTATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAAt                                                                                  <  6:356635/149‑1 (MQ=255)                                                                ccaaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATc                                                                                 <  2:71406/149‑1 (MQ=255)                                                                 caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                               >  4:227718/1‑120 (MQ=255)                                                                 caaggccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAg                                                                                                               <  3:227718/120‑1 (MQ=255)                                                                     gccaagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCg                                                                            >  1:338114/1‑149 (MQ=255)                                                                        aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                         >  4:180017/1‑149 (MQ=255)                                                                        aagCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGgcg                                                                         >  2:224237/1‑149 (MQ=255)                                                                            gcATTCCTATAAATTGAGAAAAACCAGCGCTTAACTTATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCTTCATGCTATCTGTAATCTGCCGGCGCCTc                                                                     <  8:266436/149‑1 (MQ=255)                                                                              aTTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCa                                                                                       <  6:38878/129‑1 (MQ=255)                                                                              aTTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCa                                                                                       >  5:38878/1‑129 (MQ=255)                                                                               ttCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAg                                                                  <  2:338114/149‑1 (MQ=255)                                                                                   gATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCaa                                                               >  2:30799/1‑148 (MQ=255)                                                                                        aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                                 >  4:158618/1‑65 (MQ=255)                                                                                        aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGcc                                                                                                                                                 <  3:158618/65‑1 (MQ=255)                                                                                        aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTtctc                                                                                                                               <  1:18894/81‑1 (MQ=255)                                                                                        aTGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTtctc                                                                                                                               >  2:18894/1‑81 (MQ=255)                                                                                         tGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTc                                                                                                >  6:252887/1‑109 (MQ=255)                                                                                         tGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTc                                                                                                <  5:252887/109‑1 (MQ=255)                                                                                           gCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCg                                                      <  5:133430/149‑1 (MQ=255)                                                                                              aaaaaCCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGaaa                                                   <  3:180017/149‑1 (MQ=255)                                                                                                    cAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCATCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCACCATCTGCCGAAACTCcgg                                             >  3:204566/1‑149 (MQ=255)                                                                                                           gAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCcggtaacgg                                       >  4:9860/1‑148 (MQ=255)                                                                                                                        tGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAGAGAGTCGGGGCCGTGATTGGGTGAATAATGTTCCCGCGGCTGCCGCAACATGTGCCGGAAAGCCGGTAACGGTAATGGCCgggt                           >  1:286099/1‑147 (MQ=255)                                                                                                                          cAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCcgtg                        >  8:201186/1‑148 (MQ=255)                                                                                                                           aaCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCtgtg                        >  7:168863/1‑147 (MQ=255)                                                                                                                            aCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTc                                                                                                <  3:17755/74‑1 (MQ=255)                                                                                                                            aCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTc                                                                                                >  4:17755/1‑74 (MQ=255)                                                                                                                            aCGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAg                     <  3:9860/149‑1 (MQ=255)                                                                                                                             cGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAg                     >  3:188096/1‑148 (MQ=255)                                                                                                                               aaaaaaaGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGAt                   >  2:200994/1‑148 (MQ=255)                                                                                                                                    aaGGGCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAAcc               >  4:5096/1‑149 (MQ=255)                                                                                                                                      gggCACCTCTGGAGGTGCCCTTTTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCt             <  3:5096/149‑1 (MQ=255)                                                                                                                                      gggCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCCGTGCAGATAACCCTgg           <  7:201186/149‑1 (MQ=255)                                                                                                                                            cTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTa         >  3:312134/1‑145 (MQ=255)                                                                                                                                               cAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATgg  >  7:227042/1‑149 (MQ=255)                                                                                                                                               cAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATgg  <  8:227042/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                        AGGCCGCCGGCCTGGCCAATGGCCTGCTCAGCATCGACCTGCTGCGCATCGTGCCTGAAGAAGCCAAGGCCAAGCGCATTCCGATCAATGGCGAAAAACCAGCGCTGAACTGATCGAGCTGCAACGAAAAAAAGGGCACCTCCAGAGGTGCCC‑‑TTTTTCGTGCTTCTCACGCAGTCCCGCGTAG‑GTG‑CGCGGCCGTCATGCGATCAGTAATCTGCCGGCGCCTCCAGCAACATCTGCCGAAACTCCGGTAACGGTAATGGCCGGCTGTGCAGATAACCCTGGTACAGATGG  >  NC_002947/2246099‑2246389

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