Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,529,614	C→G	6.6%	G477G (GGC→GGG)	dnaEB →	error‑prone DNA polymerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,529,614	0	C	G	6.6%	195.3 / 4.2	61	G477G (GGC→GGG)	dnaEB	error‑prone DNA polymerase
Reads supporting (aligned to +/- strand):  ref base C (29/28);  new base G (2/2);  total (31/30)
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.28e-01

ACCTGCCGATCAGGTGGATGCCCTGGCCAAATGTTGTGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTGGCCTGCT  >  NC_002947/3529467‑3529760                                                                                                                                                     |                                                                                                                                                   aCCTGCCGATCAGGTGGATGCCCTGGCCAAATGTTGTGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCt                                                                                                                                                   >  6:49744/1‑149 (MQ=255)       cGATCAGGTGGATGCCCTGGCCAAATGTTGTGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCTAACTGATCGGCTTccccc                                                                                                                                             <  4:288324/149‑1 (MQ=255)             ggTGGATGCCCTGGCCAAATGTAAAGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCAcc                                                                                                                                       >  7:204668/1‑149 (MQ=255)                      ccTGGCCAAATGTTGTGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAg                                                                                                                               >  6:132210/1‑148 (MQ=255)                                 gttgtGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGt                                                                                                                   >  2:309498/1‑149 (MQ=255)                                    gtgGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGttt                                                                                                                 <  8:173419/148‑1 (MQ=255)                                    gtgGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGttt                                                                                                                 >  7:173419/1‑148 (MQ=255)                                     tgGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTg                                                                                                               >  6:316333/1‑149 (MQ=255)                                      ggTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGt                                                                                                              >  6:294982/1‑149 (MQ=255)                                       ggcGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGt                                                                                                              <  8:2622/146‑1 (MQ=255)                                       gTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTc                                                                                                             <  5:132210/149‑1 (MQ=255)                                          gCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGtttt                                                                                                                >  3:290022/1‑143 (MQ=255)                                          gCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGtttt                                                                                                                <  4:290022/143‑1 (MQ=255)                                                gCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGGCGGCGAACTGATCGGCTTCCCCCGGCACCt                                                                                                                                      <  1:26416/115‑1 (MQ=255)                                                gCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCt                                                                                                                                      >  2:26416/1‑115 (MQ=255)                                                  gATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGc                                                                                                   >  6:309393/1‑148 (MQ=255)                                                          aCCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTg                                                                                           >  1:493/1‑148 (MQ=255)                                                            ccGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGAc                                                                                        <  5:49744/149‑1 (MQ=255)                                                            ccGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGAc                                                                                        >  8:241508/1‑149 (MQ=255)                                                                   ccAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTgg                                                                                 >  2:3620/1‑149 (MQ=255)                                                                     aGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTg                                                                               >  7:31936/1‑149 (MQ=255)                                                                      gcgcCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGc                                                                              <  5:294982/149‑1 (MQ=255)                                                                        gcCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTGCCTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGcc                                                                             >  1:246947/1‑149 (MQ=255)                                                                         ccTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCgg                                                                           <  8:204668/149‑1 (MQ=255)                                                                          cTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGt                                                                          >  4:154780/1‑149 (MQ=255)                                                                            ggCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTg                                                                         >  2:295947/1‑148 (MQ=255)                                                                            ggCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTg                                                                         >  5:306085/1‑148 (MQ=255)                                                                             gCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTg                                                                         <  4:195815/147‑1 (MQ=255)                                                                               cGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGGTTCCCCCGGCACCTGTCACAGc                                                                                                                              >  3:207713/1‑92 (MQ=255)                                                                               cGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGGTTCCCCCGGCACCTGTCACAGc                                                                                                                              <  4:207713/92‑1 (MQ=255)                                                                               cGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGGTTCCCCCGGCACCTGTCACAGc                                                                                                                              <  4:203631/92‑1 (MQ=255)                                                                               cGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGGTTCCCCCGGCACCTGTCACAGc                                                                                                                              >  3:203631/1‑92 (MQ=255)                                                                                     cGGTTTCGAAGCGGGCAGCCCTTCGCTGCGNCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCcg                                                               >  1:29257/1‑149 (MQ=255)                                                                                      ggTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCcgc                                                              <  5:316333/149‑1 (MQ=255)                                                                                          tCGAAGCGGGCAGCCCTTCGCTGGGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGGCATCTCCCAGCTACCGCTGGACCAGCTGGTGCCGGTGGAGATTGCCGCGAt                                                           <  5:164086/148‑1 (MQ=255)                                                                                             aaGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCg                                                       >  5:204517/1‑149 (MQ=255)                                                                                                      gCCCTTCGCTGCGGCGGGGG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACg                                              <  4:169249/149‑1 (MQ=255)                                                                                                          ttCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAg                                                                                                                               >  8:70032/1‑64 (MQ=255)                                                                                                          ttCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAg                                                                                                                               <  7:70032/64‑1 (MQ=255)                                                                                                          ttCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGa                                          >  3:201312/1‑149 (MQ=255)                                                                                                             gCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATc                                        <  2:208011/148‑1 (MQ=255)                                                                                                                gcggcgGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGt                                    >  1:257919/1‑149 (MQ=255)                                                                                                                 cggcggGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGt                                    <  2:246947/148‑1 (MQ=255)                                                                                                                 cggcggGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTg                                   <  1:295947/149‑1 (MQ=255)                                                                                                                    cggGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGa                                <  6:204517/149‑1 (MQ=255)                                                                                                                        tg‑ctggtgctggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAg                            <  2:29257/149‑1 (MQ=255)                                                                                                                           ctggtgctggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCATGGTGATCCAGTGGGACAAGGa                          <  8:31936/149‑1 (MQ=255)                                                                                                                            tggtgctggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGAt                         >  3:87301/1‑149 (MQ=255)                                                                                                                              gtgctggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGAt                         <  5:309393/147‑1 (MQ=255)                                                                                                                              gtgctggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGa                       >  1:57308/1‑149 (MQ=255)                                                                                                                                  tggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTg                   >  6:62699/1‑149 (MQ=255)                                                                                                                                  tggCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTg                   <  2:57308/149‑1 (MQ=255)                                                                                                                                    gCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTg                   <  5:251663/147‑1 (MQ=255)                                                                                                                                    gCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGa                 <  7:241508/149‑1 (MQ=255)                                                                                                                                             cTGATCGGCTTCCACCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTgg        >  4:37187/1‑149 (MQ=255)                                                                                                                                              tGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTGGc       >  3:237529/1‑149 (MQ=255)                                                                                                                                                aTCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTGGc       <  1:309498/147‑1 (MQ=255)                                                                                                                                                   ggCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTg         <  5:122416/142‑1 (MQ=255)                                                                                                                                                   ggCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTg         >  6:122416/1‑142 (MQ=255)                                                                                                                                                   ggCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTGGCCTGc   <  5:62699/148‑1 (MQ=255)                                                                                                                                                   ggCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTGGCCTGCt  <  1:189367/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                   ACCTGCCGATCAGGTGGATGCCCTGGCCAAATGTTGTGGTCGCTGGAGCGATCGCATACCCGATGACCAGCGCCTGGCCGAGGCCGGTTTCGAAGCGGGCAGCCCTTCGCTGCGGCGGGTG‑CTGGTGCTGGCCGGCGAACTGATCGGCTTCCCCCGGCACCTGTCACAGCACCCGGGCGGTTTTGTCATCTCCCAGCAACCGCTGGACCAGCTGGTGCCGGTGGAGAATGCCGCGATGCCGGAGCGCACGGTGATCCAGTGGGACAAGGATGACCTGGACATGGTTGGCCTGCT  >  NC_002947/3529467‑3529760

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