Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,656,012	Δ1 bp	8.3%	coding (757/4527 nt)	PP_1449 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,656,012	0	G	.	8.3%	190.7 / 11.0	60	coding (757/4527 nt)	PP_1449	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (22/33);  new base . (2/3);  total (24/36)
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.62e-01

CAGCTCGGTCGGGCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTGCTGTCCACGCGC  >  NC_002947/1655866‑1656160                                                                                                                                                   |                                                                                                                                                      cAGCTCGGTCGGGCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCa                                                                                                                                                     >  8:104878/1‑149 (MQ=255)  aGCTCGGTCGGGCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCAt                                                                                                                                                    >  6:301895/1‑149 (MQ=255)    cTCGGTCGGGCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATg‑g                                                                                                                                                 <  1:188568/149‑1 (MQ=255)           gggCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGC‑CATGCGCGCCgg                                                                                                                                           <  6:21404/148‑2 (MQ=255)            ggCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAAc                                                                                                                                         <  6:64822/149‑1 (MQ=255)              cGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAg                                                                                                                                       <  7:104878/149‑1 (MQ=255)                 tCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTTCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCt                                                                                                                                    <  6:308406/149‑1 (MQ=255)                   gCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTgg                                                                                                                                  <  5:301895/149‑1 (MQ=255)                   gCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTg                                                                                                                                   <  1:104298/148‑1 (MQ=255)                   gCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTg                                                                                                                                   <  3:66125/148‑1 (MQ=255)                   gCCACTGCAGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTg                                                                                                                                   >  4:66125/1‑148 (MQ=255)                    ccACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCCCACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGATgg                                                                                                                                  >  3:91604/1‑148 (MQ=255)                    ccACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGc                                                                                                                                 <  5:244851/149‑1 (MQ=255)                          ccGCTGACATGCAGGCCCCCGGCCGATCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCa                                                                                                                             <  3:199203/147‑1 (MQ=255)                                    gCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGC‑CATGCGCGCCgg                                                                                                                                           <  7:243710/123‑2 (MQ=255)                                    gCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGC‑CATGCGCGCCgg                                                                                                                                           >  8:243710/1‑122 (MQ=255)                                      aGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAggg                                                                                                               >  7:298379/1‑149 (MQ=255)                                        gCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGC‑CATGCGCGCCgg                                                                                                                                           >  7:182051/1‑118 (MQ=255)                                        gCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGC‑CATGCGCGCCgg                                                                                                                                           <  8:182051/119‑2 (MQ=255)                                             cGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGc                                                                                                        <  2:108376/149‑1 (MQ=255)                                               gCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAg                                                                                                      <  5:146774/149‑1 (MQ=255)                                                      tGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCATACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGgt                                                                                                    >  4:134141/1‑144 (MQ=255)                                                      tGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCATACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGgt                                                                                                    <  3:134141/144‑1 (MQ=255)                                                        aTGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGTACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGac                                                                                             <  8:319723/149‑1 (MQ=255)                                                        aTGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGac                                                                                             >  1:78584/1‑149 (MQ=255)                                                             gACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTc                                                                                                                       >  8:289322/1‑118 (MQ=255)                                                             gACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTc                                                                                                                       <  7:289322/118‑1 (MQ=255)                                                             gACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCt                                                                                        >  7:48573/1‑149 (MQ=255)                                                               cgccTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATg‑g                                                                                                                                                 <  6:147613/88‑1 (MQ=255)                                                               ccccTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATg‑g                                                                                                                                                 >  5:147613/1‑90 (MQ=255)                                                                cccTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCtt                                                                                                           <  2:131363/127‑2 (MQ=255)                                                                cccTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTg                                                                                                           >  1:131363/1‑127 (MQ=255)                                                                  cTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGgcgcgc                                                                                                                                                      <  4:259673/83‑1 (MQ=255)                                                                  cTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGgcgcgc                                                                                                                                                      >  3:259673/1‑83 (MQ=255)                                                                       gTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGGAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCg                                                                              <  3:63839/149‑1 (MQ=255)                                                                           ccAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTg                                                                          >  3:282477/1‑149 (MQ=255)                                                                           ccAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTg                                                                          >  7:137541/1‑149 (MQ=255)                                                                            cAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCgg                                                                             >  8:1754/1‑145 (MQ=255)                                                                            cAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCgg                                                                             <  7:1754/145‑1 (MQ=255)                                                                              gAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCg                                                                       <  5:250209/149‑1 (MQ=255)                                                                                aTCTGGGAGGCAACCTCTCGCAGACCGGCGGATTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCTCCGAACAGGGTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGa                                                                      <  4:91604/148‑1 (MQ=255)                                                                                     cGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATa                                                                    <  7:315005/145‑1 (MQ=255)                                                                                     cGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACgg                                                                 <  2:78584/148‑1 (MQ=255)                                                                                        ggCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACAACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGTCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGc                                                             <  8:137541/149‑1 (MQ=255)                                                                                                   gCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCg                                                                              <  4:157172/121‑1 (MQ=255)                                                                                                   gCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCg                                                                              >  3:157172/1‑121 (MQ=255)                                                                                                            cgcCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGc                                         >  5:23092/1‑149 (MQ=255)                                                                                                                    gtagtaCTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCt                                 <  8:48643/149‑1 (MQ=255)                                                                                                                       gtaCTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGc                                         >  2:310455/1‑138 (MQ=255)                                                                                                                       gtaCTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGc                                         <  1:310455/138‑1 (MQ=255)                                                                                                                           ttACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCAt                           >  5:336705/1‑148 (MQ=255)                                                                                                                              ccACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAg                        <  4:227113/148‑1 (MQ=255)                                                                                                                               cacaCGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCa                      <  6:324784/149‑1 (MQ=255)                                                                                                                                acacGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGc                       >  8:2912/1‑147 (MQ=255)                                                                                                                                 cacGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGc                       >  5:226021/1‑146 (MQ=255)                                                                                                                                  acGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGc                   >  4:332807/1‑149 (MQ=255)                                                                                                                                   cGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCaa                     <  6:226021/146‑1 (MQ=255)                                                                                                                                      aTACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCttt               <  3:332807/149‑1 (MQ=255)                                                                                                                                      aTACGCCCGGCGCGCATG‑GCGCCGAACACGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCttt               <  4:282477/149‑1 (MQ=255)                                                                                                                                       tACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTg              <  8:48573/149‑1 (MQ=255)                                                                                                                                          gCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTAGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTGCTg           >  4:256703/1‑149 (MQ=255)                                                                                                                                            ccGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTGCTGTc         <  3:145635/149‑1 (MQ=255)                                                                                                                                                   gcATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTGCTGTCCAcgcgc  >  4:317943/1‑149 (MQ=255)                                                                                                                                                   gcATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTGCTGTCCAcgcgc  >  4:127422/1‑149 (MQ=255)                                                                                                                                                   |                                                                                                                                                      CAGCTCGGTCGGGCGGTCGCCACTGCCGCTGACATGCAGGCCCCCGGCCGAGCTGATGTCGACCCCTTCGGTGCCCAGAATCTGCGAGGCCCCCACACGCACACCCGCGCCTTCGGTAGTACTTACCACACGGATACGCCCGGCGCGCATG‑GCGCCGAACAGGCTGGCATCAATGCTCGTGGAGGGGCTGGGCAGGTGTTCGACCACCTCGCCGCTGCGGCTGTCGATACGGTTGCGGCCGACAGTGATGTCGAGCGTTTCGCTGGCCATGAGCAAGCCTTTGCTGTCCACGCGC  >  NC_002947/1655866‑1656160

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