Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	990,700	T→A	10.3%	L193Q (CTG→CAG)	PP_0856 →	lipoprotein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	990,700	0	T	A	10.3%	278.4 / 21.3	98	L193Q (CTG→CAG)	PP_0856	lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (41/46);  new base A (5/5);  total (46/52)
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.90e-01

AACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCTGGTG  >  NC_002947/990552‑990849                                                                                                                                                     |                                                                                                                                                       aacaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTg                                                                                                                                                       <  1:110432/150‑1 (MQ=255) aacaacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTg                                                                                                                                                       >  1:701732/1‑150 (MQ=255)    aacGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCg                                                                                                                                                    <  1:1030430/150‑1 (MQ=255)     acGGTGACGTGGTGGTGGTGCAGACACAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCgg                                                                                                                                                   <  2:1362095/150‑1 (MQ=255)       ggTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCTTTGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGtt                                                                                                                                                 >  1:990199/1‑150 (MQ=255)         tGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGc                                                                                                                                               <  1:1428414/150‑1 (MQ=255)          gACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGcc                                                                                                                                              >  1:1368884/1‑150 (MQ=255)             gtggtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtg                                                                                                                                             <  1:1524577/148‑4 (MQ=255)             gtggtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtg                                                                                                                                             >  2:1524577/1‑145 (MQ=255)              tggtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGcc                                                                                                                                          >  1:1218682/1‑150 (MQ=255)               ggtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCAGGCCt                                                                                                                                         >  2:959524/1‑150 (MQ=255)                gtggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTg                                                                                                                                        >  2:1659612/1‑150 (MQ=255)                 tggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑t                                                                                                                                      <  2:1461904/150‑1 (MQ=255)                 tggtggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑t                                                                                                                                      >  2:296537/1‑150 (MQ=255)                     ggtgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTAc                                                                                                                                  <  2:43007/150‑1 (MQ=255)                       tgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCg                                                                                                                                <  2:1794986/150‑1 (MQ=255)                       tgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCg                                                                                                                                <  1:461238/150‑1 (MQ=255)                       tgCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCg                                                                                                                                <  1:1777359/150‑1 (MQ=255)                        gCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCgg                                                                                                                               >  1:736766/1‑150 (MQ=255)                           gACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCaa                                                                                                                            <  1:187965/150‑1 (MQ=255)                              ccAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGt                                                                                                                         >  1:208329/1‑150 (MQ=255)                               cAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTa                                                                                                                        >  2:1495666/1‑150 (MQ=255)                                aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAg                                                                                                                       >  2:1435831/1‑150 (MQ=255)                                aGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAg                                                                                                                       >  2:1338778/1‑150 (MQ=255)                                          ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGGGGCTTCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGa                                                                                                             >  1:593537/1‑150 (MQ=255)                                          ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGa                                                                                                             <  2:1939680/150‑1 (MQ=255)                                          ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCAGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGa                                                                                                             <  1:959524/150‑1 (MQ=255)                                          ccTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCATGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGa                                                                                                             >  2:344856/1‑150 (MQ=255)                                           cTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGAc                                                                                                            >  1:813384/1‑150 (MQ=255)                                           cTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGAc                                                                                                            <  2:703740/150‑1 (MQ=255)                                           cTGATCGGCCTGGATACCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGAc                                                                                                            <  2:815268/150‑1 (MQ=255)                                            tGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGt                                                                                                                                                  >  2:1482731/1‑112 (MQ=255)                                            tGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGt                                                                                                                                                  <  1:1482731/112‑1 (MQ=255)                                              aTCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAg                                                                                                                           >  1:1451327/1‑132 (MQ=255)                                              aTCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAg                                                                                                                           <  2:1451327/132‑1 (MQ=255)                                                  gCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtgg                                                                                                                                            >  2:1706368/1‑108 (MQ=255)                                                  gCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCAGGCGGTTggtgg                                                                                                                                            <  1:1706368/112‑5 (MQ=255)                                                   ccTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCAACACCGGGCCGATCGCCACCCACAGCATGGCGGCAG‑TGGCCTGCTCTACCGGCGAGGTAGTGGCGGTGGACATcaacaa                                                                                                    >  1:1399163/1‑150 (MQ=255)                                                    cTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCaacaac                                                                                                   >  1:1045258/1‑150 (MQ=255)                                                    cTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCaacaac                                                                                                   <  2:1218682/150‑1 (MQ=255)                                                     tGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACg                                                                                                  >  2:1643792/1‑150 (MQ=255)                                                        aTGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGgtg                                                                                               <  1:1941333/150‑1 (MQ=255)                                                           ccgccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGc                                                                                            <  1:1495666/150‑1 (MQ=255)                                                           ccgccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGc                                                                                            >  2:30541/1‑150 (MQ=255)                                                            cgccACTGGCGGCCGTCGTTGGCGCTCCGAAAATTCGCCGGCCGTGCTTACGCTAGGTGGCGCCGGTGCGCCGATCTCCCCCAACCGCCTGCCGGTTGGCGGCCTG‑TCTACCCGCAGGGTAGTGGCGGTGGACATCAACAACGGAGTGcc                                                                                           >  2:1730167/1‑150 (MQ=255)                                                             gccACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCg                                                                                          <  2:208329/150‑1 (MQ=255)                                                               cACTGGTGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGgt                                                                                        >  2:1474434/1‑150 (MQ=255)                                                               cACTGGTGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGgt                                                                                        <  1:1474434/150‑1 (MQ=255)                                                                   ggCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTgg                                                                                    <  1:783611/150‑1 (MQ=255)                                                                   ggCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTgg                                                                                    <  1:296537/150‑1 (MQ=255)                                                                   ggCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTgg                                                                                    <  1:1643792/150‑1 (MQ=255)                                                                    gCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTggg                                                                                   >  1:931315/1‑150 (MQ=255)                                                                       aCCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCTCCACCAACCGTCTGGCTGTTGCCTGCATG‑TCTACCGGCAAGGTAGTGGCAGTCGACATCAACCACGGTGTGTAGgtgttgcata                                                                                <  2:517981/150‑7 (MQ=255)                                                                       aCCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGgaga                                                                                <  1:1435831/150‑1 (MQ=255)                                                                       aCCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGgaga                                                                                <  2:1045258/150‑1 (MQ=255)                                                                       aCCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGgaga                                                                                <  2:931315/150‑1 (MQ=255)                                                                              ttGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGt                                                                         <  1:1532251/150‑1 (MQ=255)                                                                              ttGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGt                                                                         <  1:1739953/150‑1 (MQ=255)                                                                              ttGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCACCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGt                                                                         >  1:464487/1‑150 (MQ=255)                                                                                  aTCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGcc                                                                     >  1:1055857/1‑150 (MQ=255)                                                                                            aCTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAgg                                                           >  1:1611832/1‑150 (MQ=255)                                                                                              tCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGc                                                         >  2:961270/1‑150 (MQ=255)                                                                                               cGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGcc                                                        <  2:1368884/150‑1 (MQ=255)                                                                                                  cggccgTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGtt                                                     <  1:961270/150‑1 (MQ=255)                                                                                                      cgTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGa                                                 <  1:1338778/150‑1 (MQ=255)                                                                                                         tcggACGCTGCGGGGCACCTGTGCACCGATGGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAACGTACTGGCGGTGGACATCAGGAGCGGTGAGCCGGTGTGGGAGGGCCCCGTGGCCATCCCGAAGGGCCGTTCCGAGCt                                              <  2:813384/147‑1 (MQ=255)                                                                                                         gCTGACGTTGCGTGGCACCGGGGCGCCGATCGCCACATTCCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGCGCCGCGTGGCCATCCCGAAGGGCCGTTCTGAGCt                                              <  1:593887/150‑1 (MQ=255)                                                                                                          cTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTg                                             <  1:1582061/150‑1 (MQ=255)                                                                                                                    gTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTgg                                   >  2:924989/1‑150 (MQ=255)                                                                                                                     tGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGt                                  >  1:39839/1‑150 (MQ=255)                                                                                                                        cACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGa                               >  1:710728/1‑150 (MQ=255)                                                                                                                        cACCGGTGCCCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGCGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGACGGGCCGTTCCGAGCTGGACCGTGTGGTTGa                               <  2:1089564/150‑1 (MQ=255)                                                                                                                         aCCGGTGCGCCGATCTCCACCAACCGCCTGGCGGTTGCCGTCCAG‑TCTACCGGCAAGGTAGTGCCGATGAACATCAACAACTGTGTGCCGGCGCCGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCGGGGCCGGACCGGGTGGGTGat                              >  2:284118/1‑150 (MQ=255)                                                                                                                         aCCGGTGCGCCGATCGCCACCACCCGCCGGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGat                              <  2:593537/150‑1 (MQ=255)                                                                                                                         aCCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGat                              >  2:1570672/1‑150 (MQ=255)                                                                                                                            ggTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATc                           >  1:1828879/1‑150 (MQ=255)                                                                                                                             gTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGcc                                                        <  1:484920/120‑1 (MQ=255)                                                                                                                             gTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGcc                                                        >  2:484920/1‑120 (MQ=255)                                                                                                                             gTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCg                          >  1:1509218/1‑150 (MQ=255)                                                                                                                                   cGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCCAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAcggcg                    <  2:1399163/150‑1 (MQ=255)                                                                                                                                        gCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAc                        >  1:16561/1‑141 (MQ=255)                                                                                                                                        gCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAc                        <  2:16561/141‑1 (MQ=255)                                                                                                                                        gCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAc                        >  1:16470/1‑141 (MQ=255)                                                                                                                                        gCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGAc                        <  2:16470/141‑1 (MQ=255)                                                                                                                                         ccggcaACCGCCAGGCGGTTGGTGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCt                                              <  1:227307/114‑1 (MQ=255)                                                                                                                                         ccaccaACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGgt                                                                                                >  1:1481038/1‑68 (MQ=255)                                                                                                                                         ccaccaACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGgt                                                                                                <  2:1481038/68‑1 (MQ=255)                                                                                                                                         ccaccaACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgt              <  2:736766/150‑1 (MQ=255)                                                                                                                                         caggcaACCGCCAGGCGGTTGGTGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCt                                              >  2:227307/5‑118 (MQ=255)                                                                                                                                          cggcaACCGCCAGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGaa                                                             >  2:1082777/4‑102 (MQ=255)                                                                                                                                          cggcaACCGCCAGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGaa                                                             <  1:1082777/99‑1 (MQ=255)                                                                                                                                          caccaACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCttgtt             >  1:1301868/1‑150 (MQ=255)                                                                                                                                             caACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCt          <  2:39839/150‑1 (MQ=255)                                                                                                                                               aCCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGt        >  2:1763607/1‑150 (MQ=255)                                                                                                                                                ccGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGgt                                                                                                >  2:750780/1‑61 (MQ=255)                                                                                                                                                ccGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGgt                                                                                                <  1:750780/61‑1 (MQ=255)                                                                                                                                                    cTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCtggt   <  1:30541/150‑1 (MQ=255)                                                                                                                                                    cTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCtggt   <  2:1884799/150‑1 (MQ=255)                                                                                                                                                     tGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCtggtg  >  2:1900127/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                       AACAACGGTGACGTGGTGGTGGTGCAGACCCAGGACGATCGCCTGATCGGCCTGGATGCCGCCACTGGCGACCGCCGTTGGATCTACGAAAACTCGCCGGCCGTGCTGACGCTGCGTGGCACCGGTGCGCCGATCGCCACCAACCGCCTGGCGGTTGCCGGCCTG‑TCTACCGGCAAGGTAGTGGCGGTGGACATCAACAACGGTGTGCCGGTGTGGGAGAGCCGCGTGGCCATCCCGAAGGGCCGTTCCGAGCTGGACCGTGTGGTTGATATCGACGGCGGCTTGTTGCTGTCTGGTG  >  NC_002947/990552‑990849

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