Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	332,991	0	G	T	20.7%	165.4 / 15.8	58	L207F (TTG→TTT)	pdeL	DNA‑binding transcriptional activator/c‑di‑GMP phosphodiesterase PdeL
Reads supporting (aligned to +/- strand):  ref base G (13/33);  new base T (12/0);  total (25/33)
Fisher's exact test for biased strand distribution p-value = 5.83e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTTCAT  >  NC_000913/332852‑333129                                                                                                                                            |                                                                                                                                           ctctGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTg                                                                                                                                            >  2:1241034/1‑140 (MQ=255)    tGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCa                                                                                                                                         >  1:938316/1‑140 (MQ=255)     ggCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAg                                                                                                                                        <  2:1354836/140‑1 (MQ=255)     ggCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAg                                                                                                                                        <  1:11934/140‑1 (MQ=255)         gAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCggg                                                                                                                                    >  1:635111/1‑140 (MQ=255)         gAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCggg                                                                                                                                    <  1:457539/140‑1 (MQ=255)         gAGCCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCggg                                                                                                                                    >  1:34555/1‑140 (MQ=255)           gTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGAt                                                                                                                                  <  1:4211/140‑1 (MQ=255)                 cGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGaaa                                                                                                                            <  1:1080412/140‑1 (MQ=255)                 aggTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGaaa                                                                                                                            <  1:498142/139‑1 (MQ=255)                    tCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGaaaaaa                                                                                                                         >  1:136087/1‑139 (MQ=255)                     cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  1:712139/140‑1 (MQ=255)                     cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:34555/140‑1 (MQ=255)                     cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  1:1241034/140‑1 (MQ=255)                     cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:938316/140‑1 (MQ=255)                     cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:635111/140‑1 (MQ=255)                     cTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:787460/140‑1 (MQ=255)                     cTTATTGTAATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAaga                                                                                                                        <  2:867994/140‑1 (MQ=255)                        aTTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGa                                                                                                                              <  2:1191453/131‑1 (MQ=255)                        aTTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGa                                                                                                                              >  1:1191453/1‑131 (MQ=255)                          tGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTc                                                                                                                   <  1:1235139/140‑1 (MQ=255)                              aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCCGCGGGGTTTGAAAAAGAGTGTCTgct                                                                                                               >  2:498142/1‑138 (MQ=255)                              aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGaa                                                                                                               <  1:502588/140‑1 (MQ=255)                              aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGTAAAAAAGGGTCTGaa                                                                                                               >  2:712139/1‑140 (MQ=255)                              aTAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGAAAAAAAGGGTCTGaa                                                                                                               >  2:1235139/1‑140 (MQ=255)                                     cccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGtt                                                                                                        >  1:10619/1‑140 (MQ=255)                                      ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgtgt                                                                                                                    <  2:164439/127‑1 (MQ=255)                                      ccGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAgtgt                                                                                                                    >  1:164439/1‑127 (MQ=255)                                        gCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTtaat                                                                                                     >  1:515564/1‑140 (MQ=255)                                         ccAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAGAGGGTCTGGACCTGGTtaata                                                                                                    >  2:502588/1‑140 (MQ=255)                                           aaCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAc                                                                                                              >  1:1212001/1‑128 (MQ=255)                                           aaCTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAAc                                                                                                              <  2:1212001/128‑1 (MQ=255)                                                 tgaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGt                                                                                            <  2:515564/140‑1 (MQ=255)                                                 tgaAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGt                                                                                            <  1:1486551/140‑1 (MQ=255)                                                     aCAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGAACCTGGTTAATAAATTAGGTaaag                                                                                        >  1:678987/1‑138 (MQ=255)                                                        gACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGAACCTGGGTAATAAATTAGGTAACGATa                                                                                     >  2:1143923/1‑140 (MQ=255)                                                          cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATaaa                                                                                   >  1:178818/1‑140 (MQ=255)                                                          cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGa                                                                                       >  1:388866/1‑136 (MQ=255)                                                          cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGa                                                                                       <  2:1424261/136‑1 (MQ=255)                                                          cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGa                                                                                       >  1:1424261/1‑136 (MQ=255)                                                          cTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGa                                                                                       <  2:388866/136‑1 (MQ=255)                                                                 tatTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGc                                                                            <  2:229846/140‑1 (MQ=255)                                                                            cGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtct                                                                     >  1:1298207/1‑136 (MQ=255)                                                                            cGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTtct                                                                     <  2:1298207/136‑1 (MQ=255)                                                                              gTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCt                                                               >  1:373944/1‑140 (MQ=255)                                                                               tAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTa                                                              >  1:1302505/1‑140 (MQ=255)                                                                                  aaCATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTTAAAAAAAGGGTCTGAACCTCGTTAATAAATTAGGGAACGATAAAAACAAGCCGGGTCTCGAGCTAACg                                                           >  2:97872/1‑140 (MQ=255)                                                                                       ttGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACg                                                      <  1:1143923/140‑1 (MQ=255)                                                                                        tgctgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGGAACGATAAAATCAAACTGGTTCTCGAGCTAACGGaaagg                                                     >  2:1035278/1‑137 (MQ=255)                                                                                         gctgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGGTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTa                                                    >  2:1340036/1‑140 (MQ=255)                                                                                           tgcCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAAc                                                  <  2:1302505/140‑1 (MQ=255)                                                                                                                  aTCAACGTCTCGGCGGGTTGTTTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGc                           >  2:965466/1‑140 (MQ=255)                                                                                                                       cGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag                      <  2:373944/140‑1 (MQ=255)                                                                                                                       cGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag                      <  2:10619/140‑1 (MQ=255)                                                                                                                       cGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCagag                      <  2:678987/140‑1 (MQ=255)                                                                                                                          ctcGGCGGGTTGTTTTTTTGCAGCGGGGTTTGAAAAAAAGGGTCTGGACCTGGTTAATAAATTAGGGAACGATAAAAACAAACTGGTTCTCGAGGTAACGGGAAGTAAACCTATTCCGGTAAAGCCAGAAGCCAGAGCgg                   >  2:1337446/1‑139 (MQ=255)                                                                                                                          ctcGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGGGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGa                   >  2:636525/1‑140 (MQ=255)                                                                                                                             ggcggGTTGTTTTTTTGCAGCGGGGTTTTAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGatat                >  2:1013791/1‑140 (MQ=255)                                                                                                                             ggcggGTTGTTTTTTGTCAGCGGGATTTGGAAAAAAGGGTCTGAACCTGGTTAATAAATTAGGGAAAGATAAAATCAAGCTGGTTCTCGAGGTAAAGGAAAGTAACCCTATTTCGG‑GACGCCAAAAGCCAGAGCGCTAtt               >  2:1469106/1‑140 (MQ=255)                                                                                                                                 ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa            <  1:383535/140‑1 (MQ=255)                                                                                                                                 ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa            <  1:179366/140‑1 (MQ=255)                                                                                                                                 ggTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGa            <  1:1264215/140‑1 (MQ=255)                                                                                                                                  gttgttTTTTTGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAAACCTATTCCGGTAACGCCAAAAGCCAGAGCGATATTTTAc           >  2:179366/1‑140 (MQ=255)                                                                                                                                           tGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTtcat  <  1:640645/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                                           CTCTGGCGGAGTCATCCGGTCTTATTGTCATAATGACCCGCCAACTGATGAAACAGACTGCGGATATTCTGATGCCGGTAAAACATTTGCTGCCGGACAATTTCCATATTGGCATCAACGTCTCGGCGGGTTGTTTTTTGGCAGCGGGATTTGAAAAAGAGTGTCTGAACCTGGTTAATAAATTAGGTAACGATAAAATCAAGCTGGTTCTCGAGCTAACGGAACGTAACCCTATTCCGGTAACGCCAGAAGCCAGAGCGATATTTGACAGCCTTCAT  >  NC_000913/332852‑333129

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