Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,206,435	0	C	T	51.6%	38.6 / 95.1	64	intergenic (‑179/‑28)	ttdR/ttdA	transcriptional activator of ttdABT/L‑tartrate dehydratase, alpha subunit
Reads supporting (aligned to +/- strand):  ref base C (3/28);  new base T (33/0);  total (36/28)
Fisher's exact test for biased strand distribution p-value = 4.02e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.28e-08

TTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAA  >  NC_000913/3206293‑3206583                                                                                                                                               |                                                                                                                                                      ttCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaa                                                                                                                                                <  1:326787/150‑1 (MQ=255) ttCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACaa                                                                                                                                                <  2:877292/150‑1 (MQ=255)   cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACaaa                                                                                                                                               >  1:208970/1‑149 (MQ=255)   cGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGAAGTTTGTGGCGGTGAATTTTTCCAGACAAATa                                                                                                                                             <  1:250441/151‑1 (MQ=255)     cAAAACATTGCCATTCCCTGCTCTGGCGTAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACa                                                                                                                                           <  2:954121/151‑1 (MQ=255)      aaaaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAATa                                                                                                                                             >  1:776824/1‑148 (MQ=255)        aaCATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTTCCACAAAAACaaa                                                                                                                                         >  2:176583/1‑150 (MQ=255)           aTTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAcaaaaacaaaaa                                                                                                                                       >  1:187983/1‑149 (MQ=255)           aTTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACt                                                                                                                                     <  2:81772/151‑1 (MQ=255)            ttGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAATTTGTGGCGGTGGATTTTTTCAGGAAAAAACaaaa                                                                                                                                        >  2:93664/1‑147 (MQ=255)             tGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATg                                                                                                                                    >  2:94105/1‑150 (MQ=255)                 aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑Gtt                                                                                                                              <  1:416431/151‑1 (MQ=255)                 aTTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAATACaaaaaatcgggtt                                                                                                                               >  2:633507/1‑143 (MQ=255)                    cccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGcc                                                                                                                           <  1:854335/151‑1 (MQ=255)                     ccTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCAACAAATACAAAAAATGGA‑Gtttgc                                                                                                                           >  2:548013/1‑147 (MQ=255)                       tGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATACAAAAACttgtggt                                                                                                                               >  2:818050/1‑139 (MQ=255)                        gCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTTATACAAAAAAAAAAAATGGA‑GTTGCCat                                                                                                                         >  2:619083/1‑149 (MQ=255)                         ctctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAACGGA‑Gt                                                                                                                               >  1:271062/1‑142 (MQ=255)                          tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTTTATTTTTTTTCACAAAAAAAAAAAATAGA‑GTTgg                                                                                                                            >  1:647575/1‑143 (MQ=255)                          tctGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAaaaaaaaaaaatagggtt                                                                                                                               >  2:671190/1‑134 (MQ=255)                            tGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATAAAAAAACTGGA‑GTTGGCatgat                                                                                                                      >  1:609866/1‑148 (MQ=255)                                     attcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGtaa                                                                                                          <  1:298102/151‑1 (MQ=255)                                       tcattcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGtaat                                                                                                         <  1:916231/150‑1 (MQ=255)                                         attcTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAg                                                                                                      <  2:683869/151‑1 (MQ=255)                                             ttATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAACaa                                                                                                   >  1:270868/1‑150 (MQ=255)                                              tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCa                                                                                                    <  1:931155/148‑1 (MQ=255)                                              tATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAAc                                                                                                  <  2:576981/150‑1 (MQ=255)                                               aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATAAAAAAAATGGA‑GGTGCCATGATGAGCGGAAGGAATAAGCaa                                                                                                   >  1:193803/1‑148 (MQ=255)                                               aTTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACa                                                                                                 <  1:547668/150‑1 (MQ=255)                                                 tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAAAAAAAAACTGA‑GGTGCCAGTATGAGGGGAAGAAAAAAGCAACa                                                                                                 >  2:361501/1‑148 (MQ=255)                                                 tttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAGACAAAAACAAAAAATGGA‑GTTTCCATGATGAGGGCAAATAATAAGCAACa                                                                                                 >  1:497317/1‑148 (MQ=255)                                                  ttCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGc                                                                                              <  2:325225/150‑1 (MQ=255)                                                    cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAAACAAATACAAAAAATGGA‑GTTGCCATGATGAGCGAAAGGAATAAACAACAAGc                                                                                              >  2:420050/1‑148 (MQ=255)                                                    cgcgTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGt                                                                                           <  1:626700/151‑1 (MQ=255)                                                        ttGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAAACAAATACAAAAACTTGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGt                                                                                           <  2:87184/147‑1 (MQ=255)                                                                cGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGA‑GGTGGCATGATGAGCGAAAGTAATAAACAAAAGGCAGGGAATAAGtt                                                                                  >  1:25616/1‑148 (MQ=255)                                                                 gAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAAAAAAAAAAATGGA‑GGTGGCCTGATGAGGGGAAAGAAAAAACAAAAAGCAGGGAaaaa                                                                                     >  2:362822/1‑144 (MQ=255)                                                                  aaCTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGAATTTTTCAAACAAAAAAAAAAAATGGG‑GGTTGCATGATGAGGGGAAATAATAAACAAAAAGCAGGGaa                                                                                        >  2:48227/1‑140 (MQ=255)                                                                      gTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAACAAAAAATGGA‑GTTTCCATGATGAGCGAAAGTAATAAACAAAAAGCAGGGAAAAAATTGACa                                                                              >  1:103854/1‑146 (MQ=255)                                                                        tAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGAt                                                                         <  2:400029/149‑1 (MQ=255)                                                                          aTTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAAAAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGATGACAGAGATTGt                                                                      >  1:663384/1‑150 (MQ=255)                                                                            tACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAATACAAAAAATGGA‑GTTTGCATGGTGAGGGGAAGAAATAAGCAACAAGCAGTGAAAAAATTGACAGAGAATGt                                                                      >  2:167201/1‑148 (MQ=255)                                                                                  ttCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAATTTTTTCTCCCAAATAAAAAAAATGGA‑GGTGGCATGATGAGCGAAAATTATAAGCAAAAGGCAGTGAATAAGGTTAGAGAGAATGTTGCTAACtt                                                             >  1:578434/1‑151 (MQ=255)                                                                                   tCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACtt                                                             <  1:108969/150‑1 (MQ=255)                                                                                    cACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTa                                                           <  1:11242/151‑1 (MQ=255)                                                                                       ggTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCg                                                        <  1:324534/151‑1 (MQ=255)                                                                                       ggTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGCGAAATTTTTCACACAAATATAAAAAATGGG‑GTTTTCATGAGTAGCGAAAGAAAAAAACAAAAAGCAATGAAAAAATTGACAGAGATTGTTGCTAAATTTTcc                                                         >  2:143017/1‑150 (MQ=255)                                                                                        gTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGAAAAAAAAAAAAAATGGA‑GTTGCCATGATGAGCGAAAGAAAAAAAAAAAAGGCAGGGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGc                                                       >  1:640929/1‑151 (MQ=255)                                                                                         tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGGATTTTTCCAAAAAAATAAAAAAAATTGAGGTTGCC‑TGATGAGGGAAAAGAATAAACAAAAAGCAGGGAATAAGATGAAAGAGAATGTCGCTAACTTTTCcccc                                                      >  2:886157/1‑151 (MQ=255)                                                                                         tGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGACAAAAAAAAAAAATGGA‑GTTGTCATGATGAGGGAAAATAATAAACAAAAAGCAGTGAATAAGTTTAAAGAGATTGTCGCTAACTTTAACGcc                                                      >  2:470418/1‑151 (MQ=255)                                                                                          ggCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGcc                                                      <  2:459965/150‑1 (MQ=255)                                                                                                aaaTATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAAATTTTTCAAACAAATAAAAAAAATGGG‑GTTTGCATGATGAGCGGAAGTAATAAACAAAAAGGAGGGAATAAATTGACAGAGAATGTCGCTAACTTTACCGCCATGAttt                                               >  2:118045/1‑151 (MQ=255)                                                                                                  atatAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCt                                             <  1:615333/151‑1 (MQ=255)                                                                                                   tataAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAATACAAAAAATGGA‑GTTGGCCTGATGAGCGAAAGGAATAAACAACAGGCAGGGGATAAGTTTACAGAGAATTTCGCTAACTTTAACGCCATGAttttt                                             >  1:191431/1‑148 (MQ=255)                                                                                                        aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCCGACAAAAACAAAAAATGGA‑GGTGCCATGGTGGGCGGAAGGAATAAGGAAAAGGGAGGGGATAAGGTGACAGAGATTGTCGGTAACTTTTCCGCCCTGATTTTTAccc                                         >  1:61974/1‑149 (MQ=255)                                                                                                        aaCCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTTCAGAAAAAAACAAAAAATGGA‑GTTGTCCTGATGTGGGGAAGTAATAAAAAAAAGGCAGTGGATAAAGTGAAAGAGATTGTCGCTGAATTTTTCCCCATGCTTTTTAccc                                         >  2:38707/1‑149 (MQ=255)                                                                                                                ttttAGTGGTAGTTTGGGGCGGTAAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGAGGAGCAAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtg                                <  2:190657/150‑1 (MQ=255)                                                                                                                  ttAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgat                              <  1:776860/150‑1 (MQ=255)                                                                                                                       ggaaGTTTGTGGAGGTGAATTTTTCCAAACAAATAAAAAAACGGGA‑GTTGCAAGGATGAGAAAAAGTAAGAAACAACGGGCAGTGAAAAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCtgatga                            <  1:949202/144‑1 (MQ=255)                                                                                                                            tttGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGtggtgg                    <  2:843635/150‑1 (MQ=255)                                                                                                                                tgGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATaa                >  1:217260/1‑150 (MQ=255)                                                                                                                                  gCGGTGAATTTTTCCAGACAAATACAAAAACGGGA‑GTTGCCATGATGAGCGAAATTAATAAGCAACAGGCAGTGAAAAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACt             <  1:196327/151‑1 (MQ=255)                                                                                                                                        aaTTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAAc        >  2:537152/1‑150 (MQ=255)                                                                                                                                        aaTTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTAAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTaaa          <  2:664413/148‑1 (MQ=255)                                                                                                                                               ccAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTaa  >  1:653618/1‑149 (MQ=255)                                                                                                                                               |                                                                                                                                                      TTCGCAAAACATTGCCATTCCCTGCTCTGGCGCAATATTCATTCTTATTTTCGCGTTGATTCACGAACTGTTAATTACATATTCACGGTGGCAAAAAATATAAAACCACATTTTTAGTGGTAGTTTGTGGCGGTGAATTTTTCCAGACAAATACAAAAACTGGA‑GTTGCCATGATGAGCGAAAGTAATAAGCAACAGGCAGTGAATAAGTTGACAGAGATTGTCGCTAACTTTACCGCCATGATTTCTACCCGAATGCCTGATGACGTGGTGGATAAACTAAAACAGCTAA  >  NC_000913/3206293‑3206583

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