Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,858,489	C→G	6.4%	P542P (CCC→CCG)	PP_1661 →	dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,858,489	0	C	G	6.4%	193.7 / 3.6	62	P542P (CCC→CCG)	PP_1661	dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/28);  new base G (2/2);  total (32/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 = 9.73e-01

GATTGGGCCTGCCTCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGACCACCCCAACCTGTTCC  >  NC_002947/1858341‑1858637                                                                                                                                                     |                                                                                                                                                     gATTGGGCCTGCCTCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCg                                                                                                                                                     <  1:376182/150‑1 (MQ=255)  aTTGGGCCTGCCTCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCgg                                                                                                                                                    >  2:509483/1‑150 (MQ=255)   ttGGGCCTGCCTCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGc                                                                                                                                                   <  2:1257625/150‑1 (MQ=255)          tgcctCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCaaa                                                                                                                                            >  2:1404384/1‑150 (MQ=255)          tgcctCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCaaa                                                                                                                                            >  2:849211/1‑150 (MQ=255)              tCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCGGGCTCgagc                                                                                                                                             <  2:222210/145‑4 (MQ=255)              tCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCGGGCTCgagc                                                                                                                                             >  1:222210/1‑142 (MQ=255)                  ccGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTc                                                                                                                                    <  1:23747/150‑1 (MQ=255)                     cAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAg                                                                                                                                 <  1:850357/150‑1 (MQ=255)                          cACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTacga                                                                                                                            >  1:174584/1‑150 (MQ=255)                           aCCTACGAACTGTCCGACTACACCAAGGAGGGTTGCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGGGTTCGAGTacgac                                                                                                                           >  2:1071436/1‑150 (MQ=255)                             cTACGAACTGTCCGGCTACACCAAGGAGGGTTTCAAGCAGGCACGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACgg                                                                                                                         >  1:473548/1‑150 (MQ=255)                                 gAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGCCGCTACCGAGCTGACCGACCTGAAACATAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAg                                                                                                                     >  1:805599/1‑150 (MQ=255)                                     tGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCGGGCTCgag                                                                                                                                              >  1:870215/1‑119 (MQ=255)                                     tGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCGGGCTCgag                                                                                                                                              <  2:870215/121‑3 (MQ=255)                                     tGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACt                                                                                                                 >  2:1225974/1‑150 (MQ=255)                                       tCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCATCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCGCGAGGCCGGGTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTAc                                                                                                               <  1:248649/150‑1 (MQ=255)                                       tCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTaa                                                                                                               >  1:733082/1‑149 (MQ=255)                                       tCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTAc                                                                                                               >  2:330829/1‑150 (MQ=255)                                       tCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTAc                                                                                                               <  2:535895/150‑1 (MQ=255)                                          gACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTAcgcg                                                                                                            >  1:110586/1‑150 (MQ=255)                                                acCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTAc                                                                                                      >  1:86131/1‑150 (MQ=255)                                                  cAAGGAGGGTTTCAAGCAGGGGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACTCTAAGCTCGAGCACGGCTCGCAAACGGGGTTCGAGTGCGACGGGCGGGACGACGCGTTCTACgg                                                                                                    >  2:435738/1‑150 (MQ=255)                                                       aggGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTg                                                                                               <  1:1225974/150‑1 (MQ=255)                                                       aggGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTg                                                                                               <  2:174584/150‑1 (MQ=255)                                                       aggGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTg                                                                                               >  1:687309/1‑150 (MQ=255)                                                             tCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCAcc                                                                                         <  2:5019/150‑1 (MQ=255)                                                               aaGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTg                                                                                       >  2:1158003/1‑150 (MQ=255)                                                                gccAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCGGAACTACGCGTTCTACGGCGCTGGCCACCTGa                                                                                      >  2:443680/3‑150 (MQ=255)                                                                 gCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGTTATAACTACGCTTTCTACGGCGCTGGCCACCTGAt                                                                                     <  1:1404384/150‑1 (MQ=255)                                                                 gCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGAt                                                                                     <  2:687309/150‑1 (MQ=255)                                                                       gcgcCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCAc                                                                               >  1:848293/1‑150 (MQ=255)                                                                                cGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCAt                                                                      <  2:781665/150‑1 (MQ=255)                                                                                      gTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTAcgcg                                                                                                            <  2:73808/106‑1 (MQ=255)                                                                                      gTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTAcgcg                                                                                                            >  1:73808/1‑106 (MQ=255)                                                                                      gTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTc                                                                >  1:1411238/1‑150 (MQ=255)                                                                                      gTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGACGATCACCTACCGCATGGGCTc                                                                >  1:784076/1‑150 (MQ=255)                                                                                        aCATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGCTGCGCCCCTACCGCACggcacagc                                                              <  2:733082/150‑7 (MQ=255)                                                                                        aCATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCgg                                                              >  1:102469/1‑150 (MQ=255)                                                                                          atcatcACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGAc                                                            <  1:635518/150‑1 (MQ=255)                                                                                            catcaCCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCg                                                          >  1:890637/1‑150 (MQ=255)                                                                                             atcaCCGAGCTGATGGGCGCTACCGGGAAGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGc                                                         >  1:749210/1‑150 (MQ=255)                                                                                                  cGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCaa                                                    <  1:330829/150‑1 (MQ=255)                                                                                                         aTGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCggtg                                             <  1:127120/150‑1 (MQ=255)                                                                                                         aTGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGCTCTACGGCGCTGGCCACCCGATGGGCACCTACCGCATGGGCTCGGCCCGCCCCAAGCCgctg                                             >  2:127120/1‑150 (MQ=255)                                                                                                         aTGGGCGCTACCGAGCTGACCGACCTGAACACTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCggtg                                             <  2:473548/150‑1 (MQ=255)                                                                                                           gggCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCggtggt                                           <  2:110586/150‑1 (MQ=255)                                                                                                                 tACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACaa                                     <  1:543903/150‑1 (MQ=255)                                                                                                                 tACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACaa                                     <  1:1158003/150‑1 (MQ=255)                                                                                                                              gaccTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACTTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGc                        <  1:292811/150‑1 (MQ=255)                                                                                                                                  tGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTggg                    >  2:1328423/1‑150 (MQ=255)                                                                                                                                  tGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTggg                    <  2:848293/150‑1 (MQ=255)                                                                                                                                       ccTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGGCGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGaccac               >  1:391208/1‑150 (MQ=255)                                                                                                                                       ccTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGaccac               >  2:919465/1‑150 (MQ=255)                                                                                                                                       ccTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGaccac               <  2:102469/150‑1 (MQ=255)                                                                                                                                       ccTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGaccac               >  1:588418/1‑150 (MQ=255)                                                                                                                                        cTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATggga                                                                                 >  1:878461/1‑82 (MQ=255)                                                                                                                                        cTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGc                                                                                 <  2:878461/83‑1 (MQ=255)                                                                                                                                        cTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTATCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGaccacc              <  2:86131/150‑1 (MQ=255)                                                                                                                                         tAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGAt                                                                                     >  2:899725/1‑78 (MQ=255)                                                                                                                                         tAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGAt                                                                                     <  1:899725/78‑1 (MQ=255)                                                                                                                                             cTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGACCACCCCAAc         <  2:805599/150‑1 (MQ=255)                                                                                                                                                    ccGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGACCACCCCAACCTGTTcc  <  1:2008/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GATTGGGCCTGCCTCGGCCGCAGATCACCTACGAACTGTCCGACTACACCAAGGAGGGTTTCAAGCAGGCGCGCCTGGCCGCTACGTACATCATCACCGAGCTGATGGGCGCTACCGAGCTGACCGACCTGAACCCTAAGCTCGAGCCCGGCTCGCAAACGGTGTTCGAGTACGACGGCCAGAACTACGCGTTCTACGGCGCTGGCCACCTGATGGGCACCTACCGCATGGGCTCGGACCGCACCAAGTCGGTGGTGGACAAGGACCAGCGCAGCTGGGACCACCCCAACCTGTTCC  >  NC_002947/1858341‑1858637

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