Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	600,031	A→C	6.3%	C5G (TGC→GGC)	PP_0512 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	600,031	0	A	C	6.3%	233.8 / 4.8	75	C5G (TGC→GGC)	PP_0512	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (32/38);  new base C (3/2);  total (35/40)
Fisher's exact test for biased strand distribution p-value = 6.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.34e-01

CCAGTTCCACCTCGGCACCCGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGCCCTCCCTTTTTTT  >  NC_002947/599884‑600175                                                                                                                                                    |                                                                                                                                                  ccAGTTCCACCTCGGCACCCGCCGGCATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCAGGCCCGGGCGCCGGGTTGTCTGTCATTGGCCTGTCGATGCCGCAg                                                                                                                                                  >  2:253327/1‑149 (MQ=255)   aGTTCCACCTCGGCACCCGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAggc                                                                                                                                                >  1:117433/1‑149 (MQ=255)         aCCTCGGCACCCGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGgcaag                                                                                                                                          >  1:222985/1‑149 (MQ=255)           cTCGGCACCCGCCGGGACTCCCAGCCACGTGCCGCTCAGCTCGCGCTGGCAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCGGTCGCTGCCGCCTGCGgcagcga                                                                                                                                        <  8:208305/149‑5 (MQ=255)              ggCACCCGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTAGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAg                                                                                                                                     <  1:184358/149‑1 (MQ=255)               gCACCCGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAgc                                                                                                                                    <  3:49116/149‑1 (MQ=255)                  cccGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAgcggc                                                                                                                                 <  7:253041/149‑1 (MQ=255)                   ccgccgGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAgcggcg                                                                                                                                <  5:65258/149‑1 (MQ=255)                    cgccgGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGgcgc                                                                                                                               <  2:117433/149‑1 (MQ=255)                        gggATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCgcagca                                                                                                                           <  6:136486/149‑1 (MQ=255)                            tttCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGcagcacag                                                                                                                        >  3:151774/1‑148 (MQ=255)                                cAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCacga                                                                                                                   <  8:109206/149‑1 (MQ=255)                                cAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCacga                                                                                                                   >  6:59581/1‑149 (MQ=255)                                 aGCAAGGTGCGGCTCATCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTTTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCacgac                                                                                                                  <  6:72056/149‑1 (MQ=255)                                 aGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCacgac                                                                                                                  >  4:295917/1‑149 (MQ=255)                                           cGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCgg                                                                                                        >  6:193361/1‑149 (MQ=255)                                           cGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCgg                                                                                                        >  3:56372/1‑149 (MQ=255)                                             cTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAgg                                                                                                                                                 <  7:267139/106‑1 (MQ=255)                                             cTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAgg                                                                                                                                                 >  8:267139/1‑106 (MQ=255)                                             cTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCgg                                                                                                        >  6:150380/1‑147 (MQ=255)                                             cTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTa                                                                                                      <  3:84159/149‑1 (MQ=255)                                                aGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGt                                                                                                    <  3:295917/148‑1 (MQ=255)                                                   tCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑t                                                                                               >  2:36128/1‑149 (MQ=255)                                                    cgcgCTGGTAGGCCGGCAACGGGCCCGGTACTTTGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAgcggc                                                                                                                                 <  4:51081/115‑1 (MQ=255)                                                    cgcgCTGGTAGGCCGGCAACGGGCCCGGTACTTTGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAgcggc                                                                                                                                 >  3:51081/1‑115 (MQ=255)                                                    cgcgCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tg                                                                                              >  6:49554/1‑149 (MQ=255)                                                     gcgcTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tg                                                                                              <  6:132101/148‑1 (MQ=255)                                                     gcgcTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tg                                                                                              <  1:152025/148‑1 (MQ=255)                                                       gcTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tgggg                                                                                           <  1:36128/149‑1 (MQ=255)                                                       gcTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tgggg                                                                                           <  5:159850/149‑1 (MQ=255)                                                       gcTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tgggg                                                                                           <  4:172404/149‑1 (MQ=255)                                                           gTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCCAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tgggggctg                                                                                       >  3:25984/1‑146 (MQ=255)                                                           gTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tgggg                                                                                           >  6:117141/1‑145 (MQ=255)                                                           gTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑Tgggg                                                                                           <  5:117141/145‑1 (MQ=255)                                                           gTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGt                                                                                       >  6:49802/1‑149 (MQ=255)                                                            tAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTc                                                                                      <  7:279319/149‑1 (MQ=255)                                                              ggCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCt                                                                                    >  3:278842/1‑149 (MQ=255)                                                                      aCGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCCCTGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACggg                                                                             >  2:128924/1‑148 (MQ=255)                                                                      aCGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGc                                                                            <  5:193361/149‑1 (MQ=255)                                                                          gCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑T‑GGGGAGTACATGACGGGCaatat                                                                       >  3:151260/1‑145 (MQ=255)                                                                                   ttCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGGGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGa                                                                                         >  2:304444/1‑123 (MQ=255)                                                                                   ttCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGGGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGa                                                                                         <  1:304444/123‑1 (MQ=255)                                                                                   ttCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATg                                                                >  3:123372/1‑148 (MQ=255)                                                                                   ttCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGG‑CTGTAGCTTCCGCAGGCGGCAAGCAAAGCGGCGCATCACAGCACGACGAGCGGTAGGTTGTGCA‑TGGGGGAGTGCTTCAAGGGCATATTTACCaatga                                                              >  5:40520/1‑145 (MQ=255)                                                                                       gCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGt                                                           >  7:113829/1‑149 (MQ=255)                                                                                        ccTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGtt                                                          <  4:56372/149‑1 (MQ=255)                                                                                             gCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTTCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAAcc                                                      <  4:278842/148‑1 (MQ=255)                                                                                             gCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAAccc                                                     <  6:54013/149‑1 (MQ=255)                                                                                              ccGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAAcccc                                                    >  2:238284/1‑149 (MQ=255)                                                                                                 aGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCttt                                                 <  4:123372/149‑1 (MQ=255)                                                                                                   ggggCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTgg                                               <  8:158400/149‑1 (MQ=255)                                                                                                    gggCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGc                                              <  6:24402/149‑1 (MQ=255)                                                                                                         ggCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTc                                         >  5:174972/1‑149 (MQ=255)                                                                                                          gggcgGGGGCCGGGTTGTCGGCCTTTGGCCGGGCGCTGCCCCCGGGGCAAAGCAAAGGGGCGCAGCACCGCACGACGACCGCTCGGTAGTGCA‑T‑GGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCtt                                       <  4:25984/147‑1 (MQ=255)                                                                                                          gcgcGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGAAGCACAGCACGACGAGCGCTAGGTAGTGCA‑TGGGCGAGTCCTACACGGGCATATTTCACATGGATGTTAACACATTtgtcttgtct                                        >  8:51833/1‑149 (MQ=255)                                                                                                           cgcgGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCtt                                       <  5:59581/149‑1 (MQ=255)                                                                                                           cgcgGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCtt                                       <  4:151774/149‑1 (MQ=255)                                                                                                           cgcgGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCtt                                       <  8:113829/149‑1 (MQ=255)                                                                                                            gcgGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCATTGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCtt                                       <  1:238284/149‑1 (MQ=255)                                                                                                             cgGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGTCGAGCGCTCGGTAGTGCT‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGc                                     <  2:222985/149‑1 (MQ=255)                                                                                                             cgGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGc                                     <  5:49554/149‑1 (MQ=255)                                                                                                                  gcCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGc                                     <  4:83928/144‑1 (MQ=255)                                                                                                                  gcCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGc                                     >  3:83928/1‑144 (MQ=255)                                                                                                                        ttGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATgc                          <  2:75136/149‑1 (MQ=255)                                                                                                                           tCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAGTgcgc                        >  4:212193/1‑148 (MQ=255)                                                                                                                              gTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCtt                                       >  8:185712/1‑130 (MQ=255)                                                                                                                                   tGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGGA‑TGGGTGAGTCCTTGACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATgg                <  7:51833/148‑1 (MQ=255)                                                                                                                                       ctttgcttGCCGCCTGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGccctc           >  4:200388/8‑149 (MQ=255)                                                                                                                                       cTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGccctc           >  1:181221/1‑149 (MQ=255)                                                                                                                                       cTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGccctc           >  3:282054/1‑149 (MQ=255)                                                                                                                                          tgcttGCCGCCTGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGccctccct        >  2:152673/5‑149 (MQ=255)                                                                                                                                           gcttGCCGCCTGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGCCCTCCCtt       >  4:134900/4‑149 (MQ=255)                                                                                                                                              tGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGCCCTCCCttttt    >  1:94442/1‑149 (MQ=255)                                                                                                                                                ccGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGCGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGCCCTCCCttttttt  <  2:181221/149‑1 (MQ=255)                                                                                                                                                ccGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGCCCTCCCttttttt  <  6:174972/149‑1 (MQ=255)                                                                                                                                                    aGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTg                                <  3:189433/115‑1 (MQ=255)                                                                                                                                                    aGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTg                                >  4:189433/1‑115 (MQ=255)                                                                                                                                                    |                                                                                                                                                  CCAGTTCCACCTCGGCACCCGCCGGGATTTCCAGCAAGGTGCCGCTCAGCTCGCGCTGGTAGGCCGGCAACGGGCCCGGTACTTCGGCCTTGGCCGAGGGGGCTGGCGCGGGCGCCGGGTTGTCGGTCTTTGGCCTGTCGCTGCCGCAGGCGGCAAGCAAAGCGGCGCAGCACAGCACGACGAGCGCTCGGTAGTGCA‑TGGGGGAGTCCTTCACGGGCATATTTCCCATGGATGTTAACCCCTTTGGCTTGTCTTGCCAGTGCAATGCGCTACCATGGCCCTCCCTTTTTTT  >  NC_002947/599884‑600175

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