Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,597,888	G→A	11.6%	P463P (CCG→CCA)	PP_4923 →	outer membrane efflux protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,597,888	0	G	A	11.6%	155.0 / 9.1	61	P463P (CCG→CCA)	PP_4923	outer membrane efflux protein
Reads supporting (aligned to +/- strand):  ref base G (28/25);  new base A (4/3);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.07e-01

CAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAA  >  NC_002947/5597750‑5598036                                                                                                                                           |                                                                                                                                                      caacaaTAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACC‑GAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                             <  6:196540/148‑1 (MQ=255)    caaTAGCCGTTACGACTACATTCTCGACACCCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGgcag                                                                                                                                          <  2:36196/149‑1 (MQ=255)         gCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCagcagc                                                                                                                                      <  2:58585/148‑1 (MQ=255)         gCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCACGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCg                                                                                                                                     <  6:18530/149‑1 (MQ=255)          ccGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCgg                                                                                                                                    <  6:225487/149‑1 (MQ=255)            gTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGc                                                                                                                                   >  2:270347/1‑148 (MQ=255)             ttACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGcc                                                                                                                                  >  1:266177/1‑148 (MQ=255)                   cTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCCGCTGCTGCCGCCaggtccg                                                                                                                           <  5:113480/149‑7 (MQ=255)                     aCATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                         >  7:25532/1‑149 (MQ=255)                      cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                         <  8:25532/148‑1 (MQ=255)                      cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTc                                                                                                                         <  2:266177/148‑1 (MQ=255)                      cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                          <  4:6410/131‑1 (MQ=255)                      cATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgcag                                                                                                                                          >  3:6410/1‑131 (MQ=255)                         tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCCGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAgc                                                                                                                     >  4:127426/1‑149 (MQ=255)                         tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                           >  2:239439/1‑127 (MQ=255)                         tctcGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                           <  1:239439/127‑1 (MQ=255)                             gACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAg                                                                                                                 <  5:134370/149‑1 (MQ=255)                                  ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                       >  6:108853/1‑138 (MQ=255)                                  ccTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGa                                                                                                                       <  5:108853/138‑1 (MQ=255)                                      aGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCtt                                                                                                            >  2:153067/1‑145 (MQ=255)                                            aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACACCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGAAAgggg                                                                                                  >  3:207729/1‑149 (MQ=255)                                            aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCATGCCTTGAGATAgg                                                                                                    <  8:148856/147‑1 (MQ=255)                                            aaGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAgg                                                                                                    <  2:262123/147‑1 (MQ=255)                                             tcacggcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGc                                                                                                 >  4:95698/5‑149 (MQ=255)                                                  gcagCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGAc                                                                                            >  8:71559/1‑149 (MQ=255)                                                       cGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                             >  8:34865/1‑95 (MQ=255)                                                       cGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTggcgg                                                                                                                                             <  7:34865/95‑1 (MQ=255)                                                         gCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTcgc                                                                                     >  4:210738/1‑149 (MQ=255)                                                             ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGAc                                                                                                                                                    <  6:276018/82‑1 (MQ=255)                                                             ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGAc                                                                                                                                                    >  5:276018/1‑82 (MQ=255)                                                             ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCgg                                                                                   <  4:27662/147‑1 (MQ=255)                                                             ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCgg                                                                                   >  3:27662/1‑147 (MQ=255)                                                             ccTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGc                                                                                 >  5:98798/1‑149 (MQ=255)                                                                    ccGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCg                                                                          >  7:63629/1‑149 (MQ=255)                                                                      gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑Ct                                                                       >  8:216040/1‑149 (MQ=255)                                                                      gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                           <  4:63041/82‑1 (MQ=255)                                                                      gCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCAGGTCCGGCGgca                                                                                                                                           >  3:63041/1‑82 (MQ=255)                                                                       cAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTGCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTc                                                                      <  6:214222/149‑1 (MQ=255)                                                                                cAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACgg                                                              <  6:98798/148‑1 (MQ=255)                                                                                   gACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGcc                                                          >  4:56969/1‑149 (MQ=255)                                                                                      cTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAgg                                                       >  2:218322/1‑149 (MQ=255)                                                                                       tGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCGTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGc                                                      >  6:171836/1‑149 (MQ=255)                                                                                        gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGcc                                                     >  6:173426/1‑149 (MQ=255)                                                                                        gAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGcc                                                     >  2:266463/1‑149 (MQ=255)                                                                                                 ccTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGc                                              <  7:216040/147‑1 (MQ=255)                                                                                                  cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCGCCTCCCACACGGTGCCAGGCCTGAAGGCat                                            <  1:153067/149‑1 (MQ=255)                                                                                                  cTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCat                                            >  7:147379/1‑148 (MQ=255)                                                                                                        ccTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAc                                     >  8:68624/1‑149 (MQ=255)                                                                                                         cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                    >  4:119758/1‑149 (MQ=255)                                                                                                         cTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGACCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCAcc                                    >  2:167644/1‑149 (MQ=255)                                                                                                            aCTACAACCCCGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCTCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCtgt                                 <  3:127426/149‑1 (MQ=255)                                                                                                                cAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTga                               <  3:56969/147‑1 (MQ=255)                                                                                                                 aaCCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGagag                            >  3:65641/1‑149 (MQ=255)                                                                                                                    ccGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                         <  3:210738/149‑1 (MQ=255)                                                                                                                        aCAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCGTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCtt                     <  5:171836/149‑1 (MQ=255)                                                                                                                          aaGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGc                   >  4:109718/1‑149 (MQ=255)                                                                                                                                ttCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGa             >  1:30903/1‑149 (MQ=255)                                                                                                                                  ccTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGAt            <  1:266463/148‑1 (MQ=255)                                                                                                                                   cTGCCGCCAGGTCCGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAg          >  6:185341/1‑149 (MQ=255)                                                                                                                                    tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                         <  3:221126/133‑1 (MQ=255)                                                                                                                                    tGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTg                         >  4:221126/1‑133 (MQ=255)                                                                                                                                          cgGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCa   <  4:65641/149‑1 (MQ=255)                                                                                                                                           ggACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCaa  <  8:63629/149‑1 (MQ=255)                                                                                                                                           |                                                                                                                                                      CAACAATAGCCGTTACGACTACATTCTCGACAACCTGAGCCTGAAGCAGGCAGCCGGCACCCTGAGCCCGCAAGACCTGCAGGACCTGAAGCGCTACCTGAAACCTGACTACAACCCGGACAAGGACTTCCTGCCGCCGGACCTGGCGGCAGCAGCGGCCAAGAACTTCGAGCGCAGGCCTTGAGATAGGGGCTGGACTCTTCGCGGGCTTGCCCG‑CTCCCACACGGTGCCAGGCCTGAAGGCATGATGCACCTGTGAGAGCTGGCTTGCTGGCGATAGGCCTTCAA  >  NC_002947/5597750‑5598036

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