Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,591,577	T→A	9.2%	L258* (TTG→TAG) ‡	PP_3169 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,591,577	0	T	A	9.2%	185.9 / 8.4	65	L258* (TTG→TAG) ‡	PP_3169	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (29/30);  new base A (4/2);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.85e-01

GTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGT  >  NC_002947/3591429‑3591718                                                                                                                                                      |                                                                                                                                              gTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAATGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCtt                                                                                                                                               <  5:114158/149‑1 (MQ=255)         ggTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCAAGCGCTGGt                                                                                                                                       <  7:116491/149‑1 (MQ=255)          gTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGtt                                                                                                                                      <  1:137331/149‑1 (MQ=255)             cTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAAGGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTGAc                                                                                                                                   <  6:52322/149‑1 (MQ=255)                gCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGa                                                                                                                                <  8:210043/149‑1 (MQ=255)                gCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGa                                                                                                                                >  5:322173/1‑149 (MQ=255)                    aGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGAt                                                                                                                               >  3:348481/1‑146 (MQ=255)                    aGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGAt                                                                                                                               <  4:348481/146‑1 (MQ=255)                                gCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                >  6:66316/1‑149 (MQ=255)                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                <  5:383347/141‑1 (MQ=255)                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTAt                                                                                                                >  6:383347/1‑141 (MQ=255)                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                        >  8:283521/1‑149 (MQ=255)                                        tacGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGAt                                                                                                        >  3:234674/1‑149 (MQ=255)                                          cGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTg                                                                                                      >  6:350011/1‑149 (MQ=255)                                                       agcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCCACGCACtggt                                                                                         <  7:283521/149‑1 (MQ=255)                                                        gcGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACtggtg                                                                                        <  5:350011/149‑1 (MQ=255)                                                          gTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTggtggt                                                                                      <  5:66316/149‑1 (MQ=255)                                                           tGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTc                                                                                     <  2:41838/149‑1 (MQ=255)                                                              tCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGt                                                                                                                   <  7:88713/116‑1 (MQ=255)                                                              tCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGt                                                                                                                   >  8:88713/1‑116 (MQ=255)                                                                gACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGa                                                                                >  8:351504/1‑149 (MQ=255)                                                                   cAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCt                                                                             <  6:2637/149‑1 (MQ=255)                                                                     gCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgc                                                                           <  6:322173/149‑1 (MQ=255)                                                                      cAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgcg                                                                          <  2:219335/149‑1 (MQ=255)                                                                       aaCCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTgcgc                                                                         <  4:219081/149‑1 (MQ=255)                                                                           accCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGAAGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                     >  6:207860/1‑149 (MQ=255)                                                                            cccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                     <  6:105902/148‑1 (MQ=255)                                                                            cccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgg                                                                     >  5:105902/1‑148 (MQ=255)                                                                             ccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgggg                                                                   <  1:183044/149‑1 (MQ=255)                                                                             ccGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATggg                                                                    >  6:39765/1‑148 (MQ=255)                                                                              cGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGa                                                                  >  8:225583/1‑149 (MQ=255)                                                                              cGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGa                                                                  >  3:79847/1‑149 (MQ=255)                                                                                aaaTGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATc                                                                >  1:318192/1‑149 (MQ=255)                                                                                    gCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGCCCTGCGCATGGGGATCCCgg                                                            >  5:59181/1‑149 (MQ=255)                                                                                     cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCgg                                                            <  2:318192/148‑1 (MQ=255)                                                                                     cAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGAAGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCgg                                                            <  5:207860/148‑1 (MQ=255)                                                                                      aaCGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCggg                                                           >  7:52855/1‑148 (MQ=255)                                                                                           ttGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAAtt                                                      >  7:234508/1‑148 (MQ=255)                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTaca                                                   >  7:35247/1‑149 (MQ=255)                                                                                             gTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTa                                                     <  7:225583/147‑1 (MQ=255)                                                                                              tCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCg                                                                                 >  6:81900/1‑118 (MQ=255)                                                                                              tCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCg                                                                                 <  5:81900/118‑1 (MQ=255)                                                                                                  gcAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCttt                                              <  7:351504/149‑1 (MQ=255)                                                                                                  gcAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCttt                                              <  4:234674/149‑1 (MQ=255)                                                                                                     aCCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGc                                           <  8:52855/149‑1 (MQ=255)                                                                                                        tGGAACTGGCGCAAGCACTGAACGATGACCCGTTGGCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGCCCAGATTGGCGACGCACTggtggt                                                                                      >  1:247839/1‑104 (MQ=255)                                                                                                        tGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTggtggt                                                                                      <  2:247839/103‑1 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCAc                                                                                   <  3:260827/98‑1 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCAc                                                                                   >  4:260827/1‑98 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgggg                                                                   <  3:375275/114‑1 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATgggg                                                                   >  4:375275/1‑114 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTacac                                                  <  5:243343/131‑1 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTacac                                                  >  6:243343/1‑131 (MQ=255)                                                                                                                gcgcAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATgg                                 <  8:35247/148‑1 (MQ=255)                                                                                                                         cTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGa                       <  8:234508/149‑1 (MQ=255)                                                                                                                            aaCGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatca                    >  3:50233/1‑149 (MQ=255)                                                                                                                            aaCGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGatca                    >  4:302078/1‑149 (MQ=255)                                                                                                                                       cggtgaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAg          >  3:221653/7‑149 (MQ=255)                                                                                                                                        ggtgaaCCAGCGCAAGCGCTGGTGCAACGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGt         >  5:352248/6‑149 (MQ=255)                                                                                                                                        ggtgaaCCAGCGCAAGCGCTGGTGCAACGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGt         <  6:352248/144‑1 (MQ=255)                                                                                                                                           gaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGa      >  8:240587/3‑149 (MQ=255)                                                                                                                                           gaaCCAGCGCAAGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGa      >  6:286148/3‑149 (MQ=255)                                                                                                                                           gCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGa      >  4:315027/1‑149 (MQ=255)                                                                                                                                            cACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGAt     <  5:39765/149‑1 (MQ=255)                                                                                                                                             aCCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATg    >  5:46165/1‑149 (MQ=255)                                                                                                                                               cAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGt  >  1:80468/1‑149 (MQ=255)                                                                                                                                                     ttGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGa                                                                                >  2:113133/1‑65 (MQ=255)                                                                                                                                                      |                                                                                                                                              GTCTGGCGGGTGCTGGCCAAGACCCCGAATGGCGATTACTACGAAGGCATCAGCAGCGTGCTCGACCAGCAACCACCCGAAATGCAACGCTTGTCGCGCAACCTGGAACTGGCGCAAGCACTGAACGATGACCCGTT‑GCACCAGCGCTTGCGCTGGTTCACCGATGACTGGCTGCGTTATGACCAGATTGGCGACGCACTGGTGGTCACCGACCTGCGCATGGGGATCCCGGGGAATTACACCTTTCGCTTCAACATGGCTCAGCGCGATCATCAAGGGCAGTGGATGGT  >  NC_002947/3591429‑3591718

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