Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,949,113	T→C	10.8%	S141G (AGC→GGC)	PP_3485 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,949,113	0	T	C	10.8%	170.6 / 10.3	65	S141G (AGC→GGC)	PP_3485	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/29);  new base C (4/3);  total (33/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

CCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGACACT  >  NC_002947/3948967‑3949261                                                                                                                                                    |                                                                                                                                                      cccGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCt                                                                                                                                                       <  8:360499/147‑1 (MQ=255)     aTGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Ag                                                                                                                                                <  6:309150/149‑1 (MQ=255)       gAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTTTCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGTCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                                               <  8:192608/148‑1 (MQ=255)               cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACt                                                                                                                                      >  4:47830/1‑149 (MQ=255)               cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGCGAGGGc                                                                                                                                             <  1:156245/142‑1 (MQ=38)               cccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGCGAGGGc                                                                                                                                             >  2:156245/1‑142 (MQ=38)                ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                                               >  5:283367/1‑139 (MQ=255)                ccGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                                               <  6:283367/139‑1 (MQ=255)                   aGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                                               >  3:209519/1‑136 (MQ=255)                   aGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑Agg                                                                                                                                               <  4:209519/136‑1 (MQ=255)                     ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGCGAGGGc                                                                                                                                             >  3:197592/1‑136 (MQ=37)                     ccAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGCGAGGGc                                                                                                                                             <  4:197592/136‑1 (MQ=37)                       aaGGCCAGGCCGTGTTCGCGTATGACCATTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgag                                                                                                                                                >  1:233624/1‑128 (MQ=37)                       aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGgc                                                                                                                              >  8:12153/1‑149 (MQ=255)                       aaGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCT‑GCTCAGCGCAGTCGAC‑GCGCCGAGcgag                                                                                                                                                <  2:233624/131‑4 (MQ=38)                                     ttCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGcc                                                                                                                <  5:20423/149‑1 (MQ=255)                                     ttCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGcc                                                                                                                >  2:129290/1‑149 (MQ=255)                                        gcgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                                              <  3:47830/148‑1 (MQ=255)                                         cgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                            <  1:129290/149‑1 (MQ=255)                                         cgTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGGCCTGGTGGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                            <  7:340546/149‑1 (MQ=255)                                                  aaTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgt                                                                                                   <  7:12153/149‑1 (MQ=255)                                                   aTGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAgtg                                                                                                  <  3:289983/149‑1 (MQ=255)                                                     ggATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCGTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                                              <  6:297442/135‑1 (MQ=255)                                                     ggATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCgg                                                                                                              >  5:297442/1‑135 (MQ=255)                                                         tATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTgc                                                                                            >  5:357513/1‑149 (MQ=255)                                                           tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGc                                                                                                                                             <  6:204147/98‑1 (MQ=255)                                                           tCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGc                                                                                                                                             >  5:204147/1‑98 (MQ=255)                                                             aCTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGcc                                                                                         <  4:107650/148‑1 (MQ=255)                                                              cTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCt                                                                                       >  3:75719/1‑149 (MQ=255)                                                                        gcgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                                   <  5:190694/111‑1 (MQ=255)                                                                        gcgcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCAt                                                                                                                   >  6:190694/1‑111 (MQ=255)                                                                          gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAg                                                                                                                                                    >  6:366351/1‑77 (MQ=255)                                                                          gcTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAg                                                                                                                                                    <  5:366351/77‑1 (MQ=255)                                                                            tCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                            >  4:338775/1‑114 (MQ=255)                                                                            tCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGgt                                                                                                            <  3:338775/114‑1 (MQ=255)                                                                            tCCTTGCAGTCTGTCCTGGTTGGCATC‑AAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGCCCGCCAACTCATGACGGGACTATTCACTGTACTCAGCCCTGTCTAAGCCCGCct                                                                         >  6:398958/1‑148 (MQ=255)                                                                              cTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCctct                                                                       >  6:75031/1‑149 (MQ=255)                                                                                 gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                                     <  1:356166/148‑1 (MQ=255)                                                                                 gCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGc                                                                     <  4:75719/148‑1 (MQ=255)                                                                                         tCCTGGTTGGCATCAAAAAATGCCTCAACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCAGCCATCTCATGACGGGTGTATTCAGTGTGATGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCt                                                            >  5:393432/1‑149 (MQ=255)                                                                                              gTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                                         <  5:75031/147‑1 (MQ=255)                                                                                              gTTGGCATC‑AAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATATCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                                      <  5:398958/149‑1 (MQ=255)                                                                                                   cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                                                     >  5:344227/1‑98 (MQ=255)                                                                                                   cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCa                                                                                                     <  6:344227/98‑1 (MQ=255)                                                                                                   cATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAAc                                                  >  1:186731/1‑149 (MQ=255)                                                                                                     tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                                                >  8:364989/1‑149 (MQ=255)                                                                                                     tCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGc                                                >  6:194724/1‑149 (MQ=255)                                                                                                            aTGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                                         >  7:326682/1‑133 (MQ=255)                                                                                                            aTGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTc                                                         <  8:326682/133‑1 (MQ=255)                                                                                                            aTGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCt                                          >  3:238933/1‑148 (MQ=255)                                                                                                            aTGCCTCCACTTCGCCCGAGCTCGGCGC‑GTCGACTTCGCTGAGC‑AGGGAGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCt                                          <  4:238933/148‑1 (MQ=255)                                                                                                                cTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCa                                      >  5:16328/1‑148 (MQ=255)                                                                                                                    aCTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGCTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCACTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGc                                 >  8:157153/1‑149 (MQ=255)                                                                                                                      ttCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCg                               >  4:375205/1‑149 (MQ=255)                                                                                                                         gCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                                      <  3:222925/123‑1 (MQ=255)                                                                                                                         gCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGaa                                                      >  4:222925/1‑123 (MQ=255)                                                                                                                          cccTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc                           <  2:186731/149‑1 (MQ=255)                                                                                                                           ccctgctccGCGCAGTCGAC‑GCGCCGAGCGAGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTc                           >  7:270665/5‑149 (MQ=255)                                                                                                                            ctctctcgGCGC‑GTCTACTGCGCTGAGC‑AGGGCGTACACTTGCTGTGCCGCCATCTCATGCAGGGTGTATTCAGTGTACTGCGCCCttt                                                                                     <  8:355699/89‑3 (MQ=255)                                                                                                                            ctcgctcgGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                                     >  7:355693/1‑89 (MQ=255)                                                                                                                            ctcgctcgGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGt                                                                                     <  8:355693/89‑1 (MQ=255)                                                                                                                            atcgatcgTCGC‑GTCGACTGCGCTGAGA‑AGTGCGTACACTTGCTGGGCCGCCATATCATGCCGGGTGTATTCAGTGTACTGAGCCCTGt                                                                                     >  7:355699/2‑89 (MQ=255)                                                                                                                               gctcgGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtct                      >  8:300164/1‑149 (MQ=255)                                                                                                                               gctcgGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAt                        <  5:340439/147‑1 (MQ=255)                                                                                                                                ctcgGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGAtctc                     <  4:292790/149‑1 (MQ=255)                                                                                                                                             aCTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTgcgc         >  4:110909/1‑149 (MQ=255)                                                                                                                                                    tGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGc                                                              >  8:163803/1‑89 (MQ=255)                                                                                                                                                    tGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGc                                                              <  7:163803/89‑1 (MQ=255)                                                                                                                                                    tGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGACACt  >  4:262278/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      CCCGATGAGCGTTGCCCGAGCCAAGGCCAGGCCGTGTTCGCGTATGACCAATGGATTATCACTTTGGTTCTGCGCTCCTTGCAGTCTGTCCTGGTTGGCATCAAAAAATGCCTCCACTTCGCCCTCGCTCGGCGC‑GTCGACTGCGCTGAGC‑AGGGCGTACACTTGCTGGGCCGCCATCTCATGCCGGGTGTATTCAGTGTACTGCGCCCTGTCGAAGCCCGCCTCTGCCAGTCGCCTTTCGAATAACGCCGGGCTGCCAAATGCCGCTTCGATCTCGCCAACCTGCGCCGACACT  >  NC_002947/3948967‑3949261

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