Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,487,163	0	C	T	47.5%	48.6 / 42.8	60	intergenic (+200/‑72)	hrpA/ydcF	putative ATP‑dependent helicase/DUF218 superfamily protein, SAM‑binding
Reads supporting (aligned to +/- strand):  ref base C (7/24);  new base T (28/0);  total (36/24)
Fisher's exact test for biased strand distribution p-value = 1.23e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGG  >  NC_000913/1487016‑1487313                                                                                                                                                    |                                                                                                                                                            tttttATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTtcct                                                                                                                                                           <  1:405660/150‑1 (MQ=255)         aaTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCaga                                                                                                                                                  <  1:580218/151‑1 (MQ=255)          aTTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCagag                                                                                                                                                 <  2:421562/151‑1 (MQ=255)            tCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGt                                                                                                                                                <  2:684464/150‑1 (MQ=255)                aaaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATACGATTTTCTTCTTTTTGAGT‑GCTcc                                                                                                                                          >  1:351286/1‑151 (MQ=255)                 aaCGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCg                                                                                                                                         <  1:575732/151‑1 (MQ=255)                   cGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGATTAAGAAAAACTAATTCCTCCcttttaagagct                                                                                                                                             >  2:933655/1‑137 (MQ=255)                      tAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGGTTAGGAAAAACAACTTACTCCTTTCCGAGT‑g                                                                                                                                              >  1:72590/1‑141 (MQ=255)                       aaGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATAAAATTTCCTCCttttaatgggggccgc                                                                                                                                         >  1:37051/1‑133 (MQ=255)                       aaGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTtct                                                                                                                                    <  1:559494/150‑1 (MQ=255)                            tAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACAAACGATTTTCTCCTTTTAGAGT‑GCTCCCCctct                                                                                                                                    >  2:490021/1‑145 (MQ=255)                             aaaaaTTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATACAAATTTCTCCCTTctaagggg                                                                                                                                              >  1:798106/1‑128 (MQ=255)                                        aaTTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAAAATTTCTCCTTTTACAGT‑GCCCCCCCTCTCACCCTTAtttcccccc                                                                                                                   >  2:59275/1‑145 (MQ=255)                                           tCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAATAAGATTTTCTTCCTTCTTAGG‑GGTCCGCTTCTCCCTTTTTTCTCACCCAg                                                                                                                  >  1:772860/1‑148 (MQ=255)                                            cAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑tt                                                                                                             <  2:723450/151‑1 (MQ=255)                                             aGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGACAAAACATTTTCTCCCTTTAGAGT‑GCTCCCCCTTTCCCCATTATTTCCCCCCGTA‑tt                                                                                                             >  2:719298/1‑150 (MQ=255)                                               aaaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAGATTTTCTCCTTTTTGAGG‑GCTCCGCTTCTCCCTATTAtttcctcccgtattttt                                                                                                           >  2:94199/1‑138 (MQ=255)                                               aaaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑tt                                                                                                         <  1:540280/151‑1 (MQ=255)                                               aaaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAGAATTTCTCCTTTCATAGT‑GCTCCCCTTCTTACTATTTtcttccaccgcattttt                                                                                                           >  2:470949/1‑137 (MQ=255)                                                aaaaTACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑a                                                                                                       <  2:169274/151‑1 (MQ=255)                                                             ctgctgCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGACATAAAATTTTCTCCTTTCCGAGG‑GCTCCCCCTCTCACTATTTTTTCCCCCCGTATTTT‑TT‑AAGGGAA‑GATaa                                                                                           >  2:762100/1‑150 (MQ=255)                                                                  gcAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAggagga                                                                                     <  2:373277/151‑1 (MQ=255)                                                                     aGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGCTTAAGAAATACCATTTTCTCCCTTTCGAGG‑GGTCCCCTTCTCACTCTTTTTTCACCCAGGA‑TTT‑TT‑TAAGGAACGATAaagagggaac                                                                                  >  2:811155/1‑147 (MQ=255)                                                                       aTGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCAt                                                                                <  1:245033/151‑1 (MQ=255)                                                                       aTGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGGTTAAGACAATAGAATTTCTCCTTTTCGAGG‑GGTCCCCTTCTCCCCCTTATTTTCCGCAGTA‑TTT‑TTAAAGGGaagaaaa                                                                                            >  2:975574/1‑135 (MQ=255)                                                                         gAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑a                                                                             <  1:358729/151‑1 (MQ=255)                                                                         gAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑a                                                                             <  2:392900/151‑1 (MQ=255)                                                                             aatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACa                                                                          <  1:227940/150‑1 (MQ=255)                                                                             aatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACAt                                                                         <  2:901037/151‑1 (MQ=255)                                                                             aatGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAGATTTTCTCCCTTCCGAGG‑GGTCCGCCTCTCCCTATTATTTTACCCAGGA‑TTC‑TT‑TAAGGAAAGATAAAGAGGGAACATG‑aaaa                                                                          >  2:95700/1‑148 (MQ=255)                                                                              atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACAtt                                                                        <  1:335496/151‑1 (MQ=255)                                                                              atGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATAAGAATTTCTCCTTTCTTTGG‑GCTCCGCTTCTCCCTTTTATTTCACGCCGTA‑TTT‑TT‑TAGGGGAAGATAAAGAGGAACCATG‑AACa                                                                          >  1:805499/1‑149 (MQ=255)                                                                                 tttATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGAAATAAGATTTTCTCCTTTTTGAGG‑GCTCCGCTTCTCACTATTTTTTCACCCAGTA‑TTCTTT‑AAGGGAA‑GAAAAAGAGGaaaccagaaaatt                                                                         >  2:11642/1‑137 (MQ=255)                                                                                   tATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATAAGATTTTCTCCCTTTAGAGT‑GGTCCGCCTCTTCCTATTTTCTCACGCCGTA‑TTT‑TT‑AAAGGGACGGTAAAGAGGAACCATT‑AAAAATAcccc                                                                   >  2:26876/1‑151 (MQ=255)                                                                                    aTCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTAcccc                                                                   <  1:146070/150‑1 (MQ=255)                                                                                    aTCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTAccc                                                                    <  2:536930/149‑1 (MQ=255)                                                                                     tCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAAAAAATTTCTCCTTTTTGAGT‑GCTCCCCCTCTCCCTATTATCTCCCACAGGA‑TTT‑TT‑TATGGGACGATAAGGAGGGAACATA‑AAAATTAacccc                                                                  >  1:338773/1‑149 (MQ=255)                                                                                      cTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATAAGGTTTTCTCCCTTTACAGG‑GCTTCCCCTTTCACTATTTTTTCCCCCAGGA‑TTCTTT‑AAGGGAA‑GAAAAAGAGGGAACAAG‑AAAAATAcccc                                                                   >  2:48449/1‑148 (MQ=255)                                                                                       tACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAATACGATTTTTTCCTTTCTGAGT‑GCTCCGCCTTTCCCTATTTTCTCACACAGGA‑TTC‑TT‑TAAGGAACGAAAAAGAGGAAACATG‑AAAATTAaccccctt                                                               >  2:393761/1‑147 (MQ=255)                                                                                        aCAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGttt                                                               <  2:225907/150‑1 (MQ=255)                                                                                         cAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGCGTTGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCCCTATTATCTCACGCAGTA‑TTC‑TT‑AAGGTAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTcc                                                             <  1:284562/151‑1 (MQ=255)                                                                                            cATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAATGACAAAAAAATTTTTCCTTTTCTTGG‑GCTCCGCCTCTCCCTATTATTTCACCCAGGA‑TTT‑TT‑TAGGGGAAGATAAAGAGGGAACATGAAAAATAACCCCGtt                                                                >  1:64231/1‑146 (MQ=255)                                                                                              tttCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATAAGAATTTCTCCCTTTAGAGT‑GGTCCGCTTCTTCCTTTTATTTCACCCCGGA‑TTC‑TT‑TAGGGGACGAGAAGGAGGGAACATG‑AACATTAACCCCCTTTCGAc                                                          >  2:286551/1‑149 (MQ=255)                                                                                               ttCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGAAATACAAATTTTTTCTTTTTTAGG‑GGTCCGCTTTTCTCTATTTTTTTACCCCGTA‑TTC‑TT‑TAGGGGACGATAAGGAGGGAACATG‑AAAATTAACCCGCTTCcgccgc                                                        >  1:221385/1‑150 (MQ=255)                                                                                               ttCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGATTAGGAGCTAAGATTTTCCCCCTTTTTATT‑GGTCCGCTTTTCCATATTATTTCACTCAGGA‑TTT‑TT‑ATAGGAAAGGAAAAGAGGAAACAAG‑AACAAAAACCCCTTTTCGAcc                                                         >  1:682629/1‑148 (MQ=255)                                                                                                tCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAAGAAAAACGATTTTCTCCCTTTACAGTGGCTCC‑CTTCTCCCTATTTTTTCACGCAGGATTTT‑TT‑AAGGGAA‑GATAAAGAGGGACCATG‑AAAATTAACCCCCTTTCGACGcct                                                      >  2:342010/1‑149 (MQ=255)                                                                                                 ccTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTTAAGAAAAAAGAATTTCTCCCTTTCTAGT‑GCTCCCCTTCTCCCTATTTTTTCCCGCCGGA‑TTT‑TT‑TAAGGGACGAAAAAGAGGGAACATG‑AACAATACCCCCCTTTCGACGCctt                                                     >  2:755112/1‑151 (MQ=255)                                                                                                     aaaaGATATGTCAGGCTTGCGGAGTGGCGGTTAAGAAGAAAAGATTTTCTCCCTTTAGAGA‑GGCCCGCTTCCCACTATTATCTTACCCCGCA‑ATT‑TT‑TAAGGAAAGACAAAGAGGGAACCCG‑AAAACTACCCCCCTTCCGGCGCttttccc                                                 >  2:369781/1‑147 (MQ=255)                                                                                                        aGATATGTCAGGCTTGCGGAGTGGCGGTTTAGGAAATACGAAATTTTCCCTTTTTAGG‑GGTCCGCTTCTCCCTATTATTTCACGCAGGA‑TTT‑TT‑TAGGGGAAGATAAGGAGGAAACAAG‑AACATTACCCCCTTTTCCACGCCTTCTCCgg                                               >  1:68992/1‑150 (MQ=255)                                                                                                              gTCAGGCTTGCGGAGTGGCGGTTTAGGACATAAAAATTTCTCCTTTCTGAGG‑GGTCCGCCTCTCCCTATTTTTTCACGCAGTA‑TTT‑TT‑TAGGGGAAGGTAAAGAGGAACCATG‑AAAATTACCCCCTTTTCCACGCTTTTGCCGGCCAcca                                         >  1:150911/1‑148 (MQ=255)                                                                                                              gTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACtat                                        <  2:530063/151‑1 (MQ=255)                                                                                                              gTCAGGCTTGCGGAGTGGCGGTTAAGGAAAAAACAGTTTCTCCTTTCAAAGG‑GGTCCGCCTTTCCCTTTTATTTTACCCACTA‑GTC‑TT‑ATTGGGGAGAAGAAGAGGGAACACA‑AACAATAACCCCTTTCCGCCGCCTTCGCCgggca                                            >  1:26289/1‑144 (MQ=255)                                                                                                                       gCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCAt                               <  2:93927/151‑1 (MQ=255)                                                                                                                              ggCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGt                         <  2:751715/150‑1 (MQ=255)                                                                                                                              ggCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGt                         <  2:409322/150‑1 (MQ=255)                                                                                                                               gCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGtt                        >  2:646802/1‑150 (MQ=255)                                                                                                                                cGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTAt                      >  1:172867/1‑151 (MQ=255)                                                                                                                                 ggTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTAt                      <  2:171794/150‑1 (MQ=255)                                                                                                                                 ggTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATc                     >  1:645577/1‑151 (MQ=255)                                                                                                                                 ggTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATc                     >  1:565197/1‑151 (MQ=255)                                                                                                                                               gATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGtggctgg        >  1:320874/1‑150 (MQ=255)                                                                                                                                                    cctcctTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAgg  >  2:293559/1‑151 (MQ=255)                                                                                                                                                    cctcctTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGAAg   <  2:615266/150‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                            TTTTTATTAATTCTAAAACGGTAAGGGTAAAAATTCAGGAATTCAGAAAAATACAATTCTCTGCTGCAAGATGAATAATGTTTATCTACAGCATTTCCTTAAAAGATATGTCAGGCTTGCGGAGTGGCGGTTAAGGACATACGATTTCCTCCTTTCAGAGT‑GCTCCGCTTCTCACTATTATCTCACGCAGTA‑TTC‑TT‑AAGGGAACGATAAGGAGGAACCATG‑AACATTACCCCGTTTCCGACGCTTTCGCCGGCAACTATAGATGCCATAAATGTTATCGGACAGTGGCTGGCGCAGG  >  NC_000913/1487016‑1487313

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