Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,182,610	C→G	6.1%	R18G (CGG→GGG)	PP_4562 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,182,610	0	C	G	6.1%	262.7 / 5.8	82	R18G (CGG→GGG)	PP_4562	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (41/36);  new base G (2/3);  total (43/39)
Fisher's exact test for biased strand distribution p-value = 6.65e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.75e-01

GAGGCCGGCACTGGCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCT  >  NC_002947/5182463‑5182758                                                                                                                                                    |                                                                                                                                                     gggacgtgCACGGGCACCGCAATTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCccgg                                                                                                                                                    >  2:1090342/8‑150 (MQ=255)     ccGGCACTTGCACCTGAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACACTGCCAAGCAGCTGAGCGTCACCGCCCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCccggccgg                                                                                                                                                <  1:101189/150‑1 (MQ=255)        gCACTGGCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGc                                                                                                                                             <  2:1211113/150‑1 (MQ=255)        gCACTGGCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACCGGCCGGGCCAGTGc                                                                                                                                             <  2:834660/150‑1 (MQ=255)          aCTGGCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTAGAAAAAACCGCCGCCCTCACCCAGCTCAGCGAAGCACTTGCCCGGCCGGTGCTg                                                                                                                                           >  2:1084226/1‑150 (MQ=255)           cTGGCACCGCAAGTCACTAACCTCGCTACGCCCCGTGGCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGa                                                                                                                                          <  2:262000/150‑1 (MQ=255)           cTGGCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGa                                                                                                                                          >  1:843544/1‑150 (MQ=255)              gCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACca                                                                                                                                       >  1:925646/1‑150 (MQ=255)              gCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACca                                                                                                                                       >  2:364261/1‑150 (MQ=255)                 ccGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACcagca                                                                                                                                    <  2:6947/150‑1 (MQ=255)                   gCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCATGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCagcagc                                                                                                                                  >  2:618632/1‑150 (MQ=255)                    cAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAgcagca                                                                                                                                 >  1:827057/1‑150 (MQ=255)                          acaAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACca                                                                                                                                       <  2:1321248/138‑1 (MQ=255)                          acaAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACca                                                                                                                                       >  1:1321248/1‑138 (MQ=255)                            aaaCCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTCTCCATCGGCAACACCCTg                                                                                                                         >  2:82358/1‑150 (MQ=255)                            aaaCCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCAGCGCTGACCAGCCGCAACAccccg                                                                                                                         >  2:82227/1‑148 (MQ=255)                            aaaCCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACCGGCCGGGCCAGTGCttcgtt                                                                                                                                       >  1:277549/1‑131 (MQ=255)                            aaaCCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACCGGCCGGGCCAGTGCttcgtt                                                                                                                                       <  2:277549/136‑6 (MQ=255)                                cTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCg                                                                                                                     >  1:858965/1‑150 (MQ=255)                                     tCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATc                                                                                                                >  2:298256/1‑150 (MQ=255)                                     tCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATc                                                                                                                >  2:585940/1‑150 (MQ=255)                                     tCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGCCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATc                                                                                                                >  2:3556/1‑150 (MQ=255)                                       cGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCt                                                                                                              <  1:364261/150‑1 (MQ=255)                                        gCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAGCGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTg                                                                                                             <  1:1084226/150‑1 (MQ=255)                                          cccGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGa                                                                                                           <  2:1401627/150‑1 (MQ=255)                                                ccACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGAcaacaa                                                                                                     >  2:956063/1‑150 (MQ=255)                                                ccACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGAcaacaa                                                                                                     >  2:1308319/1‑150 (MQ=255)                                                  aCCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAAc                                                                                                   <  2:925646/150‑1 (MQ=255)                                                    cccTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCGTCCTCAGCAACAAACGc                                                                                                 >  2:747956/1‑150 (MQ=255)                                                          caATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCtg                                                                                           >  2:868461/1‑150 (MQ=255)                                                           aaTGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCtgt                                                                                          <  1:1090342/150‑1 (MQ=255)                                                              gCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTgtggt                                                                                       >  2:790026/1‑150 (MQ=255)                                                                  gcattCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCACCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGtt                                                                                   >  2:281325/4‑150 (MQ=255)                                                                  gCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGtt                                                                                   >  1:787999/1‑150 (MQ=255)                                                                   ccTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGc                                                                                                                                             <  1:931247/91‑1 (MQ=255)                                                                   ccTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGc                                                                                                                                             >  2:931247/1‑91 (MQ=255)                                                                   ccTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACCGGCCGGGCCAGTGc                                                                                                                                             >  2:293656/1‑91 (MQ=255)                                                                   ccTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACCGGCCGGGCCAGTGc                                                                                                                                             <  1:293656/91‑1 (MQ=255)                                                                      tCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGaa                                                                               <  1:298256/150‑1 (MQ=255)                                                                        aTAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGt                                                                             >  2:756163/1‑150 (MQ=255)                                                                             gTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCgg                                                                        >  1:1364582/1‑150 (MQ=255)                                                                             gTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCgg                                                                        <  1:618632/150‑1 (MQ=255)                                                                                       aGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCNGNANAACAAACGCAACCTGNGGTGGTTCGAAGTGCCGGTGGCActgct                                                              <  1:756163/150‑1 (MQ=255)                                                                                         gAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTgg                                                            >  1:1336797/1‑150 (MQ=255)                                                                                          aaCCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTggg                                                           <  1:868461/150‑1 (MQ=255)                                                                                             cTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCaa                                                        >  1:147031/1‑150 (MQ=255)                                                                                             cTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCaa                                                        >  2:635966/1‑150 (MQ=255)                                                                                               gCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAgg                                                      <  1:956063/150‑1 (MQ=255)                                                                                                cATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAggg                                                     >  1:6744/1‑150 (MQ=255)                                                                                                   gAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACACACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCa                                                  >  2:900058/1‑150 (MQ=255)                                                                                                   gAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGaa                                                                               <  2:800185/121‑1 (MQ=255)                                                                                                   gAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGaa                                                                               >  1:800185/1‑121 (MQ=255)                                                                                                    aaaCTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAg                                                 <  2:6744/150‑1 (MQ=255)                                                                                                    aaaCTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAg                                                 <  2:827057/150‑1 (MQ=255)                                                                                                     aaCTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTg                                                                                                                                              >  2:636384/1‑56 (MQ=255)                                                                                                     aaCTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTg                                                                                                                                              <  1:636384/56‑1 (MQ=255)                                                                                                     aaCTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAgg                                                >  1:257771/1‑150 (MQ=255)                                                                                                       cTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGcc                                              <  2:858965/150‑1 (MQ=255)                                                                                                             aaaCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTgg                                        <  1:1308319/150‑1 (MQ=255)                                                                                                             aaaCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTgg                                        >  2:619015/1‑150 (MQ=255)                                                                                                              aaCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTggg                                       <  1:619015/150‑1 (MQ=255)                                                                                                              aaCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTggg                                       <  2:11372/150‑1 (MQ=255)                                                                                                               aCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACaaa                                                                                                    >  1:1313876/1‑88 (MQ=255)                                                                                                               aCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACaaa                                                                                                    <  2:1313876/88‑1 (MQ=255)                                                                                                               aCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAgg                                                      >  2:378007/1‑134 (MQ=255)                                                                                                               aCCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAgg                                                      <  1:378007/134‑1 (MQ=255)                                                                                                               aCCCGCCGCCAGCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGTACAACAAACGCAACCTGTGGTGGTCCGAAGTGCCGGTGGCACTGCTGGGTAAGGGCCAGGCCGACTGGGt                                      <  2:341466/150‑1 (MQ=255)                                                                                                                cccGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTc                                     >  2:220965/1‑150 (MQ=255)                                                                                                                cccGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTc                                     <  2:787999/150‑1 (MQ=255)                                                                                                                  cgccgcCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGCCAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCaa                                   >  1:92953/1‑150 (MQ=255)                                                                                                                   gccgccATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAAc                                  <  1:82358/150‑1 (MQ=255)                                                                                                                   gccgccATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAAc                                  <  1:82227/150‑1 (MQ=255)                                                                                                                        caccGCCCGGCGCAACGAAGCACTGGCCCGGCCTGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTCGACAACAAACGCAACCCCTGGAGGATAGAACTGCCGATGGCACTGCTGGGCAAGAGCCACGCCGACCGTGTCGACctgct                             >  2:651109/4‑150 (MQ=255)                                                                                                                        cATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACctgct                             <  2:510528/150‑1 (MQ=255)                                                                                                                                     aaCGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCg                >  2:1095268/1‑150 (MQ=255)                                                                                                                                        gAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGaa             >  2:1400361/1‑150 (MQ=255)                                                                                                                                         aaGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCActgctg                                                             >  1:758627/1‑101 (MQ=255)                                                                                                                                         aaGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCActgctg                                                             <  2:758627/101‑1 (MQ=255)                                                                                                                                              cTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGAt       <  1:363973/150‑1 (MQ=255)                                                                                                                                              cTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGAt       <  1:3556/150‑1 (MQ=255)                                                                                                                                              cTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGAt       <  2:909635/150‑1 (MQ=255)                                                                                                                                              cTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGAt       <  1:790026/150‑1 (MQ=255)                                                                                                                                               tGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATg      >  1:639010/1‑150 (MQ=255)                                                                                                                                               tGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATg      <  2:1336797/150‑1 (MQ=255)                                                                                                                                                  cccGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGc   <  1:220965/150‑1 (MQ=255)                                                                                                                                                   ccggccggTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCTAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCt  <  1:1095268/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                     GAGGCCGGCACTGGCACCGCAAGTCACAAACCTTGCTCCGCCCCGTGCCACCCCTCACAATGCAGGCCTTCATAGCGTCAACGCGCAGGAACCTGCATGAAACTCGACAAACCCGCCGCCATCGCCCGGCGCAACGAAGCACTGGCCCGGCCGGTGCTGACCAGCAGCAACACCCTGTTCGCCATCCTGGACAACAAACGCAACCTGTGGTGGTTCGAAGTGCCGGTGGCACTGCTGGGCAAGGGCCAGGCCGACTGGGTCAACCTGCTGCTGCACACGCCGGAAAGCGATGAGCT  >  NC_002947/5182463‑5182758

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑