Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,621,512	T→G	6.2%	S176R (AGC→CGC)	PP_3192 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,621,512	0	T	G	6.2%	203.2 / 3.8	64	S176R (AGC→CGC)	PP_3192	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/30);  new base G (2/2);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

TCCACCATAGAATAATATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCGTTGCGG  >  NC_002947/3621370‑3621654                                                                                                                                                |                                                                                                                                               tCCACCATAGAATAATATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGccc                                                                                                                                          >  2:184482/1‑149 (MQ=255)   caccaTAGAATAATATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCAGCGGCGccga                                                                                                                                         >  1:181893/1‑146 (MQ=255)   caccaTAGAATAATATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCAGCGGCGccga                                                                                                                                         <  2:181893/148‑3 (MQ=255)      caTAGAATAATATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCg                                                                                                                                      >  8:197733/1‑148 (MQ=255)           aataatATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCa                                                                                                                                >  3:149401/1‑149 (MQ=255)                  ttCCACACACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCgg                                                                                                                         <  6:107703/149‑1 (MQ=255)                  ttCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCgg                                                                                                                         <  3:196155/149‑1 (MQ=255)                     cacaAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCATCGTAATCGGCGCCGCTGCGCCCCTCGAATGCAGCCTCGGCAt                                                                                                                      <  3:291660/149‑1 (MQ=255)                     cacaAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCAt                                                                                                                      >  2:253364/1‑149 (MQ=255)                     cacaAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCAt                                                                                                                      >  2:240280/1‑149 (MQ=255)                     cacaAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCAt                                                                                                                      <  7:197733/149‑1 (MQ=255)                      acaAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCAGCGGCGccga                                                                                                                                         <  6:64627/129‑3 (MQ=255)                      acaAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCAGCGGCGccga                                                                                                                                         >  5:64627/1‑127 (MQ=255)                        aaaCAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTTGCCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATg                                                                                                                     >  1:276923/1‑148 (MQ=255)                        aaaCAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGc                                                                                                                    >  5:311925/1‑148 (MQ=255)                         aaCAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCaa                                                                                                                  <  4:149401/149‑1 (MQ=255)                         aaCAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCaa                                                                                                                  <  4:139347/149‑1 (MQ=255)                             gCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGcc                                                                                                              <  3:235539/149‑1 (MQ=255)                                    cGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACaa                                                                                                       <  1:240280/149‑1 (MQ=255)                                        tgCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGt                                                                                                    <  1:184482/148‑1 (MQ=255)                                          ccTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGc                                                                                                 >  1:348266/1‑149 (MQ=255)                                            ttCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTAATTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGt                                                                                               >  1:252902/1‑149 (MQ=255)                                               aggcaggcAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGttt                                                                                             <  3:276229/148‑1 (MQ=255)                                                ggcaggcaATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTg                                                                                                                                               <  2:269887/97‑1 (MQ=255)                                                ggcaggcaATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTg                                                                                                                                               >  1:269887/1‑97 (MQ=255)                                                 gcaggcaATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACATAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTc                                                                                          <  3:56992/149‑1 (MQ=255)                                                   aggcaATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCtt                                                                                        >  3:255080/1‑149 (MQ=255)                                                      caATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTAAACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGAGCCCCTCGAAAGCACCATCGGCATGCAAGGCCAGGACAAGGATGCGTTTTTCTATCa                                                                                     >  6:294993/1‑149 (MQ=255)                                                           ccGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCGTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTc                                                                                <  1:253364/149‑1 (MQ=255)                                                              ccTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCa                                                                               <  7:135329/147‑1 (MQ=255)                                                              ccTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCa                                                                               >  8:135329/1‑147 (MQ=255)                                                               cTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGc                                                                                                                    >  6:206107/1‑109 (MQ=255)                                                                  ggCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTc                                                                          <  4:255080/148‑1 (MQ=255)                                                                    cccTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGc                                                                       <  7:16705/149‑1 (MQ=255)                                                                      cTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTg                                                                                            <  3:28574/126‑1 (MQ=255)                                                                      cTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTg                                                                                            >  4:28574/1‑126 (MQ=255)                                                                                ttcAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGtt                                                            >  4:237069/1‑149 (MQ=255)                                                                                ttcAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGt                                                             >  5:269731/1‑148 (MQ=255)                                                                                  cAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGtt                                                            <  6:269731/147‑1 (MQ=255)                                                                                   aGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTAt                                                          <  6:311925/148‑1 (MQ=255)                                                                                        tttGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCt                                                    >  8:194764/1‑149 (MQ=255)                                                                                        tttGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCt                                                    >  7:70914/1‑149 (MQ=255)                                                                                           gTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTg                                                   >  1:15925/1‑147 (MQ=255)                                                                                            tACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGtt                                                            >  1:325938/1‑137 (MQ=255)                                                                                            tACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGtt                                                            <  2:325938/137‑1 (MQ=255)                                                                                               acGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCa                                             <  2:164492/149‑1 (MQ=255)                                                                                                     gcccgcccCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTg                                       <  8:70914/149‑1 (MQ=255)                                                                                                       ccgcccCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTc                                                                                          <  1:353700/96‑1 (MQ=255)                                                                                                       ccgcccCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTc                                                                                          >  2:353700/1‑96 (MQ=255)                                                                                                       ccgcccCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCGGCCACGCTTGCa                                     <  2:276923/149‑1 (MQ=255)                                                                                                           cccATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAgg                                   >  6:240664/1‑147 (MQ=255)                                                                                                           cccATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGaa                                 <  2:252902/149‑1 (MQ=255)                                                                                                            ccATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGTATGCACCCTGGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAg                                <  5:294993/149‑1 (MQ=255)                                                                                                                   gAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGa                         >  8:91789/1‑149 (MQ=255)                                                                                                                           tGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTcc                 >  1:308259/1‑149 (MQ=255)                                                                                                                               aaCGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAg                                <  3:144851/130‑1 (MQ=255)                                                                                                                               aaCGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAg                                >  4:144851/1‑130 (MQ=255)                                                                                                                                    aaTCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCa                                                 >  3:109058/1‑108 (MQ=255)                                                                                                                                    aaTCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCa                                                 <  4:109058/108‑1 (MQ=255)                                                                                                                                    aaTCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCg        >  4:203679/1‑149 (MQ=255)                                                                                                                                        ggCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCGTTg     <  7:91789/148‑1 (MQ=255)                                                                                                                                        ggCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCGTTGc    <  1:1069/149‑1 (MQ=255)                                                                                                                                        ggCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCGTTGc    >  8:332478/1‑149 (MQ=255)                                                                                                                                           gccgcTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCGTTGCgg  >  1:292960/1‑148 (MQ=255)                                                                                                                                                |                                                                                                                                               TCCACCATAGAATAATATTCCACAAACAGCAAGAGCGTGTGCCTTCAGGCAGGCAATACCGCCTAGGCCCTGACTT‑CCTTCAGGTATTTGTACACGGTGGCCCGCCCCATGCTGAGGACGTTGGCAACGTAATCGGCGCCGCTGCGCCCCTCGAATGCACCCTCGGCATGCAAGGCCAGGACAAGGTTGCGTTTGTCTTTCATCTTCAGGGTCTGCAGCGACAGGTTATGCGCCTGCAGCCACGCTTGCAGGAAGGTGTTGATGCGCTCCTGCCAATCGTTGCGG  >  NC_002947/3621370‑3621654

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