Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	321,043	C→G	6.3%	C145W (TGC→TGG)	aguA →	agmatine deiminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	321,043	0	C	G	6.3%	249.8 / 7.2	79	C145W (TGC→TGG)	aguA	agmatine deiminase
Reads supporting (aligned to +/- strand):  ref base C (34/40);  new base G (3/2);  total (37/42)
Fisher's exact test for biased strand distribution p-value = 6.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.37e-01

TACCTTCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCACCTGAGCCGCGAGCAGATTGAAGCGGT  >  NC_002947/320902‑321183                                                                                                                                              |                                                                                                                                              tACCTTCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTa                                                                                                                                       <  1:1895162/150‑1 (MQ=255)   ccTTCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTAcc                                                                                                                                     >  2:938621/1‑150 (MQ=255)     ttCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAg                                                                                                                                   >  2:1621589/1‑150 (MQ=255)     ttCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAg                                                                                                                                   >  1:1724981/1‑150 (MQ=255)      tCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGa                                                                                                                                  <  2:1399988/150‑1 (MQ=255)      tCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGa                                                                                                                                  <  1:2094054/150‑1 (MQ=255)       cGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGAc                                                                                                                                 <  1:245934/150‑1 (MQ=255)        gTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGAcc                                                                                                                                >  2:1355253/1‑150 (MQ=255)         tcatcaACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGGCAGCGCtccatt                                                                                                                                   >  2:1616098/1‑141 (MQ=255)         tcatcaACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGGCAGCGCtccatt                                                                                                                                   <  1:1616098/146‑6 (MQ=255)               aCGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCt                                                                                                                         >  2:1998267/1‑150 (MQ=255)                     aTGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGc                                                                                                                   <  1:1110336/150‑1 (MQ=255)                        gCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCt                                                                                                                                        >  2:227254/1‑126 (MQ=255)                        gCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCt                                                                                                                                        <  1:227254/126‑1 (MQ=255)                        gCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTgg                                                                                                                >  2:1625678/1‑150 (MQ=255)                           aaGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGGACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAg                                                                                                             <  1:1355253/150‑1 (MQ=255)                           aaGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAg                                                                                                             <  1:592569/150‑1 (MQ=255)                           aaGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAg                                                                                                             <  1:1998267/150‑1 (MQ=255)                              tGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAggtg                                                                                                          <  1:938621/150‑1 (MQ=255)                              tGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAggtg                                                                                                          <  1:1625678/150‑1 (MQ=255)                                 gcgGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTAc                                                                                                                                      >  2:1539513/1‑119 (MQ=255)                                 gcgGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTAc                                                                                                                                      <  1:1539513/119‑1 (MQ=255)                                      gtAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATc                                                                                                  <  1:651396/150‑1 (MQ=255)                                      gtAGATTGGGGCTGCAAAGCCTGGGGCGGCTTCGCTGCCGGCCTGTAGGCAGCCTGGAGCCGTGCCGAGCGGCTGGCCGCCCAGGTGCTGGAAATGGAGCGCTGCCATCGCTACCGGACAGAATGGTGTGTACGGGCACGTGGTTCGATc                                                                                                  >  2:794631/1‑150 (MQ=255)                                         gATTGGGGCTTCAAGACCTGGGGCGGCTTCGATGCCGGCCTGTACTCACCCTGGAACCGAGAAGATCGACGAGGCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCAc                                                                                               <  1:1621589/150‑1 (MQ=255)                                         gATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGAt                                                                                                   <  2:788665/146‑1 (MQ=255)                                         gATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGAt                                                                                                   >  1:788665/1‑146 (MQ=255)                                         gATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCAc                                                                                               <  2:1858829/150‑1 (MQ=255)                                         gATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCAc                                                                                               <  1:603438/150‑1 (MQ=255)                                         gATTGGGGCTTCAACGCCTGCGGCGGCTTCGATGGCGGTCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCAc                                                                                               <  2:631963/150‑1 (MQ=255)                                               ggCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGAc                                                                                         <  1:2181563/150‑1 (MQ=255)                                                  ttCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGt                                                                                      >  2:685799/1‑150 (MQ=255)                                                    cAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCa                                                                                                                                    >  2:7683/1‑102 (MQ=255)                                                    cAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCa                                                                                                                                    <  1:7683/102‑1 (MQ=255)                                                     aaCGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGaa                                                                                   >  1:1496575/1‑150 (MQ=255)                                                       cGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAgg                                                                                 <  2:1496575/150‑1 (MQ=255)                                                        gCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAggtg                                                                                                          >  1:2356790/1‑124 (MQ=255)                                                        gCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAggtg                                                                                                          <  2:2356790/124‑1 (MQ=255)                                                          cTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCAc                                                                              <  1:1733722/150‑1 (MQ=255)                                                          cTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCAc                                                                              >  1:1612013/1‑150 (MQ=255)                                                                  gcTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGTGCGCTGCCAGCGATACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATca                                                                      >  1:172005/1‑150 (MQ=255)                                                                   cTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATcac                                                                     >  2:185971/1‑150 (MQ=255)                                                                    ttCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGt                                                                           <  1:1724292/143‑1 (MQ=255)                                                                    ttCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGt                                                                           >  2:1724292/1‑143 (MQ=255)                                                                    ttCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATcacc                                                                    >  1:1486646/1‑150 (MQ=255)                                                                     tCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATcacca                                                                   <  1:185971/150‑1 (MQ=255)                                                                     tCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATcacca                                                                   <  2:1612013/150‑1 (MQ=255)                                                                        aTGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑g                                                               >  2:1410637/1‑150 (MQ=255)                                                                         tGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑Ga                                                              >  1:1384205/1‑150 (MQ=255)                                                                         tGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑Ga                                                              >  1:852870/1‑150 (MQ=255)                                                                          ggcggcCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAg                                                             <  2:1724981/150‑1 (MQ=255)                                                                              gcCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAAt                                                         <  1:2142438/150‑1 (MQ=255)                                                                              gcCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAAt                                                         <  2:1434051/150‑1 (MQ=255)                                                                              gcCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAAt                                                         <  2:422775/150‑1 (MQ=255)                                                                                cTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGc                                                       <  2:172005/150‑1 (MQ=255)                                                                                   tACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGc                                                                                                                          >  1:470076/1‑81 (MQ=255)                                                                                   tACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGc                                                                                                                          <  2:470076/81‑1 (MQ=255)                                                                                      gaaCCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGAGCTACCAGAGCGCAGGCTTTGTGCTGGAAGGTTGTTCGTTGCACGTAGACGGTGAGGGCACCGTGATCACCACT‑GAGGAATGCCTGCTc                                                 >  2:273216/3‑150 (MQ=255)                                                                                      gCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATc                                                                                                  <  1:869252/102‑1 (MQ=255)                                                                                      gCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATc                                                                                                  >  2:869252/1‑102 (MQ=255)                                                                                      gCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGt                                                                           >  1:496146/1‑125 (MQ=255)                                                                                      gCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGt                                                                           <  2:496146/125‑1 (MQ=255)                                                                                      gCACCCTGGAACCGCGACGAGGAACTGGACGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATC‑CAAGTGGAGGAATGCCTGCTc                                                 <  1:1771497/150‑1 (MQ=255)                                                                                        aCCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTcaa                                               >  2:1618096/1‑150 (MQ=255)                                                                                         cccTGCAACCGCGACGAGGAACTGGCCGCCAAGGTGCTTGAAATGGTGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATTCACGTAGACGGGGAAGGCACCGTGATCACCACT‑GAGGAATGACTGCTcaac                                              <  1:794631/150‑1 (MQ=255)                                                                                          ccaggAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTcaacc                                             >  1:2399823/4‑150 (MQ=255)                                                                                          ccTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGGCAGCGCtcc                                                                                                                                      <  2:1974380/62‑3 (MQ=255)                                                                                          ccTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGGCAGCGCtcc                                                                                                                                      >  1:1974380/1‑60 (MQ=255)                                                                                          ccTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCt                                                     >  2:626949/1‑142 (MQ=255)                                                                                          ccTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCt                                                     <  1:626949/142‑1 (MQ=255)                                                                                            tGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTcaaccgc                                           >  1:1095400/1‑150 (MQ=255)                                                                                                        gaGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCg                               >  2:949343/1‑150 (MQ=255)                                                                                                           gAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCAc                            <  1:685799/150‑1 (MQ=255)                                                                                                               tGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCACCtaa                        >  1:594301/1‑148 (MQ=255)                                                                                                                      cAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCACCTGAGCCGCGa                 <  2:852870/150‑1 (MQ=255)                                                                                                                         ggTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGACGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGaa                                                                                   >  1:183481/1‑82 (MQ=255)                                                                                                                         ggTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGaa                                                                                   <  2:183481/82‑1 (MQ=255)                                                                                                                              tGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCACCTCAGCCCCGGGCAGACTg         >  2:173480/1‑150 (MQ=255)                                                                                                                                     gtaGCGCTGGCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCACCTGAGCCGCGAGCAGATTGAAGCGGt  >  1:1648815/3‑150 (MQ=255)                                                                                                                                              ccAGCGCTACCAGACCGAAGGCTTTGTGCTGGCAGGTGGTTCg                                                                                                     >  1:1998429/1‑43 (MQ=255)                                                                                                                                              ccAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCg                                                                                                     <  2:1998429/43‑1 (MQ=255)                                                                                                                                              |                                                                                                                                              TACCTTCGTCATCAACGACCATGGCGAAGTGCGCGGTGTAGATTGGGGCTTCAACGCCTGGGGCGGCTTCGATGGCGGCCTGTACGCACCCTGGAACCGCGACGAGGAACTGGCCGCCAAGGTGCTGGAAATGGAGCGCTGCCAGCGCTACCAGACCGAAGGCTTTGTGCTGGAAGGTGGTTCGATCCACGTAGACGGTGAAGGCACCGTGATCACCACT‑GAGGAATGCCTGCTCAACCGCAACCGCAACCCGCACCTGAGCCGCGAGCAGATTGAAGCGGT  >  NC_002947/320902‑321183

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