Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,183,003	A→T	11.0%	I244N (ATC→AAC)	PP_4563 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,183,003	0	A	T	11.0%	167.4 / 8.9	61	I244N (ATC→AAC)	PP_4563	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/26);  new base T (3/4);  total (31/30)
Fisher's exact test for biased strand distribution p-value = 7.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTC  >  NC_002947/5182855‑5183148                                                                                                                                                     |                                                                                                                                                    cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                                     <  5:345910/149‑1 (MQ=255) cGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGa                                                                                                                                                     <  3:252362/149‑1 (MQ=255)      agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCaa                                                                                                                                                 <  7:57472/148‑1 (MQ=255)      agcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAAc                                                                                                                                                <  2:272370/149‑1 (MQ=255)       gcTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACg                                                                                                                                               >  7:21781/1‑149 (MQ=255)             cTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATTGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGt                                                                                                                                         <  5:2517/149‑1 (MQ=255)                ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                                           >  7:169151/1‑144 (MQ=255)                ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAgg                                                                                                                                           <  8:169151/144‑1 (MQ=255)                ccGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCAc                                                                                                                                      >  1:12718/1‑149 (MQ=255)                  gCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCActc                                                                                                                                    <  6:332457/149‑1 (MQ=255)                    aaCTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCtt                                                                                                                                  <  2:282435/149‑1 (MQ=255)                         cGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTc                                                                                                                               <  1:4652/147‑1 (MQ=255)                              aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACAcc                                                                                                                       >  6:17491/1‑149 (MQ=255)                              aGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACAcc                                                                                                                       <  8:21781/149‑1 (MQ=255)                                   tAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCaa                                                                                                                  <  7:135729/149‑1 (MQ=255)                                        ccGTCAGGCTTTCTCGATGCGATCGTCATGGATTACGATCAGGCCCTGCTTGAACAACCAGCCAATGGCCTTCTTGAAGTTGCACTTGCTGACGTTGAACAGCTTGCTGATCAACGAAGGGTCACTCTTGTCGC‑ACAACCCCAacacac                                                                                                             >  3:157906/1‑147 (MQ=255)                                          gTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCg                                                                                                           >  8:35066/1‑149 (MQ=255)                                             aGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCt                                                                                                        >  7:313002/1‑149 (MQ=255)                                                cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCat                                                                                                                  >  7:1188/1‑135 (MQ=255)                                                cTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCat                                                                                                                  <  8:1188/136‑2 (MQ=255)                                                       gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCtt                                                                                                 <  8:313002/146‑1 (MQ=255)                                                       gATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCtt                                                                                                 <  5:17491/146‑1 (MQ=255)                                                        aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑Ag                                                                                            <  5:45611/149‑1 (MQ=255)                                                        aTGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑Ag                                                                                            <  2:132602/149‑1 (MQ=255)                                                          gCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCACCACCCCGCCTTCGGCTTCC‑AGgc                                                                                          >  8:41317/1‑149 (MQ=255)                                                            gATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGgcgc                                                                                        <  3:313289/149‑1 (MQ=255)                                                                gTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGcc                                                                                     <  5:81872/148‑1 (MQ=255)                                                                 tCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCAt                                                                                   >  7:22412/1‑149 (MQ=255)                                                                      gATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATc                                                                               <  8:22412/148‑1 (MQ=255)                                                                      gATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCt                                                                              >  4:35744/1‑149 (MQ=255)                                                                            gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                                <  7:227845/79‑1 (MQ=38)                                                                            gATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGATCAGCAAGCTGTTCAAc                                                                                                                                                >  8:227845/1‑79 (MQ=38)                                                                                    ccTGCTTGAACAACCCGCCAATGGCCTTCTTGATGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGAAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                                >  4:53250/1‑149 (MQ=255)                                                                                    ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                   >  4:290153/1‑146 (MQ=255)                                                                                    ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                   <  3:290153/146‑1 (MQ=255)                                                                                    ccTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGgtg                                                                >  3:163669/1‑149 (MQ=255)                                                                                           gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGGCTTCGGCTTCCAAGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                             <  1:41846/146‑1 (MQ=255)                                                                                           gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                             >  4:164412/1‑145 (MQ=255)                                                                                           gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                             >  2:41846/1‑145 (MQ=255)                                                                                           gAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTcc                                                             <  3:164412/145‑1 (MQ=255)                                                                                              caacCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAggg                                                      <  7:152127/149‑1 (MQ=255)                                                                                                  ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                                           >  4:37135/1‑140 (MQ=255)                                                                                                  ccGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGc                                                           <  3:37135/140‑1 (MQ=255)                                                                                                            tcttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGc                                        >  3:115663/2‑149 (MQ=255)                                                                                                              ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTg                                      >  1:323281/1‑149 (MQ=255)                                                                                                              ttcttGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACAACCCGCCTTCGGATTAC‑AGGCTCGCCATGATATGCTCTTGCAGGGTGTCCGCATCGGCTGCGCCTACAGGCt                                       >  1:227391/1‑148 (MQ=255)                                                                                                                  tGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                            >  1:54764/1‑139 (MQ=255)                                                                                                                  tGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTAc                                            <  2:54764/139‑1 (MQ=255)                                                                                                                      gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                   >  1:264821/1‑112 (MQ=255)                                                                                                                      gTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAgg                                                                   <  2:264821/112‑1 (MQ=255)                                                                                                                       ttGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAg                             >  5:278332/1‑149 (MQ=255)                                                                                                                                cctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCtt                    >  4:239915/5‑149 (MQ=255)                                                                                                                                 tGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTg                   >  7:122502/1‑149 (MQ=255)                                                                                                                                 ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAgg                                  <  4:210065/131‑1 (MQ=255)                                                                                                                                 ctgcGTTGATCAGCAAGCTGTTCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAgg                                  >  3:210065/4‑134 (MQ=255)                                                                                                                                   aCGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGcc                 <  3:35744/149‑1 (MQ=255)                                                                                                                                      ttGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTc              <  2:12718/149‑1 (MQ=255)                                                                                                                                         aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                                       >  8:6805/1‑105 (MQ=255)                                                                                                                                         aaCAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAgg                                                       <  7:6805/105‑1 (MQ=255)                                                                                                                                           cAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGCAACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTATTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTgg          >  1:51125/1‑149 (MQ=255)                                                                                                                                                tGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACt    >  6:126397/1‑149 (MQ=255)                                                                                                                                                  cTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTc  >  4:323969/1‑149 (MQ=255)                                                                                                                                                  cTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGAGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTc  <  6:150506/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CGAGCAGCTGGACTTCCGCAACTTCGTACAGGCGTAAGGCCGTCAGGCTTTCTCGATGCGATCGTCATGGATGACGATCAGGCCCTGCTTGAACAACCCGCCAATGGCCTTCTTGAAGTTGCCCTTGCTGACGTTGAACAGCTTGCTGATCAACGCAGGGTCACTCTTGTCGC‑ACACCCCCAACACCCCGCCTTCGGCTTCC‑AGGCGCGCCATGATCTGCTCTTGCAGGGTGTCCGCCAGGGCGGCGCCTACCGGCTGCAGGCTCAGTGAAATCTTGCCGTCCTGGCGCACTTC  >  NC_002947/5182855‑5183148

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