Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	664,890	A→G	7.6%	S585G (AGC→GGC)	yfhM →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	664,890	0	A	G	7.6%	186.2 / 3.7	65	S585G (AGC→GGC)	yfhM	lipoprotein
Reads supporting (aligned to +/- strand):  ref base A (29/31);  new base G (3/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.77e-01

TTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGG  >  NC_002947/664749‑665034                                                                                                                                              |                                                                                                                                                 ttCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAg                                                                                                                                          <  2:134651/150‑1 (MQ=255)        cGGTAACCAAAACCGATCTAGCCCCGCACTGGGGACCGGACGAAGTCCCCCTTGACCAGGCCGCCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                   <  1:414707/150‑1 (MQ=255)           tCAACGAAAGCGATCGGAACCAGGACCTGGAACTGGACGAAGTCACCCTGGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTg                                                                                                                                >  2:951546/1‑150 (MQ=255)           cccccGAAACCGATCTGATCCGGGCCCCGGAACTGCACGCAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTg                                                                                                                                <  1:705865/147‑1 (MQ=255)            cACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGGCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTgg                                                                                                                               >  2:848759/1‑150 (MQ=255)               cGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGt                                                                                                                                 <  2:77722/145‑1 (MQ=255)               cGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGt                                                                                                                                 >  1:77722/1‑145 (MQ=255)               cGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTggccg                                                                                                                            >  1:1058638/1‑150 (MQ=255)                     cGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGAt                                                                                                                      >  2:1054015/1‑150 (MQ=255)                                  gACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTg                                                                                                         <  1:816480/150‑1 (MQ=255)                                   aCCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGGCGATATCGAAAGCCGTTGGGGCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTgg                                                                                                        >  1:470864/1‑150 (MQ=255)                                    ccTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGa                                                                                                       <  2:1008227/150‑1 (MQ=255)                                     cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                      >  2:803095/1‑150 (MQ=255)                                     cTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGag                                                                                                      >  1:946428/1‑150 (MQ=255)                                      tGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                   <  2:992941/120‑1 (MQ=255)                                      tGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCgg                                                                                                                                   >  1:992941/1‑120 (MQ=255)                                      tGGAACTGGACGAAGGCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagc                                                                                                     >  1:1102912/1‑150 (MQ=255)                                            tGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTa                                                                                               >  2:674076/1‑150 (MQ=255)                                              gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGt                                                                                                                                 <  2:96097/114‑1 (MQ=255)                                              gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGt                                                                                                                                 >  1:96097/1‑114 (MQ=255)                                              gACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATc                                                                                             >  1:509960/1‑150 (MQ=255)                                                 gAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTgg                                                                                          >  1:470509/1‑150 (MQ=255)                                                       aCCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCt                                                                                    <  1:803095/150‑1 (MQ=255)                                                         ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagca                                                                                                    <  2:381480/132‑1 (MQ=255)                                                         ccTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGagca                                                                                                    >  1:381480/1‑132 (MQ=255)                                                                 aGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGc                                                                          >  1:876987/1‑150 (MQ=255)                                                                    ccGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCgg                                                                       <  2:470509/150‑1 (MQ=255)                                                                      ggCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGc                                                                     <  2:1047587/150‑1 (MQ=255)                                                                         aaGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTg                                                                  >  2:1066106/1‑150 (MQ=255)                                                                            gCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGa                                                               <  1:1066106/150‑1 (MQ=255)                                                                            gCGGTGGGCGATATCGAAAGCCGTGGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCGGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGa                                                               <  1:154611/150‑1 (MQ=255)                                                                                   tCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGc                                                                                     <  1:913444/121‑1 (MQ=255)                                                                                   tCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGc                                                                                     >  2:913444/1‑121 (MQ=255)                                                                                         tCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAg                                                  <  2:509960/150‑1 (MQ=255)                                                                                          cGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAgg                                                 <  2:456814/150‑1 (MQ=255)                                                                                           gcaggCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGTGTGGCCGGCCGATCAGCCGGCGCATGGAGCAGCGTCGCTGGCCGGCTGAGCGCCTGCCGTGCCTGCGAGGCCTGTTCTAAggg                                                <  2:470864/146‑2 (MQ=255)                                                                                           gAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGt                                                <  1:951546/150‑1 (MQ=255)                                                                                            aaaGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTg                                               <  2:946428/150‑1 (MQ=255)                                                                                              aGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCTGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                             <  2:1058638/150‑1 (MQ=255)                                                                                                     gggCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCg                                      >  2:1071954/1‑150 (MQ=255)                                                                                                       gCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATc                                                                                                                     >  2:1029960/1‑69 (MQ=255)                                                                                                       gCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATc                                                                                                                     <  1:1029960/69‑1 (MQ=255)                                                                                                              cccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                              <  2:630893/133‑1 (MQ=255)                                                                                                              cccGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTga                                              >  1:630893/1‑133 (MQ=255)                                                                                                                    ctcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCTCGTGCTGGACCTGTTAATCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCg                       <  2:645035/150‑1 (MQ=255)                                                                                                                     tcCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGt                      >  1:726904/1‑150 (MQ=255)                                                                                                                      ccGCTGCAACTGAACGTGCAGGCCAGCCTGCAGGAGTCGGGTGGACGGTCGATCACCAGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTg                     <  1:848759/150‑1 (MQ=255)                                                                                                                          tGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                             >  2:242745/1‑122 (MQ=255)                                                                                                                          tGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                             <  1:242745/122‑1 (MQ=255)                                                                                                                          tGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGc                                 >  1:640846/1‑134 (MQ=255)                                                                                                                          tGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGc                                 <  2:640846/134‑1 (MQ=255)                                                                                                                             aaCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                        <  2:151279/124‑1 (MQ=255)                                                                                                                             aaCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                        >  1:151279/1‑124 (MQ=255)                                                                                                                              aCTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGa             <  1:1071954/150‑1 (MQ=255)                                                                                                                                tGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGt           >  1:92192/1‑150 (MQ=255)                                                                                                                                   tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                            >  1:142910/4‑98 (MQ=255)                                                                                                                                   tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGc                                                            <  2:142910/95‑1 (MQ=255)                                                                                                                                   tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                        >  2:60773/4‑118 (MQ=255)                                                                                                                                   tcctGCAGGCTGGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGAc                                        <  1:60773/115‑1 (MQ=255)                                                                                                                                     ctGCAGGCTGGCCTGCACGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGGGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGGGAGCGATAGCGACGGGCCCGTGGGATTCGAGGTCCtg      >  2:969463/2‑150 (MQ=255)                                                                                                                                       gCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                             <  2:132148/109‑1 (MQ=255)                                                                                                                                       gCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTgag                                             >  1:132148/1‑109 (MQ=255)                                                                                                                                         aGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCtggtgg  >  1:824416/1‑150 (MQ=255)                                                                                                                                             cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAGTTCGAGTTCCtgg     >  1:451612/1‑143 (MQ=255)                                                                                                                                             cAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAGTTCGAGTTCCtgg     <  2:451612/143‑1 (MQ=255)                                                                                                                                              |                                                                                                                                                 TTCGGCTCGGTCACCGAAACCGATCTGAACCAGGACCTGGAACTGGACGAAGTCACCCTTGACCAGGCCGGCAAGGCGGTGGTCGATATCGAAAGCCGTTGGGCCGAAGCCCGCTCTCCGCTGCAACTGACCGTGCAGGCCAGCCTGCAGGAGTCGGGTGGCCGGCCGATCACCCGGCGCCTGGAGCAGCCTATCTGGCCGGCTGAGCGCCTGCCGGGCCTGCGAGGCCTGTTCGAAGGTGAGGAGACCGATAGCGACGGGCCCGTGGAATTCGAGTTCCTGGTGG  >  NC_002947/664749‑665034

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