Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,080,837	0	G	T	41.7%	61.0 / 29.2	60	intergenic (+24/‑519)	putP/efeU	proline:sodium symporter/ferrous iron permease (pseudogene)
Reads supporting (aligned to +/- strand):  ref base G (10/25);  new base T (25/0);  total (35/25)
Fisher's exact test for biased strand distribution p-value = 3.83e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAC  >  NC_000913/1080688‑1080987                                                                                                                                                       |                                                                                                                                                         ggCAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑Ttttgt                                                                                                                                                         <  2:743471/151‑1 (MQ=255) ggCAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTGGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGCCTTAGCCTGCGGCGG‑Ttttgt                                                                                                                                                         <  2:774948/151‑1 (MQ=255)  gCAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑Ttttttt                                                                                                                                                        >  1:1032876/1‑151 (MQ=255)  gCAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑Ttttgtt                                                                                                                                                        <  2:865622/151‑1 (MQ=255)     gTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTgg                                                                                                                                                     <  1:190468/151‑1 (MQ=255)       aTTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑Cttttttt                                                                                                                                                       >  2:389489/1‑147 (MQ=255)         tGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCtt                                                                                                                                                  <  2:656175/150‑1 (MQ=255)         tGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCtt                                                                                                                                                  <  1:572170/150‑1 (MQ=255)           ggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑tttttttttggttt                                                                                                                                                 >  1:593567/1‑143 (MQ=255)           ggATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGGTTccg                                                                                                                                               >  1:882757/1‑149 (MQ=255)             aTTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTcagc                                                                                                                                              >  1:1092053/1‑150 (MQ=255)                gTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTGGGttttgcaac                                                                                                                                           >  2:1132383/1‑143 (MQ=255)                gTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCg                                                                                                                                          <  2:106890/151‑1 (MQ=255)                 tAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑GTTTGTTTTGCTTTAGCAGCgg                                                                                                                                         >  2:27077/1‑151 (MQ=255)                   gtgtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGGTTCTGCAGCggg                                                                                                                                        >  1:456072/1‑150 (MQ=255)                   gtgtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGt                                                                                                                                       <  1:1184185/151‑1 (MQ=255)                    tgtTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGGTTTTTTTTTGGCTT‑AACAGCggggt                                                                                                                                      >  1:116984/1‑149 (MQ=255)                      tttAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTGGTTTcttaaacggg                                                                                                                                        >  2:1050654/1‑138 (MQ=255)                      tttAGTTTGCTGGGTAAAGAGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCCCCGCCGTCACGTTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑AGTTGCTTGCCTTTCGCACCgg                                                                                                                                         >  2:412060/1‑146 (MQ=255)                            ttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTGGCTTCCGCAGCGGGGTTCGCT‑ccc                                                                                                                             >  1:1041215/1‑151 (MQ=255)                            ttGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑cc                                                                                                                              <  2:121963/150‑1 (MQ=255)                             tGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑ccc                                                                                                                             <  1:218645/150‑1 (MQ=255)                               cTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑CATTGTTTGGGTTCAGCAGAGGGGTGCGCC‑ACCtt                                                                                                                           >  1:285752/1‑150 (MQ=255)                                tGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑CTTTTTTTGGCTTTAGCAACggggt                                                                                                                                      >  2:745769/1‑137 (MQ=255)                                   gTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTTGGCTCCGCGACGGGGTGCGCG‑CCCCTAAtg                                                                                                                       >  1:1182959/1‑150 (MQ=255)                                     aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAAtgtg                                                                                                                     <  1:654436/150‑1 (MQ=255)                                     aaaGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAAtgtg                                                                                                                     <  2:929055/150‑1 (MQ=255)                                       aGCGCCGTCCGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑CTTTTTTTTGCTTCCCACCCGGGATGTGCT‑CCCTTAAtgtcgc                                                                                                                   >  2:321934/1‑147 (MQ=255)                                       aGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCt                                                                                                                  <  1:173684/151‑1 (MQ=255)                                           ccGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑CTTTTTTTTGCGTTTGCAGCGGGGTTCGCTCCCCTTAAttggcccccc                                                                                                               >  2:430435/1‑142 (MQ=255)                                              tCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTAAGCCTGCGGCGG‑TTTTTTTTGGCCTCCGCTGCGCGATGGGCA‑TCCCTGATGTGGCTCAccc                                                                                                             >  2:343893/1‑148 (MQ=255)                                               cAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTCAGCAGCGGGGTGCGCT‑CCCTTTATTTGCCTCGcc                                                                                                              >  1:858530/1‑147 (MQ=255)                                                aGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGTTTCAGCGGCGGGGTTTGCT‑CCCTTAATGTGGCTCGCCac                                                                                                            >  1:792656/1‑147 (MQ=255)                                                aGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTCAGCAGCGGGGTGCGCT‑CCCCTAATGTGCCTCGCCctat                                                                                                          >  1:761528/1‑150 (MQ=255)                                                   ggcgATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTTAGCAGCGGGGTGCGCT‑CCCTTAATTTGCCTCGCCCTATaaa                                                                                                       >  1:69626/1‑150 (MQ=255)                                                           aaaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTCAGCCGCGGGGTGCGCT‑CCCCTAATGTGCCCCCCCCTATAAAATGAAttg                                                                                               >  2:667100/1‑148 (MQ=255)                                                            aaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTg                                                                                             <  2:354576/151‑1 (MQ=255)                                                            aaaaCGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTg                                                                                             <  2:493240/151‑1 (MQ=255)                                                                 gCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTCAGCAACGGGGTTCGCC‑CCCTTAATGTGCCCCCCCATATAAAATTAAAGGGGCAggg                                                                                        >  1:860159/1‑151 (MQ=255)                                                                 gCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTGGCTTTAGCAACGGGGTGCGCC‑CCCCTAATTTGGCTCGCCCTATAAAATTAATTGGGCCggg                                                                                        >  2:241569/1‑151 (MQ=255)                                                                    ttGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGa                                                                                       <  2:157619/149‑1 (MQ=255)                                                                     tGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAgc                                                                                     <  1:1180946/150‑1 (MQ=255)                                                                         gAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTTAGCAGCGGGGTGTGCC‑CCCCTAATGTGCCTTGCCCTCTTAATTGAATGGTggggggggggcg                                                                                  >  2:1088666/1‑138 (MQ=255)                                                                               gATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGggc                                                                         <  1:113426/151‑1 (MQ=255)                                                                                aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGggcg                                                                        <  2:116977/151‑1 (MQ=255)                                                                                aTGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGggcg                                                                        <  1:277578/151‑1 (MQ=255)                                                                                  gcgcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGggcggc                                                                      <  2:1031080/151‑1 (MQ=255)                                                                                   cgcACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGCTTCAGCAGCGGGGTGCGCT‑CCCTTTATGTGCCCCGCCCTATAAATTGAATTGGGCAGGGGGGGCG‑CAgggggggg                                                                       >  1:19973/1‑149 (MQ=255)                                                                                        tATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGGTGCGCC‑CCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGGGGGCCATcg                                                                >  1:316104/1‑151 (MQ=255)                                                                                        tATCATTCGGCTCCGCCGTAACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTTGCTTCCGCAACGGGGTGCGCT‑CCCTTAATGTGCCTCGCCCTATAAATTGAATGGTGCAGGGAGCGCG‑CCGGGGGCGGCCaaac                                                                 >  1:552834/1‑148 (MQ=255)                                                                                               cGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTTGCTTCAGCAGCGGGGTGCGCC‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGGGGCCAATcgccgcc                                                           >  1:341739/1‑149 (MQ=255)                                                                                                    ccgccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGccccc                                                     <  2:1032807/150‑1 (MQ=255)                                                                                                      gccgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTGGGTTTAGCAGCGGGGTGCGCT‑CCCCTAATGTGGCTCGCCCTATAAATTGAATTGTGGAGGGGGCGCGCCAGGGGGGGGCC‑ATCGCcccccccccc                                                     >  2:1060664/1‑148 (MQ=255)                                                                                                        cgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTTGCTTCAGCAGCGGGTTGCGCT‑CCCTTTATGTGCCTCGCCATATAAATTGAATGGTGCAGGGGGCGCG‑CAGGGGGGGGCCAATCGCCGcccccccc                                                     >  1:680394/1‑146 (MQ=255)                                                                                                        cgTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCctgctg                                                <  2:1160303/151‑1 (MQ=255)                                                                                                            aCGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTTTTTGGCTTTAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCTCCCTATAAAATGAATGGTGCGGGGGGAGCG‑CAGGGGGCGGCCAATCGCCCCCGCCCCCCGCTGTccc                                            >  1:263159/1‑151 (MQ=255)                                                                                                                  gCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCt                                       <  2:947930/150‑1 (MQ=255)                                                                                                                  gCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCt                                       <  2:1129661/150‑1 (MQ=255)                                                                                                                    aGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTTGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTTAATGGGGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCCCCCCCCGCTGTCCCCGcc                                        >  1:799011/1‑147 (MQ=39)                                                                                                                      gAAAGCTAAGGGACTTAGCCTGCGGCGG‑GTTTTTTTTGTCGCAGCAGCGGGGTTCGGT‑ACCCTTATGTGTCTCGCCCCATAAATTTAAATGTTGTGGGAGCGCG‑TAGGGGGGGGGGAATCACCCCCCCCCCCTGCGGACCCCGGCTTCggg                                  >  1:396027/1‑151 (MQ=255)                                                                                                                        aaGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGa                                <  2:901973/151‑1 (MQ=255)                                                                                                                            tAAGGGACTTAGCCTGCGGCGG‑TTTTTGTTGGCTTCAGCCGAGGGGTGCGCC‑CCCCTTATGTGCCTCCCCCTATAAAATTAATGGGGGAGGGGGGGGG‑CAGGGGGGGGGCCATACCCCCCCCCCCCCGCCCTCCCGGccctttggg                                 >  2:857611/1‑139 (MQ=255)                                                                                                                             aaGGGACTTAGCCTGCGGCGG‑TTTTTTTTTGGTTCCGCAGCGGGGTGCGCT‑CCCCTTATGTGCCTCGCCATATAAAATTAATTGTGCAGGGAGGGCG‑CCGGGGGGGGCCCATTCCCGCCGCCCCCCGCCGTCCCCGCCCTCcggg                                 >  1:180610/1‑145 (MQ=255)                                                                                                                                           gcggcgG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTg               <  2:606283/149‑1 (MQ=255)                                                                                                                                            cggcgg‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACg            <  1:1189647/151‑1 (MQ=255)                                                                                                                                             ggcgg‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGc           >  1:1021552/1‑151 (MQ=255)                                                                                                                                                    tttgtttgGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGccacca      >  1:465780/1‑150 (MQ=255)                                                                                                                                                      tgtttgGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAac    <  2:724476/150‑1 (MQ=255)                                                                                                                                                       gtttgGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAacac  <  2:144456/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                         GGCAGTATTGGGATTGTAGTGTTTAGTTTGCTGGGTAAAGCGCCGTCAGCGGCGATGCAAAAACGCTTTGCCGAGGCCGATGCGCACTATCATTCGGCTCCGCCGTCACGGTTGCAGGAAAGCTAAGGGACTTAGCCTGCGGCGG‑TTTTGTTTGGCTTCAGCAGCGGGTTGCGCT‑CCCTTAATGTGCCTCGCCATATAAATTGAATGGTGCAGGGAGCGCG‑CAGGGGGCGGCCAATCGCCGCCGCCCCCTGCTGTCCCGGCCTTCGGGGAACGCTTCAGCGATTTTGACGCCACCAACAC  >  NC_000913/1080688‑1080987

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