Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,970,038	T→G	6.8%	T84P (ACT→CCT)	fliC ←	flagellin

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,970,038	0	T	G	6.8%	206.2 / 3.6	75	T84P (ACT→CCT)	fliC	flagellin
Reads supporting (aligned to +/- strand):  ref base T (37/32);  new base G (3/3);  total (40/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAG  >  NC_002947/4969890‑4970187                                                                                                                                                     |                                                                                                                                                       ccGAAAGTGGTGGTGTTAGCAATACTGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTc                                                                                                                                                       <  1:921627/150‑1 (MQ=255)       gtggtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGCGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATg                                                                                                                                                 >  1:1691907/1‑150 (MQ=255)       gtggtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATg                                                                                                                                                 >  2:1828831/1‑150 (MQ=255)         ggtggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgc                                                                                                                                               <  2:254305/150‑1 (MQ=255)           tggtgTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGgcgc                                                                                                                                             <  1:395657/150‑1 (MQ=255)             gtgtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCt                                                                                                                                           >  1:1816558/1‑150 (MQ=255)               gtTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTc                                                                                                                                         <  1:1732530/150‑1 (MQ=255)                ttAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCttaa                                                                                                                                        >  1:1925365/1‑148 (MQ=255)                ttAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTca                                                                                                                                        >  1:748775/1‑150 (MQ=255)                ttAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTca                                                                                                                                        >  2:1071162/1‑150 (MQ=255)                ttAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCATTATGGCCCCATGGCTCCTTca                                                                                                                                        <  2:1418583/150‑1 (MQ=255)                ttAGCAATACGGGTGAGATCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCTGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTca                                                                                                                                        >  2:962884/1‑150 (MQ=255)                   gCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTGCTGCTGGACGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTca‑gca                                                                                                                                    >  1:38600/1‑150 (MQ=255)                            gTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACTCTTCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCa                                                                                                                           >  1:1367269/1‑150 (MQ=255)                                 cTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATaga                                                                                                                      >  1:1141926/1‑150 (MQ=255)                                    gCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGAt                                                                                                                   >  2:1320590/1‑150 (MQ=255)                                      ccGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATAc                                                                                                                 <  1:1336200/150‑1 (MQ=255)                                        gACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCa                                                                                                               >  1:575139/1‑150 (MQ=255)                                        gACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCa                                                                                                               <  1:303119/150‑1 (MQ=255)                                          cAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATc                                                                                                             <  1:1011680/150‑1 (MQ=255)                                           aGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCg                                                                                                            >  1:1909709/1‑150 (MQ=255)                                                ggTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGc                                                                                                       >  1:1075901/1‑150 (MQ=255)                                                 gTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCg                                                                                                      >  1:1494299/1‑150 (MQ=255)                                                  tAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGt                                                                                                     <  2:776186/150‑1 (MQ=255)                                                    aaCTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTc                                                                                                   <  2:575139/150‑1 (MQ=255)                                                      cTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCtt                                                                                                 <  2:1925365/150‑1 (MQ=255)                                                        ccggctgcATGGGGGCGCGGTCGTCAGCGCTGTTCGAGCCGTGCGCCGATTGCAGGGCCAGTGCACGCATACGCTTCATCATGTTGCCTGAAGGCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGa                                                                                               <  1:69305/147‑1 (MQ=255)                                                        ccTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGa                                                                                               <  2:1816558/150‑1 (MQ=255)                                                         ctgctgCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAc                                                                                              >  2:1625912/1‑150 (MQ=255)                                                          tgctgcATGGAGGCGCGGTCGTCATCGCTGTGCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGAcc                                                                                             <  2:138254/150‑1 (MQ=255)                                                              gcATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCaa                                                                                         >  1:780347/1‑150 (MQ=255)                                                                 tGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGc                                                                                      >  2:51599/1‑150 (MQ=255)                                                                       cgcgGTCGTCCTCGCTGTTCGACACGTTGATCCACTACATTTCCATGTCAGGCGTACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGCTACCATCGTTGGCGTTCTTGACCGCAACGCTCAGAc                                                                                <  2:1646518/150‑1 (MQ=255)                                                                           gtcgtcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCttt                                                                            <  2:1717383/150‑1 (MQ=255)                                                                            tcgtcATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttg                                                                           >  1:577274/1‑150 (MQ=255)                                                                                 aTCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttgatttg                                                                      >  2:545545/1‑150 (MQ=255)                                                                                  tCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttgatttga                                                                     <  2:1092252/150‑1 (MQ=255)                                                                                    gCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGt                                                                   <  1:51599/150‑1 (MQ=255)                                                                                      tGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTgg                                                                 <  2:780347/150‑1 (MQ=255)                                                                                        ttCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACTCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTc                                                               >  1:1650349/1‑150 (MQ=255)                                                                                         tCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCa                                                              <  2:1691907/150‑1 (MQ=255)                                                                                                gTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGtt                                                       >  2:724263/1‑150 (MQ=255)                                                                                                    gCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTGGATTTGCGTGGTCAGACGGTTGgag                                                   >  2:714511/1‑150 (MQ=255)                                                                                                     ccGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTg                                                      <  2:774294/146‑1 (MQ=255)                                                                                                     ccGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTg                                                      >  1:774294/1‑146 (MQ=255)                                                                                                     ccGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGgaga                                                  >  2:1819959/1‑150 (MQ=255)                                                                                                            cAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAg                                           <  1:454997/150‑1 (MQ=255)                                                                                                            cAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAg                                           >  1:1236842/1‑150 (MQ=255)                                                                                                                gacAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCt                                       >  2:651167/3‑150 (MQ=255)                                                                                                                gCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttga                                                                          >  1:46920/1‑115 (MQ=255)                                                                                                                gCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttga                                                                          <  2:46920/115‑1 (MQ=255)                                                                                                                    gTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTgcgg                                   <  1:962884/150‑1 (MQ=255)                                                                                                                        aCGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGgag                                                   >  1:1665096/1‑130 (MQ=255)                                                                                                                        aCGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGgag                                                   <  2:1665096/130‑1 (MQ=255)                                                                                                                          tcATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtcgt                             <  2:1650349/149‑1 (MQ=255)                                                                                                                          gCANACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtcgt                             >  2:1191719/1‑150 (MQ=255)                                                                                                                           cATACGCTGCATAATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCgtcgtc                            >  2:1715703/1‑150 (MQ=255)                                                                                                                             tACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCtt                          <  2:748775/150‑1 (MQ=255)                                                                                                                                gCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACAATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTgg                       >  1:1163061/1‑150 (MQ=255)                                                                                                                                       gatgTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGtt                <  2:177951/150‑1 (MQ=255)                                                                                                                                             tgcTGAAGGCGCCATGGCGACTTCAGGCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGAttt           >  2:1666433/3‑150 (MQ=255)                                                                                                                                              ccTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCt                                       <  2:210198/120‑1 (MQ=255)                                                                                                                                              ccTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCt                                       >  1:210198/1‑120 (MQ=255)                                                                                                                                              ccTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                               >  2:780520/1‑128 (MQ=255)                                                                                                                                              ccTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGcgtc                               <  1:780520/128‑1 (MQ=255)                                                                                                                                               cTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCt                                                                                                  <  2:84851/60‑1 (MQ=255)                                                                                                                                               cTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCt                                                                                                  >  1:84851/1‑60 (MQ=255)                                                                                                                                               cTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCa        >  2:1145073/1‑150 (MQ=255)                                                                                                                                                tGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCTCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg       <  2:38600/150‑1 (MQ=255)                                                                                                                                                tGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg       <  1:545545/150‑1 (MQ=255)                                                                                                                                                tGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg       >  1:589095/1‑150 (MQ=255)                                                                                                                                                tGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg       <  2:1141926/150‑1 (MQ=255)                                                                                                                                                tGAAGGCGCCATGGCGACTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg       >  1:1577499/1‑150 (MQ=255)                                                                                                                                                tGAAGGCGCCATGGCGACTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAg       <  2:1577499/150‑1 (MQ=255)                                                                                                                                                  aaGGCGCCATGGCGACTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCtttga                                                                          >  1:897681/1‑81 (MQ=255)                                                                                                                                                     tCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCag  <  2:1075901/150‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                       CCGAAAGTGGTGGTGTTAGCAATACGGGTGAGCTCGCCCGACAGCGCGGTAAACTCCTGCTGCATGGAGGCGCGGTCGTCATCGCTGTTCGAGCCGTTCGCCGATTGCAGGGCCAGTTCACGCATACGCTGCATGATGTTGCCTGAAGTCGCCATGGCGCCTTCA‑GCAGTTTGCGCAATAGAGATACCATCGTTGGCGTTCTTGACCGCAACGCTCAGACCTTTGATTTGAGTGGTCAGACGGTTGGAGATCTGCAGGCCTGCGGCGTCGTCTTTGGCACTGTTGATTTTCAGGCCAG  >  NC_002947/4969890‑4970187

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