Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	322,133	G→A	18.0%	intergenic (+420/+38)	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,133	0	G	A	18.0%	130.6 / 19.3	61	intergenic (+420/+38)	aguA/PP_0267	agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base G (24/26);  new base A (4/7);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 5.26e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

TGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGC  >  NC_002947/321992‑322274                                                                                                                                               |                                                                                                                                              tGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCt                                                                                                                                        >  1:7979/1‑149 (MQ=255)    aaGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCcgc                                                                                                                                     >  3:79173/1‑149 (MQ=255)     aGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTcgc                                                                                                                                     >  4:335048/1‑148 (MQ=255)        aTCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTcgcg                                                                                                                                    >  6:311411/1‑146 (MQ=255)         tccgggTCAGGACGCCCGACCCACGACGTCTGCTTGCCGGGCGATCCCCCCAAGGGTTGATCCGCTAAGCAGAGATGCAGGCGCTGATGGGGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCTAaga                                                                                                                                <  6:9106/146‑1 (MQ=255)         tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAaga                                                                                                                                <  5:311411/149‑1 (MQ=255)         tCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                <  4:79173/149‑1 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCTAGTACCTGCCTCTTCGCGAGCATGCCCGCGAaga                                                                                                                                <  1:217700/148‑1 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAagag                                                                                                                               <  6:124814/148‑1 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAgg                                                                                                                              >  2:146645/1‑149 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAgg                                                                                                                              <  4:333254/149‑1 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCGAagag                                                                                                                               <  5:144739/148‑1 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCGAag                                                                                                                                 >  4:221037/1‑146 (MQ=255)           aGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCGAag                                                                                                                                 <  3:221037/146‑1 (MQ=255)            gggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGc                                                                                                                             <  8:321234/149‑1 (MQ=255)             ggTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCa                                                                                                                            <  2:7979/149‑1 (MQ=255)              gTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGAGCATGCCCGCGAAGAGGCAg                                                                                                                           <  7:275917/149‑1 (MQ=255)               tCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAgg                                                                                                                          <  1:146645/149‑1 (MQ=255)                        ccgaccCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCaa                                                                                                                 <  3:240949/149‑1 (MQ=255)                             ccACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCACGAAGAGGCAGGTACAGGCAACTTa                                                                                                             >  1:232260/1‑148 (MQ=255)                              cACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTaa                                                                                                            >  1:75624/1‑148 (MQ=255)                               acgacgTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCACGAAGAGGCAGGTACAGGCAACTTaaaa                                                                                                          <  2:232260/149‑1 (MQ=255)                                             gcGGGCGATACCCCCAAGTGTTGATGAGCGATGCAGAGTTGCAGTCGCTGATGGCGGTGTGGCACTGATCGTTTGCACCT‑GTACCTGCCTGTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAATGTGTAg                                                                                             <  3:233266/147‑1 (MQ=38)                                             ccGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGc                                                                                            <  2:75624/149‑1 (MQ=255)                                                       cccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTca                                                                                   <  3:335048/148‑1 (MQ=255)                                                        ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACt                                                                                                 >  8:258243/1‑133 (MQ=255)                                                        ccccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACt                                                                                                 <  7:258243/133‑1 (MQ=255)                                                         cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                           >  2:310975/1‑106 (MQ=255)                                                         cccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAg                                                                                                                           <  1:310975/106‑1 (MQ=255)                                                          ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                >  8:2038/1‑148 (MQ=255)                                                          ccAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTcacca                                                                                >  7:229648/1‑148 (MQ=255)                                                            aaGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCCTGCTCGCGAAGAGGCAGGTACAGGCAAGTTAAAATCAGATACTGTAGCTACCCGTcaccac                                                                               >  2:200397/1‑147 (MQ=255)                                                             aGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAg                                                                            <  7:2038/149‑1 (MQ=255)                                                              gggTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCa                                                                                                       <  4:97995/121‑1 (MQ=255)                                                                gTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCAGTCACCACAAGGc                                                                          >  8:146915/1‑148 (MQ=255)                                                                     tCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTAcgcg                                                                    >  1:214781/1‑149 (MQ=255)                                                                       aGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGa                                                                  <  8:229648/149‑1 (MQ=255)                                                                         cGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATc                                                                >  8:7133/1‑149 (MQ=255)                                                                             gCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGccc                                                            >  5:251613/1‑149 (MQ=255)                                                                               agagATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACAACAAGGCTACGCGGATCGCCCAc                                                          >  8:248594/1‑149 (MQ=255)                                                                                gagaTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCAcc                                                         >  6:340249/1‑149 (MQ=255)                                                                                   aTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTgg                                                      <  6:251613/149‑1 (MQ=255)                                                                                   aTGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTgg                                                      <  7:7133/149‑1 (MQ=255)                                                                                             tGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATc                                                                                                        >  1:194473/1‑89 (MQ=255)                                                                                             tGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATc                                                                                                        <  2:194473/89‑1 (MQ=255)                                                                                               aTGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAg                                                                            <  6:277824/115‑1 (MQ=255)                                                                                               aTGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAg                                                                            >  5:277824/1‑115 (MQ=255)                                                                                               aTGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATc                                                   >  6:145582/1‑140 (MQ=255)                                                                                               aTGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATc                                                   <  5:145582/140‑1 (MQ=255)                                                                                                    ggTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGt                                     <  7:146915/149‑1 (MQ=255)                                                                                                     gtgtGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTc                                    <  2:214781/149‑1 (MQ=255)                                                                                                           cACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTCTAGCTACCCGTCACCACAAGGCTACGAGGATCGCCCACCTGGATCTGcgcc                                             >  1:138689/1‑134 (MQ=255)                                                                                                                aTCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCg                                   <  6:259825/139‑1 (MQ=255)                                                                                                                aTCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCg                                   >  5:259825/1‑139 (MQ=255)                                                                                                                aTCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTgta                         <  5:340249/149‑1 (MQ=255)                                                                                                                 tCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTg                           >  7:81796/1‑146 (MQ=255)                                                                                                                 tCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTg                           <  8:81796/146‑1 (MQ=255)                                                                                                                                    gCCTCTTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATg      >  5:35834/1‑149 (MQ=255)                                                                                                                                       tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa        <  5:325292/144‑1 (MQ=255)                                                                                                                                       tctTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCaa        >  6:325292/1‑144 (MQ=255)                                                                                                                                        ctTCGCGAGCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGc  >  5:227595/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                              TGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGC  >  NC_002947/321992‑322274

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