Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,652,968	A→T	7.0%	intergenic (‑69/+92)	PP_4961 ← / ← PP_4962	lipoprotein/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,652,968	0	A	T	7.0%	157.7 / 3.4	67	intergenic (‑69/+92)	PP_4961/PP_4962	lipoprotein/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/27);  new base T (2/3);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.86e-01

CCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCC  >  NC_002947/5652821‑5653113                                                                                                                                                     |                                                                                                                                                  ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCATGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                                   <  3:31644/148‑1 (MQ=255) ccGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCa                                                                                                                                                   <  7:85539/148‑1 (MQ=255)     ggCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGtctc                                                                                                                                              <  4:8116/149‑1 (MQ=255)         ggTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGCACCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCtgt                                                                                                                                           <  1:278301/149‑1 (MQ=255)          gTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                         <  7:70832/149‑1 (MQ=255)               gcCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCg                                                                                                                                    <  8:112879/149‑1 (MQ=25)                        agcagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGccc                                                                                                                           >  8:39161/1‑149 (MQ=255)                         gcagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCcg                                                                                                                          >  5:199905/1‑149 (MQ=255)                          cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                                                             >  4:73487/1‑145 (MQ=255)                          cagTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGc                                                                                                                             <  3:73487/145‑1 (MQ=255)                           agTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                                                 >  4:95663/1‑124 (MQ=255)                           agTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGt                                                                                                                                                 <  3:95663/124‑1 (MQ=255)                                gCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                         >  8:50733/1‑127 (MQ=255)                                gCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTgg                                                                                                                                         <  7:50733/127‑1 (MQ=255)                                gCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAggg                                                                                                                   >  7:223909/1‑149 (MQ=255)                                 cATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGc                                                                                                                  >  8:142272/1‑149 (MQ=255)                                     cTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGca                                                                                                              >  6:288572/1‑149 (MQ=255)                                      tGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGca                                                                                                              >  5:307628/1‑148 (MQ=255)                                         cAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCGGTGCGGC‑CCTTCGCGGGCTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGAGAAGGGCTGCACAGc                                                                                                          <  6:4517/149‑1 (MQ=14)                                          aGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGGGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCg                                                                                                         <  3:95436/149‑1 (MQ=255)                                              tCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGccc                                                                                                     >  7:367341/1‑149 (MQ=255)                                               ccccAGGGCCGCCACGGCAAAACCCCGTTCCATTCCTTTTCCTCCATTGGGGCCGGCGGCGCCGC‑CCTTCCCGGCCGCCCCCGCCCCCCCCGGTCCAGTCCCGGCTCTGTGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGcccc                                                                                                    <  8:342102/147‑1 (MQ=255)                                                gCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCct                                                                                                   <  6:307628/149‑1 (MQ=255)                                                 ccAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCctt                                                                                                  >  5:293744/1‑149 (MQ=255)                                                  cAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACCGTGCAGGCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCcttc                                                                                                 <  5:146477/149‑1 (MQ=255)                                                       ccGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCaa                                                                                            <  6:293744/149‑1 (MQ=255)                                                         tacccGGCCAAACGGGGTTTCATTCCGTTCCCTCGTTTGGGGGCTGCGGGGCTGC‑CGTTCTCGGGCTGGCCCGGTCCCACTGTAACTTTCCCGTCGGTGTGGGAGGGGGCGTGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGa                                                                                          <  2:8702/145‑1 (MQ=255)                                                                   aaaGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTGGTCAAGACCACAAATc                                                                                 >  4:258592/1‑148 (MQ=255)                                                                   aaaGCGATTTAATACAGTTTCCTCGATTGGGGGCTGGTGTCCAGC‑CCTTCTCGGGCGCCCCCGCACCCACTGGTACCGGCCAGTCTCTGTGGTAGCGGGCTAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCa                                                                                <  7:142272/149‑1 (MQ=25)                                                                       cgcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGt                                                                            <  5:288572/149‑1 (MQ=255)                                                                        gcTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGtt                                                                           <  6:199905/149‑1 (MQ=255)                                                                            tcattcaGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAgg                                                                       >  3:306604/1‑149 (MQ=255)                                                                             cattcaGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAACGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGt                                                                      >  6:379772/1‑149 (MQ=255)                                                                                      tCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTAttt                                                             >  3:187654/1‑148 (MQ=255)                                                                                       ccTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCGTCTTCAAGACGACAAATCAACGTTCAGGTTGCCTATTAc                                                            >  5:228518/1‑149 (MQ=255)                                                                                         tCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTAccc                                                          >  3:7547/1‑149 (MQ=255)                                                                                             tttGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGAt                                                       >  2:1299/1‑148 (MQ=255)                                                                                             tttGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGGGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATc                                                      <  8:220010/149‑1 (MQ=255)                                                                                                ggggCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGGCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                                   <  8:367341/149‑1 (MQ=255)                                                                                                 gggCTGCTGTGCAGC‑CTTCCTCGGGCGCGCCCGCTCCCACAGGAACAGGCCAGTCTCTGTGGTTGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                                   <  7:49821/148‑1 (MQ=255)                                                                                                 gggCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACCGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAgg                                                  <  4:256457/149‑1 (MQ=255)                                                                                                 gggCTGCTGTGCAGC‑CCGTGTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                                   <  5:297216/148‑1 (MQ=255)                                                                                                 ggcctggtgTGCTGC‑CCTTGGGGGGCCGGCCGGCCCCCACCGTTCCTGTGCAGTCTCGGGGGCTGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAAGGTTCAGGTTGCCTATTACCCGATCAAgg                                                  <  4:262536/146‑1 (MQ=255)                                                                                                  gcctgcggGGCCGC‑CCTTCTCGGGCGCCCCCGCCCCCACCGCTACAGTGCCGCCTCTGGGGGCGGGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAggg                                                 <  6:228518/147‑1 (MQ=255)                                                                                                    ctgctgTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGc                                                >  3:41351/1‑148 (MQ=255)                                                                                                        tgtgCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                                   >  4:348117/1‑141 (MQ=255)                                                                                                        tgtgCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAg                                                   <  3:348117/141‑1 (MQ=255)                                                                                                          tgctgc‑CCTTCTCGGGCTCGCCCGCTCCCACAGTTACTGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCgg                                         <  5:322404/149‑1 (MQ=255)                                                                                                           gCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgc                                        >  8:123930/1‑149 (MQ=255)                                                                                                            cAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGgcg                                       >  1:244203/1‑149 (MQ=255)                                                                                                             aGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGACGGGCCGCCCAGCGGCCCCGTCTTCAAGACCACAAAGCAACGTTCAGGTTGGCGATTTCCCGCTCAAGGGCAACCCGgcgc                                      >  5:233400/1‑149 (MQ=38)                                                                                                             aGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGAGCAAGGGCAACCCGgcgc                                      >  2:334216/1‑149 (MQ=255)                                                                                                                  cTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGAAGTCTCTGTGGGAGCGGGCGAGCACGCGCAGGGCCGCAAAGCGGCCCCGTATTCAAGACCACAAAGCAACGTTCAGGTTGCCTATTTCGCGATCAAGGGCAACCCGGCGCCGTg                                  >  4:337751/1‑149 (MQ=255)                                                                                                                   ttCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTgc                                 >  5:290206/1‑149 (MQ=255)                                                                                                                         ggCGCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCTAGGGCAACCCGGCGCCTTGCGCtggt                           <  1:11631/149‑1 (MQ=255)                                                                                                                         ggCGCGCCCGCTCCCACAGGAACAGTGCAGTCTCTGGGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggt                           <  8:223909/149‑1 (MQ=255)                                                                                                                          gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACCGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTTTTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtg                          >  3:51105/1‑149 (MQ=255)                                                                                                                          gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                                   <  3:225340/140‑1 (MQ=255)                                                                                                                          gCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCtt                                   >  4:225340/1‑140 (MQ=255)                                                                                                                           cTCGCCCGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCtggtgg                         >  4:13148/1‑149 (MQ=25)                                                                                                                              gCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCACGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCt                      >  1:281917/1‑149 (MQ=255)                                                                                                                               cccGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCtt                     >  5:245833/1‑149 (MQ=255)                                                                                                                               cccGCTCCCACAGAGACTGCACTGTACCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCtt                     >  4:386470/1‑149 (MQ=25)                                                                                                                                ccGCTCCCTCTGCTACAGTGCTGGGTCTGTGGGTGGGGGCGAGCCCGCGAAGGGCCGCACAGCTGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtc                    <  4:51105/149‑1 (MQ=255)                                                                                                                                      ccACCGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTa              <  4:41351/149‑1 (MQ=255)                                                                                                                                      ccACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTa              >  4:246419/1‑149 (MQ=255)                                                                                                                                            gTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                   >  4:311665/1‑138 (MQ=255)                                                                                                                                            ctACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTtcg                   <  3:311665/137‑1 (MQ=255)                                                                                                                                                  tGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGcc  <  1:1299/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                  CCGGGGCCGGTTGCGCCGGCTGCAGCAGTGAGCATCCTGCCAGTGTCGCCAGGGCCGTCACGGCCAAAAGCGCTTTCATTCAGTTTCCTCGATTTGGGGCTGCTGTGCAGC‑CCTTCTCGGGCTCGCCCGCTCCCACAGGTACAGTGCAGTCTCTGTGGGAGCGGGCGAGCCCGCGAAGGGCCGCACAGCGGCCCCTTCTTCAAGACCACAAATCAACGTTCAGGTTGCCTATTACCCGATCAAGGGCAACCCGGCGCCTTGCGCTGGTGGTCTTCGTCGTACTTTTCCAGGCC  >  NC_002947/5652821‑5653113

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