Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	925,900	0	A	G	45.4%	156.2 / 30.7	98	noncoding (72/88 nt)	serW	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (7/46);  new base G (44/0);  total (52/46)
Fisher's exact test for biased strand distribution p-value = 1.38e-20
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925751‑926049                                                                                                                                                       |                                                                                                                                                          gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGa                                                                                                                                                           <  1:155846/150‑1 (MQ=255)  aGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGtt                                                                                                                                                         >  1:958427/1‑151 (MQ=255)        cTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAcc                                                                                                                                                   >  2:440840/1‑151 (MQ=255)           ccGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAAccccc                                                                                                                                                <  1:1199679/151‑1 (MQ=255)            cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGGAcccccc                                                                                                                                               >  2:256519/1‑150 (MQ=255)              tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCAt                                                                                                                                             >  2:89604/1‑151 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGAt                                                                                                                                             <  2:714362/150‑1 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGA‑GGGGGGGGGTGGAAccccccaac                                                                                                                                            >  1:1254848/1‑147 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAacccccctac                                                                                                                                           >  2:940968/1‑146 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATAc                                                                                                                                           <  1:1574221/151‑1 (MQ=255)                   gCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGtt                                                                                                                                        <  2:266141/151‑1 (MQ=255)                    cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAACCCCCCATACGGTg                                                                                                                                       >  1:48176/1‑151 (MQ=255)                    cAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTg                                                                                                                                       <  1:637831/151‑1 (MQ=255)                      gCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAACCCCCCCTACGGTggc                                                                                                                                     >  2:438065/1‑149 (MQ=37)                       cTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCCTAAGTTTCCg                                                                                                                                    >  2:318856/1‑151 (MQ=37)                          aGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGta                                                                                                                                  <  1:796072/150‑1 (MQ=255)                             aaTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTGAACCCCCCCAACGTTGCCCtataa                                                                                                                               >  2:1471174/1‑149 (MQ=25)                                 cacaGGTAACATACTCCA‑‑‑‑CCCACCATGTTCACTGCCGTACAGACAGATAAAATGCG‑AAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTcc                                                                                                                     <  1:1116264/150‑1 (MQ=18)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGGTTGGAACCCCC‑ATACGTTTGCGTATACACA‑Cttt                                                                                                                       >  2:1488883/1‑151 (MQ=25)                                       tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGAACCCCCCCTACGGTTCCCTATACACA‑CTTTTCa                                                                                                                    >  1:1454807/1‑150 (MQ=25)                                        aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAg                                                                                                                   <  1:1235042/150‑1 (MQ=255)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGGTCGAAACCCCCCTACGGTTCCGGATACACACCTTTcccgcgg                                                                                                                >  2:1339634/1‑146 (MQ=21)                                           aTACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGt                                                                                                               <  2:319801/151‑1 (MQ=255)                                            tACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCATACGTTGCCGGATACACA‑CTTTCCAAGCggg                                                                                                              >  2:825404/1‑149 (MQ=21)                                             aCTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCGG‑GGGGGGGGGTTGCACCCCCCCCAACTTGCCGCATTCACA‑CCTTTCAggagcgg                                                                                                             >  1:130989/1‑146 (MQ=21)                                                ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAACCCCCCACACCTTTCCCGTTAAAAA‑CTTTCCCGGCGGGCTcc                                                                                                          >  1:1098942/1‑151 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTTCCGTATACACA‑CTTTCCCGGCGGGCTCCTTc                                                                                                       >  2:436572/1‑151 (MQ=18)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCa                                                                                                      <  1:659800/151‑1 (MQ=255)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAg                                                                                                     <  2:777436/151‑1 (MQ=255)                                                      gCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCCTACGGTGCCGTATACACA‑CCTTTCAGGCGGGCTCCTTCCGc                                                                                                    >  1:1251033/1‑151 (MQ=18)                                                       cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCATACCTTGCCGTATACACA‑CTTTCCAGGCGCGCTCCTTCAGcc                                                                                                   >  2:48008/1‑151 (MQ=255)                                                       cccACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCATACGTTGCCGTATACACA‑CTTTTCAGGCGGGCTCCTTCTGcc                                                                                                   >  2:148107/1‑151 (MQ=18)                                                        ccaccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACACCCCATACGTTGGCGTATAAACA‑CCTTTCAGGCGTGCTCCTTCAGccc                                                                                                  >  2:1321142/1‑150 (MQ=18)                                                           ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAAACCCCCATACGTTGCCGGATACCAA‑CTTTTCCGGCGGGCTCCCTCCGCCAccc                                                                                               >  2:30789/1‑149 (MQ=18)                                                             aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCgg                                                                                             <  2:470615/151‑1 (MQ=255)                                                                ttCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGGACCCCCCCTAAGGTTGCGTATACACA‑CTTTTCCGGGGGGCTCCCTCAGCCACCCCGac                                                                                           >  2:618998/1‑150 (MQ=17)                                                                gtCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCGTCTTTAAATATGGCGGTGA‑GGGGGCGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGaca                                                                                          <  2:1297474/150‑1 (MQ=255)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCGATACGTTTCCGTATACACA‑CCTTCCAGGCGGGCTCCTTCCGCCACTCCGacac                                                                                         >  1:259301/1‑151 (MQ=17)                                                                      gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCa                                                                                     <  2:596149/150‑1 (MQ=255)                                                                      gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCAc                                                                                    <  1:153005/151‑1 (MQ=255)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGTTGCCGTATACACA‑CTTTTCAGGCGTGCTCCTTCAGCCACCCGGACACCCCAcc                                                                                   >  2:231359/1‑151 (MQ=17)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AAt                                                                             <  1:426591/151‑1 (MQ=255)                                                                            cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTGAAACCCCCCTACGCTGCCGGATAAAAA‑CCTTCCCGGGGGGGTCCCTCAGCCCCCCCGGCACCCCACC‑c‑aaat                                                                            >  1:196752/1‑149 (MQ=9)                                                                                cagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGAAACCCCCATAAGTTGCCGTATACACA‑CTTTTCAGGCGGGCTCCTTCAGCCCCCCGGACAACTCCCC‑A‑AATTGttt                                                                        >  2:1219657/1‑151 (MQ=14)                                                                                 agaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGtttt                                                                       <  1:1258484/151‑1 (MQ=255)                                                                                   aTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGGTGCCGTATACACA‑CTTTTCAGGCGGGCTCCCTCAGCCACCCGGACACCCCACC‑C‑AAATGtttttt                                                                     >  2:357361/1‑149 (MQ=14)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGTTGCCGTATACACA‑CTTTTCAGGCGGGCTCCTTCAGCCACCCCGACACCCCACC‑A‑AAATGTTTTTCt                                                                    >  2:169913/1‑150 (MQ=14)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGG‑GGGGGGGGGTCGAACCCCCCACACGCTGCCGTATACACA‑CCTTTCACGCGGGCCCCCTCACCCACCCCGACACCCCACC‑ACAATttttttttt                                                                    >  1:1434608/1‑148 (MQ=9)                                                                                      aaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGAACCCCCCTACGCTGCCGGATACACA‑CTTTTCCCGGGGGCTCCCTTACCCCCTCGGGCACCCCCCC‑A‑AAATGttttttttc                                                                  >  1:1565264/1‑146 (MQ=9)                                                                                           cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGAACCCCCCATACGTTTCCGTATACCCA‑CCTTTCCCGCGTGCTCCCTCTGCCACCCGGACACCCCCCC‑A‑AAAttttttttttcccaac                                                             >  1:531696/1‑141 (MQ=2)                                                                                           cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGGACCCCCCCTACGTTGCCGTATACACA‑CTTTTCAGGCGGGGTCCTTCAGCCACCCGGACACCCCACC‑A‑AAATTTTTTTTTACCAAA‑c                                                            >  1:1398562/1‑151 (MQ=11)                                                                                           cGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑c                                                            <  1:1538975/151‑1 (MQ=255)                                                                                            gAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑cc                                                           <  2:433168/151‑1 (MQ=255)                                                                                             aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCt                                                          <  1:540609/151‑1 (MQ=255)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGTTGGCGGATACACA‑CTTTTCCGGCGTGCTCCTTCAGCCACCCGGGCACCCCACC‑C‑AAATGTTTTTCTACCAAA‑CCTCCTg                                                      >  2:573175/1‑151 (MQ=14)                                                                                                     gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTg                                                  <  2:537116/151‑1 (MQ=255)                                                                                                       tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGAACCCCCCCTACGTTGCCGTATACACA‑CTTTTCAGGCGTGCTCCTTCCGCCCCTCGGACACCCCACC‑A‑AATTGTTTTTTTTCCAAA‑ACTCCTggggggg                                                >  2:201713/1‑150 (MQ=11)                                                                                                       tCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTTAAACCCCCCTCCGTTGCCGTATAAAAA‑CTTTCCAGGCGGGGTCCTTCCGCCACCCCGGCACCCCACC‑C‑AAATGTTTTTTTTCCCAA‑CCTCATtggtggg                                                >  1:985972/1‑150 (MQ=9)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTgg                                                 <  1:964468/149‑1 (MQ=255)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTgg                                                 <  2:518526/149‑1 (MQ=255)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCa                                               <  2:513565/151‑1 (MQ=255)                                                                                                         gCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGGACCCCCCATACGTTGCCGTATAAACA‑ATTTTCCGGGGGGGTCCTTCCGCCACCCCGGCACCCCACC‑C‑AAATGTTTTTTTTCCAAA‑ACTCATggggggga                                               >  1:388560/1‑148 (MQ=2)                                                                                                          cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTTGAACCCCC‑CTACGTTGTCGTATAAACA‑CTTTCCAAGCGGGGTCCTTCCGCCACCCCGGCAACCCCCC‑A‑AAAATTTTTTTTTCCCAA‑CCACATgggggggcac                                             >  1:1415810/1‑147 (MQ=2)                                                                                                          cACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCCTACGTTTTCGTATTCACA‑CTTTTCAGGCGGGGTCCCTCCGCCACCCCGGCACCCCCCC‑C‑AATTGTTTTTTTTCCAAA‑CCCCCTGGGTGGGAAc                                             >  1:8999/1‑151 (MQ=2)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAACCCCCCATACGTTTCCGTATACACA‑CTTTCCCGGCGGGCCCCTTCAGCCCCCCGGGCACCCCCCC‑C‑AAATTTTTTTCTACCCAA‑ACCCATGGGGGGGAACg                                            >  2:429633/1‑151 (MQ=2)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGAAACCCCCATACGTTGTCGTATACACA‑CTTTCCCGGCGGGCTCCTTCAGCCCCCCCGACACCCCCCC‑C‑AATTGTTTTGTTACCCAA‑CCCCCTGGGTGGCAACg                                            >  1:233380/1‑151 (MQ=14)                                                                                                            cTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAACCCCCCCTACGTTGCCGTATACACA‑CTTTTCAGGCGGGCTCCCTCAGCCACCCGGACACCCCACC‑A‑AAATGTTTTTCTACCAAA‑CCTCCTGGGGGGCAAAgg                                           >  2:1359712/1‑151 (MQ=14)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGgc                                        <  1:1031473/150‑1 (MQ=255)                                                                                                                cGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGgc                                        <  2:140814/150‑1 (MQ=255)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGgcgc                                      <  1:706221/151‑1 (MQ=255)                                                                                                                  tACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCCTACCTTTTCGGTTACACA‑CTTTTCCGGGGGGCTCCTTCCGCCCCCCCGGCACCCCCCCAA‑AATTTTTTTTTTTCCAAA‑CCCCTTgggggggacggggggcc                                      >  1:925752/1‑140 (MQ=2)                                                                                                                   aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGcta                                    <  2:871321/151‑1 (MQ=255)                                                                                                                   aCGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGct                                     <  1:407279/150‑1 (MQ=255)                                                                                                                      aGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGGAACCCCCATAAGTTTCCGGATAAACA‑CTTTTCAGGCGGGCTCCTTCAACCACTCCGAAACCCCACC‑A‑AAATTTTTTTCTTCCAAA‑ACCCATTGGGGGCAACGGGGGGctactt                                 >  2:1153683/1‑150 (MQ=9)                                                                                                                          cttcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACtata                               <  1:282222/149‑1 (MQ=255)                                                                                                                            tcttTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAggg                            <  2:1568298/150‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                           <  1:1525585/148‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGa                           <  1:358886/148‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                         <  1:855726/150‑1 (MQ=255)                                                                                                                                tAAATATGGCGGTGA‑GGGGGGGGGTTGAAACCCCCCCACGTTGTCGTATACACA‑CTTTCCACGCGGGCTCCTTCCGCCCCTCCGGCACCCCCCC‑C‑AAATTTTTTTCTTCCAAACCCCCTTGGGGGGAAACGGGG‑GCGACTATAAGGAGGTg                       >  2:1405455/1‑151 (MQ=2)                                                                                                                                  aaTATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa                     <  1:1608825/151‑1 (MQ=255)                                                                                                                                  aaTATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGa                     <  1:318839/151‑1 (MQ=255)                                                                                                                                     atGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                  <  1:820571/151‑1 (MQ=255)                                                                                                                                      tGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                  >  2:778000/1‑150 (MQ=255)                                                                                                                                       ggCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaaa                <  1:1607771/151‑1 (MQ=255)                                                                                                                                       ggCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                 >  2:789932/1‑150 (MQ=255)                                                                                                                                       ggCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa                 <  1:78111/150‑1 (MQ=255)                                                                                                                                         cGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAAc              >  1:744651/1‑151 (MQ=255)                                                                                                                                          ggTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg             <  2:1099517/151‑1 (MQ=255)                                                                                                                                             gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt           <  2:735530/150‑1 (MQ=255)                                                                                                                                             gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc          >  2:1157301/1‑151 (MQ=255)                                                                                                                                              a‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa         <  1:656079/151‑1 (MQ=255)                                                                                                                                                gggggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaa        <  2:1163162/151‑1 (MQ=255)                                                                                                                                                 ggggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa         >  1:137471/1‑149 (MQ=255)                                                                                                                                                   ggggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa     <  1:606138/151‑1 (MQ=255)                                                                                                                                                    gggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga      >  1:1116816/1‑149 (MQ=255)                                                                                                                                                     ggATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaaga      >  2:954647/1‑148 (MQ=255)                                                                                                                                                      gATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa  <  1:1241167/151‑1 (MQ=255)                                                                                                                                                      gATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAgaagaa  <  1:468887/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                          GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACA‑CTTTCCAGGCGTGCTCCTTCAGCCACTCGGACACCTCACC‑A‑AATTGTTTTGCTACCAAA‑CCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAGAA  >  NC_000913/925751‑926049

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