Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,605,915	C→A	5.5%	G283C (GGC→TGC)	PP_3181 ←	site‑specific recombinase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,605,915	0	C	A	5.5%	236.9 / 3.6	73	G283C (GGC→TGC)	PP_3181	site‑specific recombinase
Reads supporting (aligned to +/- strand):  ref base C (31/38);  new base A (2/2);  total (33/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

TATGCCATCCGCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCCCACCAGGCTGTTG  >  NC_002947/3605768‑3606060                                                                                                                                                    |                                                                                                                                                     tATGCCATCCGCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGcc                                                                                                                                                      >  5:190076/1‑148 (MQ=255) tATGCCATCCGCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGcc                                                                                                                                                      >  3:84612/1‑148 (MQ=255)     ccATCCGCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCcaca                                                                                                                                                  <  6:168275/148‑1 (MQ=255)           gCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCtca                                                                                                                                           >  5:103971/1‑149 (MQ=255)           gCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCtc                                                                                                                                            >  4:254919/1‑148 (MQ=255)             ccGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCtcatc                                                                                                                                         <  6:103971/149‑1 (MQ=255)              cGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCtcat                                                                                                                                          >  1:243103/1‑147 (MQ=255)                  aaCGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGa                                                                                                                                    <  3:254919/149‑1 (MQ=255)                     gggCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACtt                                                                                                                                 >  4:29856/1‑149 (MQ=255)                       gcggcgTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCa                                                                                                                               <  3:273212/149‑1 (MQ=255)                        cggcgTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAg                                                                                                                              <  3:29856/149‑1 (MQ=255)                          gcgTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGtt                                                                                                                            <  4:261462/149‑1 (MQ=255)                                 cAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAAcc                                                                                                                     <  4:84612/149‑1 (MQ=255)                                 cAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAAc                                                                                                                      >  3:221015/1‑148 (MQ=255)                                     cacaGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTTGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTg                                                                                                                 <  3:128747/149‑1 (MQ=255)                                     cacaGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAg                                                                                                                                               >  6:38205/1‑119 (MQ=255)                                     cacaGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAg                                                                                                                                               <  5:38205/119‑1 (MQ=255)                                     cacaGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTg                                                                                                                 <  4:234053/149‑1 (MQ=255)                                      acaGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTgg                                                                                                                <  1:266181/149‑1 (MQ=255)                                           ttCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTg                                                                                                           >  7:73004/1‑149 (MQ=255)                                             cAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTgg                                                                                                          <  4:221015/148‑1 (MQ=255)                                             cAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTgg                                                                                                          <  6:190076/148‑1 (MQ=255)                                             cAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTg                                                                                                           <  3:107525/147‑1 (MQ=255)                                              aaGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGt                                                                                                                             <  3:285454/128‑1 (MQ=255)                                              aaGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGt                                                                                                                             >  4:285454/1‑128 (MQ=255)                                               aGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGAc                                                                                                                                   <  8:97851/121‑1 (MQ=255)                                               aGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGAc                                                                                                                                   >  7:97851/1‑121 (MQ=255)                                                ggAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCtc                                                                                                                                            <  4:118386/111‑1 (MQ=255)                                                         ggTTGCAACTGACTGGCCGGTAGCACGGCCAGATCGCCAAATCGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGaa                                                                                           <  4:261883/149‑1 (MQ=255)                                                                ccTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGc                                                                                    <  5:314328/149‑1 (MQ=255)                                                                 cTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACg                                                                                     >  8:275458/1‑147 (MQ=255)                                                                 cTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGcc                                                                                   >  4:207776/1‑149 (MQ=255)                                                                 cTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCAACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGcc                                                                                   >  3:212350/1‑149 (MQ=255)                                                                  tGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCACCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGaa                                                                                           <  6:177701/140‑1 (MQ=255)                                                                  tGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCACCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGaa                                                                                           >  5:177701/1‑140 (MQ=255)                                                                        gCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGa                                                                            <  7:275458/149‑1 (MQ=255)                                                                         ccGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGa                                                                            <  3:289962/148‑1 (MQ=255)                                                                           ggTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGa                                                                            >  5:282247/1‑146 (MQ=255)                                                                            gTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTgg                                                                        >  2:137539/1‑149 (MQ=255)                                                                               gCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGa                                                                            <  1:166019/142‑1 (MQ=255)                                                                               gCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGa                                                                            >  2:166019/1‑142 (MQ=255)                                                                                   ggCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGta                                                                  >  4:192097/1‑148 (MQ=255)                                                                                      cAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGtata                                                                >  5:270057/1‑147 (MQ=255)                                                                                           cGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGcgcc                                                         >  2:93581/1‑149 (MQ=255)                                                                                             ccAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGcgccgc                                                       <  6:74758/149‑1 (MQ=255)                                                                                             ccAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGcgccgc                                                       >  1:262905/1‑149 (MQ=255)                                                                                             ccAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCACACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGcc                                                                                                        <  3:46701/102‑1 (MQ=255)                                                                                             ccAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCACACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGcc                                                                                                        <  1:67392/102‑1 (MQ=255)                                                                                             ccAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCACACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGcc                                                                                                        >  4:46701/1‑102 (MQ=255)                                                                                             ccAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCACACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGcc                                                                                                        >  2:67392/1‑102 (MQ=255)                                                                                                      tGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCAcccc                                              >  5:291240/1‑149 (MQ=255)                                                                                                          ccGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGc                                                                       <  4:110910/120‑1 (MQ=255)                                                                                                          ccGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGc                                                                       >  3:110910/1‑120 (MQ=255)                                                                                                           cGAGTCGATGTAGCGCAGCGCAGACTTGGCGTCTTTCCCGCCCACATAGGTCATCAGCGACTTCAGTTAGCAACCATTGGCCGTGGCACCG‑‑CTGGCGAAGCCACGCGGTAGCGCGTGGCCGGTATACAGCGCCGCCGGCACCCCGgcgc                                         >  8:176940/1‑149 (MQ=255)                                                                                                            gAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGtt                                                                                                                            <  7:104833/67‑1 (MQ=255)                                                                                                            gAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGtt                                                                                                                            >  8:104833/1‑67 (MQ=255)                                                                                                               tCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGTCGC‑GTTc                                    <  2:243103/149‑1 (MQ=255)                                                                                                                 gATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCa                                  <  2:262905/149‑1 (MQ=255)                                                                                                                 gATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCa                                  <  3:207776/149‑1 (MQ=255)                                                                                                                 gATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCa                                  <  6:291240/149‑1 (MQ=255)                                                                                                                  aTGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAg                                 >  5:9916/1‑149 (MQ=255)                                                                                                                   tGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAGTCGTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCa                                  <  1:253007/149‑1 (MQ=255)                                                                                                                   tGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGc                                <  1:93581/149‑1 (MQ=255)                                                                                                                          cagcgcCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGcgcc                         >  1:17420/1‑149 (MQ=255)                                                                                                                                      ggCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGcc               <  1:137539/147‑1 (MQ=255)                                                                                                                                       gCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCccac            >  8:154720/1‑149 (MQ=255)                                                                                                                                          tCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGACCCACCAg         <  2:17420/149‑1 (MQ=255)                                                                                                                                            tttCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCCCACCAGgc       >  2:216646/1‑149 (MQ=255)                                                                                                                                             ttCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCg                                                                     <  1:188556/87‑1 (MQ=255)                                                                                                                                             ttCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCg                                                                     >  2:188556/1‑87 (MQ=255)                                                                                                                                              tCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCCCACCAGgctg     >  3:33441/1‑149 (MQ=255)                                                                                                                                                cAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCCCACCAGgctgtt   <  6:160087/149‑1 (MQ=255)                                                                                                                                                cAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCCCACCAGgctgtt   <  4:212350/149‑1 (MQ=255)                                                                                                                                                 aGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACa                                                              >  2:14525/1‑90 (MQ=255)                                                                                                                                                 aGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACa                                                              <  1:14525/90‑1 (MQ=255)                                                                                                                                                  tcccACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGGATACAGCGCCGCCGGCACCCCGGCGCGGTTCCAGCATGCGCCGCAGCAGGCCCACCAGgctgttg  <  1:216646/148‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     TATGCCATCCGCCCGGAAACGGGCGGCGTCGCCAGGCACAGCTTCAAGGAATAAGCGGTTGCACCTGACTGGCCGGTAGCACGGCCAGATCGCCAAAACGCTGCGCCGAGTCGATGTAGCGCAGCGCCGACTTGGCGTCTTTCCAGCCCACATAGCTCATCAGCGACTTCAGTTCCCAACCATTGGCCGTGGCCCAG‑‑GTGGCGAAGCCACGCCGTAGCGAGTGGCCGGTATACAGCGCCGCCGGCACCCCGGCGC‑GTTCCAGCATGCGCCGCAGCAGGCCCACCAGGCTGTTG  >  NC_002947/3605768‑3606060

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