Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	786,855	0	A	C	45.3%	116.9 / 21.5	75	G247G (GGT→GGG)	gpmA	phosphoglyceromutase 1
Reads supporting (aligned to +/- strand):  ref base A (36/5);  new base C (0/34);  total (36/39)
Fisher's exact test for biased strand distribution p-value = 3.40e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786705‑787005                                                                                                                                                        |                                                                                                                                                       ccGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa                                                                                                                                                        <  1:44765/151‑1 (MQ=255) ccGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa                                                                                                                                                        <  2:823270/151‑1 (MQ=255) ccGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTa                                                                                                                                                        <  2:599427/151‑1 (MQ=255)         gTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGtt                                                                                                                                                >  1:24737/1‑151 (MQ=255)            gCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGc                                                                                                                                             <  1:80009/151‑1 (MQ=255)            gCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGc                                                                                                                                             <  1:181146/151‑1 (MQ=255)             cGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGc                                                                                                                                             >  1:258934/1‑150 (MQ=255)             cGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCa                                                                                                                                            <  2:613626/151‑1 (MQ=255)                aTCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGTCTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAcc                                                                                                                                         <  2:716661/151‑1 (MQ=255)                 tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccg                                                                                                                                        >  1:73522/1‑151 (MQ=255)                 tCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccg                                                                                                                                        >  2:1072901/1‑151 (MQ=255)                  ccGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgc                                                                                                                                       >  2:208620/1‑151 (MQ=255)                   cGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAAccgcc                                                                                                                                      >  1:181432/1‑151 (MQ=255)                         aTCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCgctgct                                                                                                                                >  1:537021/1‑151 (MQ=255)                           ccccccaGGCCTGATAAGTCGCGCGGCGTCGCATCGGGCATTGTGCTCCTTTTTTAGAAAAAAAAAAGCCGCCTCAC‑TTGCAGTCGGTTTTCTCATTTAAAACAAATGCCTTTTCCTTCGCTTTCCCCTGGTTTGCAACCGCCgctg                                                                                                                                  <  2:656487/142‑1 (MQ=255)                            accacAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTc                                                                                                                             >  2:466312/1‑151 (MQ=255)                                 aGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGc                                                                                                                        >  2:822958/1‑151 (MQ=255)                                   gCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCg                                                                                                                       >  1:145819/1‑150 (MQ=255)                                        aaaagcagcgcGGCGTCGCATGAGGCAATGTGCTCAATGGTAAGCAAAAAAAAGGCCGACTCCT‑TTGCAGTCGGCTTTCTTTTTTAAAGAGAATGGGGTTTACTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGAtctc                                                                                                                  <  2:680484/143‑1 (MQ=255)                                          acgTCGCGCGGCGTCGCATCAGGCAATGTGCTCCATTTTTACCAACAAAAAAGCCGCCTCTC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTGCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTc                                                                                                               <  1:741363/149‑1 (MQ=255)                                          aaGCCGCGCGGCTTCGCATCAGGAAATGTCCTCCTTTTTTACAAAAAAAAAACCCGACTCCT‑TGGCGGTCGTTTTTCTCATTTAAAACGAATGACTTTTACTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCg                                                                                                                 <  2:978264/149‑1 (MQ=255)                                            gTCGCGCAGCGTCGCACCAGGCATTGTGCTCCATTGTTAGCAAAAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTAAAACAAATGACTTTTTCTCCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                             <  1:58289/151‑1 (MQ=255)                                                 gcAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAT‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGTTTTCCCCTGTTTTACAACCGCCGCTGCTTTCGCTGCAATCTCGTCa                                                                                                              <  1:542699/145‑1 (MQ=255)                                                   aGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAcc                                                                                                      >  2:677935/1‑151 (MQ=255)                                                    gggTCGCAACAGGCATTGTCCCCCATTTTTAGCAACAAAAAAGCCGACTCTC‑TTGCAGCCGTTTTTCTCTTTTAAAACGAAGGCCTTTTTCTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     <  1:1131344/149‑1 (MQ=255)                                                    gcggcgCATCAGGCATTGTGCTCCATTGTTAACAAAAAAAAAGCCGCCTCCT‑TTGCGGTCGTTTTCCTCATTTAAAACAAATGCCGTTTCCTTCGTTTTCCCCTGGTTTCAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAg                                                                                                             <  2:937991/143‑1 (MQ=255)                                                      gTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTAccc                                                                                                     >  2:1029026/1‑149 (MQ=255)                                                        cGCATCAGGAAATGTGCTCCATTGTAAGAAAAAAAAAAGCCGACTCAT‑TTGCAGTCGGTTTTCTCATTTAAAACAAATGACTTTTACTTCGCTTTCCCCGGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCa                                                                                                    <  1:585246/148‑1 (MQ=255)                                                           aTCAGGAAATGTGCTCCATTGTTAGAAAAAAAAAGGCCGACTCAC‑TTGCAGTCGCCTTTCTCATTTTAAACGAATGCCTTTTACTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataat                                                                                              <  1:200428/151‑1 (MQ=255)                                                               ggaattgtgCTCCATTGTAAGAAAAGAAAAAGCCGACTCAC‑TTGCGTTCGGTTTTCTCTTTTTAAACGAATGCTTTTTACTCCGTTTTCCCTTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataa                                                                                               <  1:652242/141‑1 (MQ=255)                                                               ggCAATGTGCCCCATTGTAAGCAACAAAAAAGCCAACTCAT‑TTGCAGTCGGCTTTCTCATTTAAAACGAAGGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGataata                                                                                             <  1:754841/148‑1 (MQ=255)                                                                  attgtgCTCCATTGTTAGCAAAAAAAAAGCCGACTCCC‑TTGCGGTCGGTTTTCTCATTTAAAACGAATGACGTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGt                                                                                         <  2:956754/147‑1 (MQ=255)                                                                  aaTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGttt                                                                                       >  1:179850/1‑151 (MQ=255)                                                                  aaTGTCCTCCATTTTTAGCAAAAAAAAGGCCGCCTCCC‑TTGCGGTCGGTTTTCTCATTTAAAACGAATGACGTTTCCTCCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGtt                                                                                        <  2:481173/150‑1 (MQ=255)                                                                   aTGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                      >  2:1171371/1‑151 (MQ=255)                                                                     gtgCCCCATTTTTAGAAAAAAAAAAGCCGACCCCC‑TTGCAGTCGGTTTTCTCATTTTAAACGAATGTCGTTTCCTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTc                                                                                      <  2:48440/149‑1 (MQ=255)                                                                          cctttGTTAGCAAAAAAAAAGCCTACTCAC‑TTGCAGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTCCTCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGtt                                                                               <  2:258929/148‑1 (MQ=255)                                                                          caaTTGTTAGAAAAAAAAAACCCG‑CTCCCTTTGCAGGCGTTTTTCTCATTTTAAACGAATTCCTTTTCCTTCGCTTTCCCCTGTTTTGAAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGc                                                                                           <  2:464179/137‑1 (MQ=255)                                                                           cATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGttt                                                                              >  2:804216/1‑151 (MQ=255)                                                                             ttGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTg                                                                             >  2:702091/1‑150 (MQ=255)                                                                              ttttAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGa                                                                            <  1:1127640/148‑1 (MQ=255)                                                                                ttAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGaaa                                                                          >  1:984586/1‑150 (MQ=255)                                                                                  agaaacaaAAAGGCCGACTCAC‑TTGCGGTCGGTTTTCTCATTTTAACGGAATGCCGTTTACTTCGTTTTCCCCTGGTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAAt                                                                         <  1:144727/146‑1 (MQ=255)                                                                                  aGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtc                                                                       >  1:1056612/1‑151 (MQ=255)                                                                                  aGCAAAAAAAACGCCGCCTCTT‑TTGCAGTCGGCTTTCTCTTTTAAAACGAATGACTTTTACTTCGTTTTCCCCGGTTTTGCACCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATtc                                                                       <  1:207835/151‑1 (MQ=255)                                                                                     aaaaaaaaaGCCGACTCCC‑TTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                    <  1:943815/148‑1 (MQ=255)                                                                                     aaaaaaaaaGCCGACTCCC‑TTGCAGTCGGTTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcg                                                                    <  1:222633/148‑1 (MQ=255)                                                                                      acaaAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACTTTTACTTCGCTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   <  1:197621/151‑1 (MQ=255)                                                                                      acaaAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   >  1:891773/1‑151 (MQ=255)                                                                                      aaaaaaaaGCCGACTCCT‑TTGCGGTCGGTTTTCTCATTTTAAACGAATGACGTTTACTTCTTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   <  2:881445/149‑1 (MQ=255)                                                                                      aaaaaaaaGCCCACCCCC‑TTGGCGACGGCTTTCTCATTTTAAACGAATGCCTTTTTCTTCGTTTTCCCCTGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTCCCCAGATAATAGCGTTTCAGCGGTTTGAAATTCtcgt                                                                   <  2:106549/149‑1 (MQ=255)                                                                                                cGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCAt                                                          >  1:282703/1‑150 (MQ=255)                                                                                                  aCTCAT‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCACCATTACCCAGATAATAGCGTTTCAGCGGATTGAAATTCTCGTCGAACTCATac                                                        >  2:518717/1‑150 (MQ=255)                                                                                                    tCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACAcc                                                     >  2:934261/1‑151 (MQ=255)                                                                                                          tGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                >  1:486680/1‑151 (MQ=255)                                                                                                          tGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTCTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCgg                                                >  1:555326/1‑151 (MQ=255)                                                                                                           gCAGTGGGCTTTCTCATTTAAAACGAATGACTTTTCCTTCGTTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGc                                                  <  1:958525/148‑1 (MQ=255)                                                                                                             aGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCAc                                             >  2:188400/1‑151 (MQ=255)                                                                                                               tCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGc                                           >  1:372007/1‑151 (MQ=255)                                                                                                                cgtttttCTCATTTAAAACGAATGTCGTTTTCTTCGCTTTCCCCTGTTTTGAAGCCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCa                                              <  2:825134/143‑1 (MQ=255)                                                                                                                    tttCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGt                                       >  1:420505/1‑150 (MQ=255)                                                                                                                     ttttctttttAAACGAATGCTTTTTCCTTCTCTTTCCCCTGTTTTTAAACCGCGGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTc                                      <  1:1029078/144‑1 (MQ=255)                                                                                                                     ttCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTACGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCg                                     >  2:468547/1‑151 (MQ=255)                                                                                                                         cATTTTAAACGACTGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                                 >  2:757312/1‑151 (MQ=255)                                                                                                                         cATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGat                                 >  1:1018808/1‑151 (MQ=255)                                                                                                                           ttttAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                               >  1:957229/1‑151 (MQ=255)                                                                                                                             taaaaCGAATGACTTTTACTTCTTTTTCCCCTGTTTCACAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGatat                               <  2:1124947/147‑1 (MQ=255)                                                                                                                                       gACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    >  2:141257/1‑150 (MQ=255)                                                                                                                                       gACGTTTACTTCGCTTCCCCTTGTTTTCCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    <  2:3971/150‑1 (MQ=255)                                                                                                                                       aacttttATTTCTCTTTCCCCTGGTTCGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGa                    <  1:672313/146‑1 (MQ=255)                                                                                                                                          gTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAtct                >  1:883265/1‑151 (MQ=255)                                                                                                                                           tttCCTTCCCTTTCCCCTGGTTGGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAAt                  <  1:677959/148‑1 (MQ=255)                                                                                                                                                ttCGTTTTCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct           <  1:1195531/150‑1 (MQ=255)                                                                                                                                                ttCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTcttctt          >  1:530364/1‑151 (MQ=255)                                                                                                                                                 tctcTTTCCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCtcttct           <  2:149777/146‑1 (MQ=255)                                                                                                                                                   gCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGGTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGTATATTAAGCTCAAGAATCTCTTCTTCGc       >  2:558961/1‑151 (MQ=255)                                                                                                                                                    cTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt      >  1:513701/1‑151 (MQ=255)                                                                                                                                                     tttCCCCTGTTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCt      <  2:432906/150‑1 (MQ=255)                                                                                                                                                        aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  2:851902/1‑151 (MQ=255)                                                                                                                                                        aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  1:895532/1‑151 (MQ=255)                                                                                                                                                        aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  1:455795/1‑151 (MQ=255)                                                                                                                                                        aCCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATg  >  2:386879/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                       CCGTCAGAGTGGCGTATCCGATGAATCACCACAGGCCTGATAAGTCGCGCAGCGTCGCATCAGGCAATGTGCTCCATTGTTAGCAACAAAAAAGCCGACTCAC‑TTGCAGTCGGCTTTCTCATTTTAAACGAATGACGTTTACTTCGCTTTACCCTGGTTTGCAACCGCCGCTGCTTTCGCTGCGATCTCGTCAGCATTACCCAGATAATAGCGTTTCAGCGGTTTGAAATTCTCGTCGAACTCATACACCAGCGGCACGCCAGTCGGGATATTAAGCTCAAGAATCTCTTCTTCGCTCATG  >  NC_000913/786705‑787005

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