Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	322,140	T→C	19.3%	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	19.3%	119.4 / 24.5	62	intergenic (+427/+31)	aguA/PP_0267	agmatine deiminase/ferric siderophore receptor
Reads supporting (aligned to +/- strand):  ref base T (24/26);  new base C (4/8);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 5.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.61e-01

TGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCG  >  NC_002947/321992‑322283                                                                                                                                                      |                                                                                                                                                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                                                            <  5:145582/140‑1 (MQ=255)                                                                                               aTGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATc                                                            >  6:145582/1‑140 (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                                            >  5:259825/1‑139 (MQ=255)                                                                                                                aTCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCg                                            <  6:259825/139‑1 (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)                                                                                                                                                gCATGCCCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgcc   <  6:35834/149‑1 (MQ=255)                                                                                                                                                 cATGATCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACAACAAGGCTACTCGGATCTCCCACCTGGATCTGCGCAGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATgccg  >  5:247551/1‑149 (MQ=255)                                                                                                                                                      |                                                                                                                                                TGGAAGCATCAGGGTCAGGACGACCGACCCACGACGTCTGCTTGCCGGGCGATACCCCCAAGGGTTGATCAGCGAAGCAGAGATGCAGGCGCTGATGGCGGTGTGGCACTGATCGATTGCACCT‑GTACCTGCCTCTTCGCGGGCATGCTCGCGAAGAGGCAGGTACAGGCAACTTAAAATCAGAAACTGTAGCTACCCGTCACCACAAGGCTACGCGGATCGCCCACCTGGATCTGCGCCGCACTGGTCGCCGAGCTGTAGTAGGTCTTGTCAGCAATGTTGCTCAATGCCG  >  NC_002947/321992‑322283

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