Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,816,589	G→A	6.0%	D115N (GAC→AAC)	PP_5092 →	NLP/P60 protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,816,589	0	G	A	6.0%	221.3 / 2.7	76	D115N (GAC→AAC)	PP_5092	NLP/P60 protein
Reads supporting (aligned to +/- strand):  ref base G (35/36);  new base A (3/2);  total (38/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGACCG  >  NC_002947/5816446‑5816737                                                                                                                                                |                                                                                                                                                      ccgccATGAAGCTTTATCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCg                                                                                                                                                 >  2:176179/1‑150 (MQ=255)   gccATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGAc                                                                                                                                               <  2:1477776/150‑1 (MQ=255)      aTGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCAACGTCGACCTg                                                                                                                                            >  2:254712/1‑150 (MQ=255)      aTGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCAACGTCGACCTg                                                                                                                                            >  2:257744/1‑150 (MQ=255)           gCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCGTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCcgcg                                                                                                                                       >  1:613664/1‑150 (MQ=255)                       cagcGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGc                                                                                                                           >  1:639743/1‑150 (MQ=255)                        agcGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAATGCTCTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCTACCTGCCGCGCACCTCCAATGc‑c                                                                                                                         <  2:945049/150‑1 (MQ=255)                           gATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CCTg                                                                                                                      >  1:570774/1‑150 (MQ=255)                           gATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATg                                                                                                                      <  2:654298/150‑1 (MQ=255)                            aTGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCAACGTCGACCTGCCGCGCACCTCCAATGC‑CATgg                                                                                                                     <  1:257744/150‑1 (MQ=255)                            aTGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCAACGTCGACCTGCCGCGCACCTCCAATGC‑CATgg                                                                                                                     <  1:254712/150‑1 (MQ=255)                             tgtgCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGc                                                                                                                    <  2:420443/150‑1 (MQ=255)                                   gAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCgggg                                                                                                              <  2:931485/150‑1 (MQ=255)                                         tGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCAGGGGCCACgg                                                                                                        >  1:904551/1‑150 (MQ=255)                                          gCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACggg                                                                                                       <  1:1737886/150‑1 (MQ=255)                                          gCGGTAAATGTGCTCGGTACGCCTTCTGTGGGGGGCGGCAGACGCCCGAAAAAAGGCTTACACTTCCGATCTCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACggg                                                                                                       <  1:512350/150‑1 (MQ=255)                                           cGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACgggg                                                                                                      >  2:1136046/1‑150 (MQ=255)                                            ggTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGTCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGt                                                                                                     <  1:1567007/150‑1 (MQ=255)                                            ggTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACgggg                                                                                                      >  2:1149722/1‑149 (MQ=255)                                            ggTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACgggg                                                                                                      <  1:1149722/149‑1 (MQ=255)                                             gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTg                                                                                                    >  1:793859/1‑150 (MQ=255)                                             gTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTAGAAAAAATTGGCCGTCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCCCGGGGTg                                                                                                    >  2:49884/1‑150 (MQ=255)                                                 aTGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGGGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAgg                                                                                                <  2:93122/150‑1 (MQ=255)                                                 aTGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCT‑TTATGCACATGGCCCGGGGCCACGGGGTGAAgg                                                                                                <  1:49884/150‑1 (MQ=255)                                                   gtgCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTg                                                                                              <  2:793859/150‑1 (MQ=255)                                                       tCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCa                                                                                          >  1:1411741/1‑150 (MQ=255)                                                        cGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCaa                                                                                         <  1:1136046/150‑1 (MQ=255)                                                         ggTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAg                                                                                        >  2:424274/1‑150 (MQ=255)                                                          gTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACggg                                                                                                       <  2:1146427/134‑1 (MQ=255)                                                          gTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACggg                                                                                                       >  1:1146427/1‑134 (MQ=255)                                                          gTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAgg                                                                                       <  1:673342/150‑1 (MQ=255)                                                           tACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAggg                                                                                      >  1:483258/1‑150 (MQ=255)                                                              tCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGACGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACggg                                                                                                       >  1:53664/1‑130 (MQ=255)                                                              tCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGACGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACggg                                                                                                       <  2:53664/130‑1 (MQ=255)                                                               ccTTAGGTTTGGGGCGGAAGAAGACCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGAt                                                                                  <  2:743094/150‑1 (MQ=255)                                                                gctATGTTGGGCGCGGCAGCAGCCCGAAAAAAGGATTTGACTGCAGTGGGCTGGTCGAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATc                                                                                 <  2:1755401/148‑1 (MQ=255)                                                                   aTGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTgg                                                                                             <  2:237882/135‑1 (MQ=255)                                                                   aTGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTgg                                                                                             >  1:237882/1‑135 (MQ=255)                                                                      tttGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGc                                                                           >  2:916475/1‑150 (MQ=255)                                                                                cagcagCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCt                                                                 <  2:127197/150‑1 (MQ=255)                                                                                  gcagcCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGa                                                               <  1:916475/150‑1 (MQ=255)                                                                                     gcCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCAGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtct                                                            <  2:904551/150‑1 (MQ=255)                                                                                      cccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGAtctt                                                           >  2:1727871/1‑150 (MQ=255)                                                                                       ccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCgcgc                                                                                                                                      <  2:87587/75‑1 (MQ=255)                                                                                       ccGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCgcgc                                                                                                                                      >  1:87587/1‑75 (MQ=255)                                                                                         gAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                        <  2:570774/150‑1 (MQ=255)                                                                                         gAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                        <  1:1727871/150‑1 (MQ=255)                                                                                          aaaaaaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCttcttc                                                       <  2:613664/150‑1 (MQ=255)                                                                                           aaaaaGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGCTGGCAAAGTGTTATGTGAATCCTTGCGACCTAATCTTCTTCa                                                      >  2:310553/1‑150 (MQ=255)                                                                                               aGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                      >  2:1626168/1‑130 (MQ=255)                                                                                               aGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGc                                                                      <  1:1626168/130‑1 (MQ=255)                                                                                                ggCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                        <  2:1318707/143‑1 (MQ=255)                                                                                                ggCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATcttctt                                                        >  1:1318707/1‑143 (MQ=255)                                                                                                      gACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGc                                           <  2:1083485/150‑1 (MQ=255)                                                                                                        cTGCAGCGGGCTGGTCAAATACGCCTTCAACGACAGGGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCg                                         >  1:1893253/1‑150 (MQ=255)                                                                                                         tGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACa                                                   <  1:714074/139‑1 (MQ=255)                                                                                                         tGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACa                                                   >  2:714074/1‑139 (MQ=255)                                                                                                           cAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCg                                      >  1:1422110/1‑150 (MQ=255)                                                                                                            aGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCgg                                     <  2:1411741/150‑1 (MQ=255)                                                                                                                    ggTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCa                             >  1:665441/1‑150 (MQ=255)                                                                                                                     gTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCAc                            >  1:1611396/1‑150 (MQ=255)                                                                                                                        aaaTACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTa                         >  1:492374/1‑150 (MQ=255)                                                                                                                           tACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAag                                              >  1:1507236/1‑126 (MQ=255)                                                                                                                           tACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACACCAag                                              <  2:1507236/126‑1 (MQ=255)                                                                                                                                cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTa                 >  1:277059/1‑150 (MQ=255)                                                                                                                                cTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTa                 >  2:1576883/1‑150 (MQ=255)                                                                                                                                 ttCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc                >  1:1163190/1‑150 (MQ=255)                                                                                                                                 ttCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc                >  1:633132/1‑150 (MQ=255)                                                                                                                                 ttCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAc                >  2:307153/1‑150 (MQ=255)                                                                                                                                  tCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTAcc               <  2:1422110/150‑1 (MQ=255)                                                                                                                                   ggtcgacgTCGGCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCt              >  1:103943/4‑150 (MQ=255)                                                                                                                                     acgacgTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTTATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTgg            <  1:292917/150‑1 (MQ=255)                                                                                                                                             gCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGAc    <  1:1824464/150‑1 (MQ=255)                                                                                                                                              ccGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAg                                                                                                 >  1:943829/1‑56 (MQ=255)                                                                                                                                              ccGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAg                                                                                                 <  2:943829/56‑1 (MQ=255)                                                                                                                                               cGACGTCGACCTGCCTCGCACCCCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGACCg  >  2:72819/1‑150 (MQ=255)                                                                                                                                                |                                                                                                                                                      CCGCCATGTCGCCTAAGCAAAGCAGCGATGTGCTGAGCCGTGCGGTAAATGTGCTCGGTACTCCTTATGTTTGGGGCGGCAGCAGCCCGAAAAAAGGCTTTGACTGCAGCGGGCTGGTCAAATACGCCTTCAACGACGTCGCCGACGTCGACCTGCCGCGCACCTCCAATGC‑CATGGCCCGGGGCCACGGGGTGAAGGTGGCCAAGGGTGATCTGAAGCCTGGCGACCTGATCTTCTTCAACATCAAGAGCCGTCGGGTCAATCACGTAGCCATCTACCTGGGTAACGACCG  >  NC_002947/5816446‑5816737

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