Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,235,720	G→A	5.8%	A352T (GCC→ACC)	fecA →	outer membrane ferric citrate porin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,235,720	0	G	A	5.8%	205.3 / 2.1	69	A352T (GCC→ACC)	fecA	outer membrane ferric citrate porin
Reads supporting (aligned to +/- strand):  ref base G (32/33);  new base A (2/2);  total (34/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.65e-01

CCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGCTTCTTCACCAAGACCCTGCG  >  NC_002947/5235572‑5235862                                                                                                                                                     |                                                                                                                                                ccGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCAccg                                                                                                                                                 >  5:299326/1‑149 (MQ=255)       gggTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGAtt                                                                                                                                           <  6:297348/149‑1 (MQ=255)        ggTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGa                                                                                                                                             <  1:130598/146‑1 (MQ=255)        ggTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGa                                                                                                                                             >  2:130598/1‑146 (MQ=255)            aCTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGAt                                                                                                                                      <  6:299326/149‑1 (MQ=255)                gcgcGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGcc                                                                                                                                   >  5:115716/1‑148 (MQ=255)                      acacAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCg                                                                                                                             >  4:50031/1‑148 (MQ=255)                          aGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTAc                                                                                                                          >  4:171718/1‑147 (MQ=255)                          aGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTAc                                                                                                                          >  4:114453/1‑147 (MQ=255)                          aGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCa                                                                                                                        <  4:302328/149‑1 (MQ=255)                              acacGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGCCATGCCAGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTcc                                                                                                                    <  8:237002/149‑1 (MQ=255)                                acGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCAc                                                                                                                  >  7:206448/1‑149 (MQ=255)                                      aTCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTAATACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGAATTCGATGCCGACCCGTACCAGTCAACGCTCCt                                                                                                            >  4:231615/1‑149 (MQ=255)                                      aTCGACGACCTGATCCCCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCt                                                                                                            >  2:75426/1‑149 (MQ=255)                                         gacgacCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGa                                                                                                          >  4:260279/1‑148 (MQ=255)                                          acgacCTGATCCTCACGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGaaa                                                                                                        <  1:170534/149‑1 (MQ=255)                                          acgacCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGaaa                                                                                                        <  3:308294/149‑1 (MQ=255)                                          acgacCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGa                                                                                                          >  2:197780/1‑147 (MQ=255)                                              ccTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACa                                                                                                    <  3:50031/149‑1 (MQ=255)                                               cTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACa‑a                                                                                                  <  3:190441/149‑1 (MQ=255)                                                 gATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGaaa                                                                                                        >  2:240983/1‑142 (MQ=255)                                                  aTCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑Ag                                                                                                 <  5:89304/147‑1 (MQ=255)                                                  aTCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGtt                                                                                               >  4:281671/1‑149 (MQ=255)                                                      tCAAGGGCAAGCTGCATCTGGACGAAGCCAACAGCCTGCACTCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGACCGGCGGCCTGAGAACCGCCGATTTCGATGCCGACCAGTACGATTCGACGCGCGTGAAAGTCT‑AGGTCTgg                                                                                           >  6:14638/1‑149 (MQ=255)                                                      tCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTgg                                                                                           >  7:193389/1‑149 (MQ=255)                                                      tCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTgg                                                                                           >  8:97026/1‑149 (MQ=255)                                                        aaGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTgggg                                                                                         >  8:63476/1‑149 (MQ=255)                                                        aaGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTgg                                                                                           <  6:115716/147‑1 (MQ=255)                                                         aGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGgc                                                                                        >  8:266198/1‑149 (MQ=255)                                                            gCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCACAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGgccgc                                                                                     >  5:308526/1‑149 (MQ=255)                                                                gCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCa                                                                                 <  8:193389/149‑1 (MQ=255)                                                                 cTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCAc                                                                                <  3:231615/149‑1 (MQ=255)                                                                  tGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACg                                                                               >  8:43762/1‑149 (MQ=255)                                                                  tGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACg                                                                               >  2:287789/1‑149 (MQ=255)                                                                    cAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGAccc                                                                                                                              >  2:160321/1‑101 (MQ=255)                                                                    cAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGAccc                                                                                                                              <  1:160321/101‑1 (MQ=255)                                                                    cAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCt                                                                             >  4:210447/1‑149 (MQ=255)                                                                     aaCTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTg                                                                            >  8:60892/1‑149 (MQ=255)                                                                      aCTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGATAGACA‑AGTTCTGGGGCCGCCGCACGCTGt                                                                           <  1:197780/149‑1 (MQ=255)                                                                          gACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTc                                                                         <  7:97026/147‑1 (MQ=255)                                                                             gAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCAGCACGCTGTTCAACtt                                                                    >  3:85583/1‑149 (MQ=255)                                                                             gAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATCTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACtt                                                                    >  8:149933/1‑149 (MQ=255)                                                                              aaGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTc                                                                   <  3:114453/149‑1 (MQ=255)                                                                               aGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCg                                                                  <  8:73598/149‑1 (MQ=255)                                                                                 ccAACAGCCTGTACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGc                                                                <  1:287789/149‑1 (MQ=255)                                                                                 ccAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATCTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGc                                                                <  7:149933/149‑1 (MQ=255)                                                                                      aGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCTACTTCGGCTACGa                                                           <  6:308526/149‑1 (MQ=255)                                                                                           gCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTAc                                                       >  1:303755/1‑148 (MQ=255)                                                                                            cACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCACCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTAc                                                       >  1:206149/1‑147 (MQ=255)                                                                                             aCGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACACAGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACaa                                                     <  1:75426/149‑1 (MQ=255)                                                                                               gCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAg                                                    >  7:307599/1‑147 (MQ=255)                                                                                                ccATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAg                                                 <  7:60892/149‑1 (MQ=255)                                                                                                ccATGGACCAGTACTACGATCGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCCTACCAGTCCACTCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTCTTCAACTTCGGCTACGACTACAACCAg                                                 <  5:14638/149‑1 (MQ=255)                                                                                                 cATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAgg                                                <  7:249979/149‑1 (MQ=255)                                                                                                        cAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCg                                         >  6:36004/1‑149 (MQ=255)                                                                                                         aGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGc                                        <  3:210447/149‑1 (MQ=255)                                                                                                          gTACTACGAGGGCGAGGCCGCCATGCCCCGCGGCATGAGCACCACCGATTTCGATGCCGACCCGTACAAGTCCACGCGCCTGAAATACA‑AGTTCTGGGGCCGCCGAACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGc                                        >  1:189824/1‑148 (MQ=255)                                                                                                             ctacGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAg                                     <  3:281671/148‑1 (MQ=255)                                                                                                             ctacGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGACGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGt                                    >  3:197927/1‑149 (MQ=255)                                                                                                                cGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCa                                 <  3:260279/149‑1 (MQ=255)                                                                                                                     gCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTc                            <  5:36004/149‑1 (MQ=255)                                                                                                                       gAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCaa                          <  8:206448/149‑1 (MQ=255)                                                                                                                          gCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAg                       <  3:171718/149‑1 (MQ=255)                                                                                                                              aCATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGcttc                   <  7:43762/149‑1 (MQ=255)                                                                                                                              aCATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGcttc                   <  7:266198/149‑1 (MQ=255)                                                                                                                              aCATGCCCGGCGGCCTGAGCACCACCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACATcttc                   <  2:189824/149‑1 (MQ=255)                                                                                                                                    ccGGCGGCCTGAGCACCACCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGCTTCTTCAcc             <  2:206149/149‑1 (MQ=255)                                                                                                                                     cggcggCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGCTTCTTCAcc             >  4:240943/1‑148 (MQ=255)                                                                                                                                           ccTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGCTTCTTCACCAAGAccc      >  1:309380/1‑149 (MQ=255)                                                                                                                                                 cACCGCCGATTTCGATGCCGAACCGTACCAGTCCACGAGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGAATACAAGAAGGATGACAGAAAGTTCAGAGTCAACAGATTATTAACCAATACCCTgcg  >  4:22592/1‑147 (MQ=255)                                                                                                                                                     |                                                                                                                                                CCGCGGGGGTGACTGGCGCGAACACAGCGACACGCAGATCGACGACCTGATCCTCAAGGGCAAGCTGCAACTGGACGAAGCCAACAGCCTGCACGCCATGGCCCAGTACTACGAGGGCGAGGCCGACATGCCCGGCGGCCTGAGCACCGCCGATTTCGATGCCGACCCGTACCAGTCCACGCGCCTGAAAGACA‑AGTTCTGGGGCCGCCGCACGCTGTTCAACTTCGGCTACGACTACAAGCAGGATGACCGCCAGTTCAGCGTCAACAGCTTCTTCACCAAGACCCTGCG  >  NC_002947/5235572‑5235862

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