Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	322,140	T→C	18.0%	intergenic (+427/+31)	aguA → / ← PP_0267	agmatine deiminase/ferric siderophore receptor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	322,140	0	T	C	18.0%	129.3 / 22.0	64	intergenic (+427/+31)	aguA/PP_0267	agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (24/27);  new base C (7/4);  total (32/32)
Fisher's exact test for biased strand distribution p-value = 5.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.94e-01

CATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGC  >  NC_002947/321998‑322288                                                                                                                                                |                                                                                                                                                     cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                <  1:106253/149‑1 (MQ=255) cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                <  5:36294/149‑1 (MQ=255) cATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGaa                                                                                                                                                <  5:12227/149‑1 (MQ=255)   tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGTGCGATACGCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                              <  2:139903/149‑1 (MQ=255)   tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                >  4:127363/1‑147 (MQ=255)   tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                              <  8:37436/149‑1 (MQ=255)       ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCa                                                                                                                                          <  2:303531/149‑1 (MQ=255)       ggTCAGCTCGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCATGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                              <  1:166382/145‑1 (MQ=255)            ggACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTAc                                                                                                                                     <  3:179087/149‑1 (MQ=255)             gacgacCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTa                                                                                                                                      >  1:40949/1‑147 (MQ=255)             gacgacCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTa                                                                                                                                      >  3:80409/1‑147 (MQ=255)                 accgaccCTCGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGGAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgc                                                                                                                                                   <  5:184928/130‑1 (MQ=255)                 accgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTcgc                                                                                                                                                   >  6:184928/1‑130 (MQ=255)                    gaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACt                                                                                                                             >  4:347536/1‑149 (MQ=255)                        cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCCCCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaaa                                                                                                                         <  6:295918/149‑1 (MQ=255)                             cgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                                    <  2:209349/149‑1 (MQ=255)                              gTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGa                                                                                                                   <  3:127363/149‑1 (MQ=255)                                   cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                                    >  2:20450/1‑143 (MQ=255)                                   cTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAg                                                                                                                    <  1:20450/143‑1 (MQ=255)                                    ttGTCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGt                                                                                                             <  4:273385/149‑1 (MQ=255)                                      gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                            <  3:347536/148‑1 (MQ=255)                                      gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                            <  5:307499/148‑1 (MQ=255)                                      gCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTa                                                                                                            <  4:80409/148‑1 (MQ=255)                                          ggCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTAc                                                                                                       <  6:106678/149‑1 (MQ=255)                                                  ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                >  7:290619/1‑100 (MQ=255)                                                  ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGaa                                                                                                                                                <  8:290619/100‑1 (MQ=255)                                                  ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                              <  3:267186/102‑1 (MQ=255)                                                  ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                              >  4:267186/1‑102 (MQ=255)                                                    ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                              >  4:268179/1‑148 (MQ=255)                                                     cAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcaccac                                                                                             >  2:16373/1‑148 (MQ=255)                                                         ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                              <  7:109559/95‑1 (MQ=255)                                                         ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                              >  8:109559/1‑95 (MQ=255)                                                         ggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCACGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGc                                                                                        >  7:32387/1‑149 (MQ=255)                                                                         aaccaggCAGGCGCGGATGCGGGTGTGGCACTTATC‑GATTGCACCTGTACCTGCCGCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTAAAAATCACAAACTGAAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAc                                                                        <  4:329125/143‑1 (MQ=255)                                                                          gagaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGGATTGCACCTGTACCTGCCTCTTCGCGGGCATGATCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGGTACCCTACACCACAAGGCTACGCGGAACGCCCAc                                                                        >  1:250130/1‑149 (MQ=38)                                                                                   gcgcTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGATACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgc                                                              >  8:454/1‑149 (MQ=255)                                                                                   gcgcTGATGGCGGGGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTgc                                                              <  2:250130/149‑1 (MQ=255)                                                                                          tGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCAc                                                       <  2:40949/149‑1 (MQ=255)                                                                                                      aCTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGACAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGc                                           >  7:223824/1‑149 (MQ=255)                                                                                                      aCTGATC‑GAGCGCACCTGCACCTGCCTCTCCGCGGGCATGCACCCGAAGAGGCAGGTACAGGCAACTTAAAATCACAAACTGTAGCTACCCGTCACCACAAGGCTCCGCGGATCGCCCACCTGCATCTGCGCCGCACTGCGCGCCGAGc                                           <  8:32387/149‑1 (MQ=255)                                                                                                              gATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgtagta                                    <  4:225833/149‑1 (MQ=255)                                                                                                               aTTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGtagtag                                   <  3:268179/149‑1 (MQ=255)                                                                                                                      cTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGt                            <  3:173/149‑1 (MQ=255)                                                                                                                      cTGTACCTACCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGt                            >  6:39533/1‑149 (MQ=255)                                                                                                                       tGTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTc                           >  3:118757/1‑149 (MQ=255)                                                                                                                           ccTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGc                        >  1:191039/1‑148 (MQ=255)                                                                                                                            cTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGCTCGCACACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa                      >  4:85614/1‑149 (MQ=255)                                                                                                                             tGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAAt                     >  6:356059/1‑149 (MQ=255)                                                                                                                              gCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATg                    >  7:253529/1‑149 (MQ=255)                                                                                                                                 tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                 >  3:22877/1‑149 (MQ=255)                                                                                                                                 tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                 >  3:191740/1‑149 (MQ=255)                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGTAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                 >  8:310157/1‑149 (MQ=255)                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTc                           <  6:132805/139‑1 (MQ=255)                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTc                           >  5:132805/1‑139 (MQ=255)                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                 >  6:253375/1‑149 (MQ=255)                                                                                                                                 tctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                 >  6:14809/1‑149 (MQ=255)                                                                                                                                  ctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc                >  3:98648/1‑149 (MQ=255)                                                                                                                                   ttCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGTTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTg                 <  7:454/147‑1 (MQ=255)                                                                                                                                   ttCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc                >  1:342315/1‑148 (MQ=255)                                                                                                                                     cgcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt               <  1:51224/147‑1 (MQ=255)                                                                                                                                     cgcgGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCt               >  2:51224/1‑147 (MQ=255)                                                                                                                                        gggCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATg          >  4:254192/1‑149 (MQ=255)                                                                                                                                            aTGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgccgc      >  8:120710/1‑149 (MQ=255)                                                                                                                                               cTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGccgccct   <  4:98648/149‑2 (MQ=255)                                                                                                                                                tCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGccgcccgc  >  7:292493/1‑149 (MQ=255)                                                                                                                                                |                                                                                                                                                     CATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATC‑GATTGCACCTGTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCGCCCGC  >  NC_002947/321998‑322288

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