Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,969,794	0	A	C	41.7%	102.0 / 16.3	62	G397G (GGT→GGG)	tap	methyl‑accepting protein IV
Reads supporting (aligned to +/- strand):  ref base A (28/7);  new base C (0/25);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 7.75e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGC  >  NC_000913/1969644‑1969944                                                                                                                                                         |                                                                                                                                                       gcgGGTCACTGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑Ca                                                                                                                                                        <  2:604135/151‑1 (MQ=255) gcgGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑Ca                                                                                                                                                        <  1:296673/151‑1 (MQ=255)       cACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGcc                                                                                                                                                  <  2:416237/151‑1 (MQ=255)           gAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACt                                                                                                                                               >  2:462689/1‑150 (MQ=255)               tGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGc                                                                                                                                          >  1:816813/1‑151 (MQ=255)                         gATCATGGCCGCTGCGGCGTTATTCACCAGTTTCGCACCCTGCTGGACACGATTGCCTGACTCTTCGATCAGCCCTTTC‑ATCTCTTTTGCTGCCTGCGTGCTACGGCTGGCAAGATTGCGCACTT‑CTCCTCCCACTACCGTAAAAcga                                                                                                                                   >  1:1213922/1‑146 (MQ=255)                             aTGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                            >  2:1187476/1‑151 (MQ=255)                               ggccgcgGCGGCGTTATTCACCTGTTCCGACCCCGGCTGAACACGTTTGACTGCCTCTTCTATCACCCCTTTG‑ATCTCTTTGGCGGCCGGCGCGCTACGCCTGGGAAGTTGGCCCTCTT‑CACCTGCCACTACCGCAAAACCACGCCCCt                                                                                                                            <  2:627158/146‑1 (MQ=255)                               ggTCGGGGGGGGGTTATTCACCAGTTTCAACCCCTGCTGGACACGTTTGCCTACCTCTTCGACCGGCCTTTTG‑ATCTCTTTTGCTGCCTGCGCGCTCCGGCTGGCAAAATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGc                                                                                                                          <  2:110621/151‑1 (MQ=255)                               ggTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGc                                                                                                                          >  1:453141/1‑151 (MQ=255)                                gTCGCGGCGGCTTTATTCCCCAGTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCACCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGct                                                                                                                         <  2:909953/151‑1 (MQ=255)                                    cggcggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTcc                                                                                                                     >  2:294033/1‑151 (MQ=255)                                     ggcggctTTATTCACCATTTTCGACCCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCTTTTG‑ATCTTTTTTGCTGCCGGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCg                                                                                                                    <  1:481416/151‑1 (MQ=255)                                        ggcgTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgcg                                                                                                                 >  1:1035560/1‑151 (MQ=255)                                          ctttATCCACAATTTTCGAACCCTGCGGGACACGATTGACTGCCTCTCCGACCCGCCCTTTG‑ACCTCTTTTGCTGCCTGCGCGCTCCGGCTGGAAAGATTGCGCTCTT‑CTCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGgc                                                                                                                  <  1:1238828/146‑1 (MQ=255)                                            ttATCCACCAGTTTCGACCCCTGCGGGACACGTTTGCCTTACTCTTCGATCAGCCTTTTG‑ATCTTTTTGGCTCCCTGCGCGCTACGCCGGGCAAGTTTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGcgcg                                                                                                               <  2:129153/149‑1 (MQ=255)                                              aTTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGcc                                                                                                           >  1:185366/1‑151 (MQ=255)                                               ttCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCg                                                                                                          >  2:511512/1‑151 (MQ=255)                                                            accccTGCTGAACACGTTTGACTGACTCTTCGATCAGCCTTTTG‑ATCTTTTTTGCGGCCGGCGCGCTACGGCTGGCAAGTTGCCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCg                                                                                              <  2:110582/148‑1 (MQ=255)                                                              cccTGCTGGACACGTTTGACTGACTCTTCGATCACCCTTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCGGGAAAGTTTGCGCACTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                           <  2:677152/151‑1 (MQ=255)                                                              cccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                           >  2:122994/1‑151 (MQ=255)                                                              cccTGCTGGACACGATTGACTGACTCTTCGACCACCCTTTTG‑ATCTTTTTTGCTGCCTGCGCGCTACGGCTGGAAAGTTTGCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAt                                                                                           <  2:762735/151‑1 (MQ=255)                                                               ccTGTTGGACACAATTGCCTGCCTTTCCGATCAGCCTTTTG‑ATCTTTTTTGCTCCCGGCCCGCTACGGTGGGCAAGTTTGCCCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                          <  2:92892/151‑1 (MQ=255)                                                               ccTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCAtt                                                                                          >  2:278345/1‑151 (MQ=255)                                                                             ttGCCTGCCTCTCCGACCA‑CCCTTTGTATCTTTTTTGCTGCCGGCGCCCTCCGGCTGGCAAGATTGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCa                                                                                  <  1:1214882/145‑1 (MQ=255)                                                                               gactgcctCTTCGATCAGCCCTTTG‑ATCTTTTTTGCGGCCTGCGCGCTACGGCGGGCAAGTTTGCGCCCTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTg                                                                          <  2:45250/151‑1 (MQ=38)                                                                                  tgactCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGt                                                                        >  2:1257757/1‑150 (MQ=255)                                                                                      tctTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTgg                                                                    >  1:776415/1‑150 (MQ=255)                                                                                        ttCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGaa                                                                  >  1:551080/1‑150 (MQ=255)                                                                                        ttCGATCACCCCTTTG‑CTCTTTTTTGCTCCCGGCCCGCTCGGGCTGGAAATATTGCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTg                                                                     <  2:280702/147‑1 (MQ=21)                                                                                          cGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGc                                                               >  1:765071/1‑151 (MQ=255)                                                                                          cGACCACCCTTTTG‑ATCTTTTTTGCTGCCTGCGCCCTAGGGGTGGAAAGATTGCGCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAg                                                                <  2:1114072/150‑1 (MQ=25)                                                                                              ccgCCCTTTG‑TTCTTTTTTGCTGCCTGCCCCCTACGGCTGGCAAGTTTGCCCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                                              <  2:999203/146‑1 (MQ=25)                                                                                              cAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAtt                                                           >  1:47063/1‑151 (MQ=255)                                                                                               cgCCTTTTT‑ATCTTTTTCGCCGCCGGCGCCCTACGGCCGGAAAG‑TTGCCCAATTCCACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCa                                                              <  1:695958/146‑1 (MQ=21)                                                                                               aGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGTCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTc                                                          >  1:648177/1‑151 (MQ=255)                                                                                                 ccttttG‑TTCTTTTTTGCTCCCTGCGCCCTCCGGCTGGCAAAATGGCGCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAAt                                                            <  2:1080885/144‑1 (MQ=21)                                                                                                 ccttttG‑ATCTCTTTTGCGCCCGGCGCCCGCCGGCGGGCAATTTTGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCaa                                                             <  1:1023496/143‑1 (MQ=21)                                                                                                         aTCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                  >  1:884864/1‑150 (MQ=255)                                                                                                           ctctTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                               >  2:521070/1‑151 (MQ=255)                                                                                                            tctTTTGCTCCCGGCGCCCTCCGGCGGCAAAGATGGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGAt                                                   <  1:629091/146‑1 (MQ=21)                                                                                                             tttttGCTCCCTGCGCGCTCCGGCGGGAAAGTTTGCGCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                               <  2:692925/148‑1 (MQ=25)                                                                                                              ttttcctcccgccgctcTCCGTCTGGCAGGATTGCCACCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATa                                                  <  2:181830/133‑1 (MQ=21)                                                                                                              ttttGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTg                                            >  2:410930/1‑151 (MQ=255)                                                                                                              ttttGCGGCCGGCGCCCTCGGGCTGGAAAGTTTGCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                             <  2:962200/150‑1 (MQ=25)                                                                                                              ttttCCTCCCTGCGCGCTACGGCTGGAAAGATGCCGCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACg                                               <  2:634406/148‑1 (MQ=25)                                                                                                              tttgGCTCCCGGCGCGCTACGGCGGGAAAGATTGCGCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCt                                             <  2:795536/150‑1 (MQ=37)                                                                                                                   ctgcctgcGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAAt                                       >  1:692719/1‑151 (MQ=255)                                                                                                                     ccctccGCCCTCCGGCGGACAAGTTGGCCCCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGa                                           <  1:205525/144‑1 (MQ=25)                                                                                                                      cctgcGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTc                                    >  1:181940/1‑151 (MQ=255)                                                                                                                         gcgcgcTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                                 >  2:930582/1‑151 (MQ=255)                                                                                                                          cgcgcTACGGCTGGCAAGATTGCGCCCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                                 <  2:747494/150‑1 (MQ=255)                                                                                                                          cgcgcTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGcc                                 >  1:825873/1‑150 (MQ=255)                                                                                                                          cgcgcAACGGCTGGCAGGTTTGCCCATTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCg                                   <  1:1072229/148‑1 (MQ=25)                                                                                                                             gcTGCGGCGGGCAAGTTTGCCCTCTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCa                                <  1:498651/148‑1 (MQ=25)                                                                                                                                 cGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAAtttttt                         >  2:1114420/1‑151 (MQ=255)                                                                                                                                     tGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                      >  1:1145497/1‑150 (MQ=255)                                                                                                                                      ggaaagttggcgcCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                      <  2:204466/140‑1 (MQ=25)                                                                                                                                      ggaaagttggcgcCCTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCg                      <  2:1053049/140‑1 (MQ=25)                                                                                                                                      gccAAGTTTGCCCACTC‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGc                       <  1:1149975/146‑1 (MQ=255)                                                                                                                                       acAAGATTGCGCACTT‑CCCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGa                     <  2:430310/148‑1 (MQ=255)                                                                                                                                          aGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCtgg                >  1:429050/1‑151 (MQ=255)                                                                                                                                                 gcACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGAt         >  2:784153/1‑151 (MQ=255)                                                                                                                                                     tt‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTcc     >  1:665113/1‑151 (MQ=255)                                                                                                                                                         aCCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGc  >  2:376440/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                       GCGGGTCACCGAACTGACAATATCGATCATGGTCGCGGCGGCGTTATTCACCAGTTTCGAACCCTGCTGGACACGATTGACTGACTCTTCGATCAGCCCTTTG‑ATCTCTTTTGCTGCCTGCGCGCTACGGCTGGCAAGATTGCGCACTT‑CACCTGCCACTACCGCAAAACCACGCCCCTGCTCTCCGGCGCGAGCCGCTTCCACTGCCGCATTCAGGGCCAGAATATTGGTCTGGAAAGCAATTCCGTCGATAACGCTGATAATGTCGCCAATTTTTTGCGAGCTGGTGGCGATCTCCTGC  >  NC_000913/1969644‑1969944

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