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

GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925751‑926047                                                                                                                                                       |                                                                                                                                                        gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCGCCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGAt                                                                                                                                                        <  2:49108/151‑1 (MQ=255) gAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGa                                                                                                                                                         <  2:416073/150‑1 (MQ=255)   gTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTc                                                                                                                                                      <  1:1027381/151‑1 (MQ=255)         tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTCGa‑a                                                                                                                                                  >  1:551795/1‑148 (MQ=255)         tGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑Ac                                                                                                                                                 <  2:551781/149‑1 (MQ=255)            cGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGC‑Accccc                                                                                                                                             >  2:650889/1‑150 (MQ=255)              tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTGGA‑Accccc                                                                                                                                             >  2:375820/1‑149 (MQ=255)              tACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGA‑Aaccccc                                                                                                                                            >  2:260855/1‑148 (MQ=255)               aCAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑a                                                                                                                                        <  1:1307477/151‑1 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTTA‑GGGGGGGGTTTGT‑Acacccccaa                                                                                                                                         >  1:1364191/1‑146 (MQ=38)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGA‑Acccccccta                                                                                                                                         >  1:1242135/1‑146 (MQ=255)                cAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑Ac                                                                                                                                       <  2:1420597/151‑1 (MQ=255)                 aGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACg                                                                                                                                      <  2:1121036/151‑1 (MQ=255)                           gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGtata                                                                                                                            <  2:479937/151‑1 (MQ=255)                           gAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GAGGGGGGGTCGA‑Accccccgcaccttggccttt                                                                                                                              >  2:459349/1‑135 (MQ=37)                              aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATacac                                                                                                                         <  1:648940/151‑1 (MQ=255)                              aTTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATacac                                                                                                                         <  1:419828/151‑1 (MQ=255)                                tCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGTTTCGA‑AACCCCGTA‑AAGTTGTCGTCTacac                                                                                                                         >  1:236285/1‑149 (MQ=25)                                   caGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGC‑ACCCCCGCT‑ACGTTGCCGTATAAAC‑ACttt                                                                                                                   >  1:519601/1‑151 (MQ=255)                                     ggTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTcc                                                                                                                 <  2:1134491/151‑1 (MQ=255)                                       tAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGA‑ACCCCCCAA‑ACGGTGGCGTTTACAA‑ACctttc                                                                                                                 >  1:412988/1‑147 (MQ=25)                                        aaCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAgg                                                                                                              <  2:126909/151‑1 (MQ=255)                                         aCATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑ACCCCCCCT‑ACGTTGCCGTATACAC‑ACTTTcccg                                                                                                               >  2:704070/1‑147 (MQ=255)                                          cATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCTA‑ACCCCCCCT‑ACGTTGCCGTATAAAC‑ACTTTCCCgg                                                                                                              >  2:433466/1‑149 (MQ=21)                                              cTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑ACCCCCCCT‑ACGTTGCCGTATACAC‑ACCTTCCCgg                                                                                                              >  2:1193227/1‑145 (MQ=21)                                               tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑AACCCCCAT‑ACGTTGCCGGATACAC‑ACTTTCCAGGCGGGCTc                                                                                                       >  1:456199/1‑151 (MQ=21)                                               tCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGA‑ACCCCCCCT‑ACGTTGCCGTATACAC‑ACCTTCCAGGCggg                                                                                                          >  2:134907/1‑146 (MQ=255)                                                ccACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGG‑GGGGGGGGATCGA‑AACCCCCGC‑ACGTTGCCGTTTATAC‑ACTTTTCAGGCGTGc                                                                                                         >  1:544846/1‑148 (MQ=21)                                                   cccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTc                                                                                                   <  1:161222/151‑1 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGA‑ACCCCCCCT‑AAGTTTGCGTATACAC‑ACTTTTCCGGCg                                                                                                            >  2:1359324/1‑141 (MQ=21)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCa                                                                                                  <  2:1143344/151‑1 (MQ=255)                                                    ccGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCa                                                                                                  <  2:35924/151‑1 (MQ=255)                                                     cGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGA‑ACCCCCCAT‑ACGTTGCCGTATACAC‑ACTTTTCCGGCGTGGTCCtt                                                                                                    >  2:337482/1‑148 (MQ=18)                                                         caacaTGTTCACTGCCGTTCAGGCAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCAc                                                                                             <  2:719634/148‑1 (MQ=255)                                                           ccaTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGA‑AACCCCCCT‑ACGTTTCCGTATACAC‑ACTTTTCCGGCGGGGTCCTTCCGCCAccc                                                                                           >  2:590049/1‑149 (MQ=18)                                                             aTGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCgg                                                                                         <  2:1242107/151‑1 (MQ=255)                                                              tGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGG‑AACCCCCCT‑ACCGTGCCGGATACAC‑ACTTTTCTGGCGTGCTCCTTCAACCACCCgg                                                                                         >  1:42725/1‑150 (MQ=18)                                                                 tCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGA‑ACCCCCCAT‑ACGTTTTCGTATTAAC‑ACATTTCCGGCGGGCTCCCTCACCCCCCCCGaca                                                                                      >  2:667975/1‑150 (MQ=17)                                                                      gCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCAc                                                                               <  1:225018/151‑1 (MQ=255)                                                                       ccGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑ACCCCCCAT‑ACGTTTCCGTATAAAC‑ACCTTTCAGGCGTGCTCCTTCCGCCACCCCGGCA‑ccccccc                                                                              >  1:767232/1‑146 (MQ=17)                                                                          tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑AACCCCCAT‑ACGTTGCCCTATACAC‑ACTTTTCCGGCGTGCTCCTTCAGCCCCCCCGGCA‑CCCCACCaa                                                                            >  1:497034/1‑150 (MQ=17)                                                                          tACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGG‑ACCCCCGAT‑ACGGTTCCGTATAAAC‑ACCTTTCACGCGTGCGCCCTCAGCCACTCCGACA‑cccccccca                                                                            >  1:997472/1‑143 (MQ=17)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCaaa                                                                           <  2:835710/150‑1 (MQ=255)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAAt                                                                          <  1:158180/151‑1 (MQ=255)                                                                           acagacagATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAAt                                                                          <  1:843671/151‑1 (MQ=255)                                                                            cagacagaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGAGA‑GGGGGGGGGTTGA‑ACCCCCCCT‑AAGTTTCCGTATAAAC‑ACTTTTCAGGCGTGCTCCTTCCGCCACCCCGGCA‑CCCCCCCaaa                                                                           >  2:111731/1‑149 (MQ=14)                                                                                  gaTAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTg                                                                   <  2:511549/151‑1 (MQ=255)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGA‑AACCCCCAC‑AAGGTGGCGCATACAC‑ACCTTTCCGGGGGGCTCCTTTCGCCACCCGGGCA‑CCCCCCCAAAATGTTTTTCt                                                                 >  2:107823/1‑150 (MQ=9)                                                                                     aaaaTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTCGA‑AACCCCCAT‑ACGTTGCCGTATACAC‑ACCTTCCATGCGGGCTCCTTCAGCCACCCGGGCA‑CCCCACCCAATttttttttt                                                                 >  2:627414/1‑147 (MQ=14)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGA‑ACCCCCCAT‑ACGTTTCCGTATACAC‑ACTTTCCCAGCGTGCTCCTTCAGCCACTCGGGCA‑CCCCCCCAAAAtttttt                                                                    >  1:1090799/1‑144 (MQ=11)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGG‑ACCCCCCATAACGTT‑TCGTATAAAC‑ACTTTTCAGGCGTGGTCCTTTAGCCCCCCCGGCA‑CCCCACCCAAATGttttttttcc                                                              >  2:1367617/1‑144 (MQ=9)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCt                                                                 <  2:1090776/147‑1 (MQ=40)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCa                                                             <  1:1071182/151‑1 (MQ=255)                                                                                        aTGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCa                                                             <  2:475793/151‑1 (MQ=255)                                                                                          gCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGAGGGGGGGGTTTTGA‑ACCCCC‑AC‑ACGTTGCCGTATACAC‑ACTTTTCCGGCGGGCTCCTTCAGCCCCCCCGGCA‑CCCCCCCCAAATGttttttttac                                                              >  1:1136045/1‑142 (MQ=2)                                                                                             aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCt                                                        <  2:1319936/151‑1 (MQ=255)                                                                                             aaaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGTTTGG‑ACCCCCGAA‑AAGTTGGCGTATATAC‑ACATTTCAAGCGGGGCCCTTTCCCCCCTCGGGCA‑CCCCACCCAAATTTTTTTCTTCCAAAccc                                                        >  1:100520/1‑150 (MQ=9)                                                                                              aaaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTc                                                       <  1:486881/151‑1 (MQ=255)                                                                                               aaaaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCa                                                      <  2:667183/151‑1 (MQ=255)                                                                                                 aaaaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGATA‑GGGGGGGATTCGG‑AACCCCCCT‑ACGTTTCCGTATACACAACTTTCC‑GGGGGGGTCCTTCAGACCCCCCGACA‑CCCCACCAAAAtttttttttttcaaa                                                           >  2:62586/1‑135 (MQ=2)                                                                                                   aaGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGG‑ACCCCCCAC‑AAGGTGCCGTATACAC‑ACTTTTCAGGGGGGCTCCTTCAGCCCCCCCGGCA‑CCCCCCCCAATTTTTTTTCTACCCAACCCCACggg                                                  >  2:779544/1‑151 (MQ=9)                                                                                                     gCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTTGA‑AACCCCCCT‑ACGTTTCCGGTTACAC‑ACTTTTCCGGCGTGGTCCTTCACCCACTCGGACA‑CCTCACCAAAATTTTTTTCTACCCAAACACATTGGTg                                                >  1:515621/1‑151 (MQ=9)                                                                                                      cTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGG‑ACCCCCCAT‑ACGTTTCCGTATACAC‑ACTTTTCAGGCGGGCTCCCTCCGCCCCCCCGACA‑CCCCACCCAAATTTTTTTCTACCAAACCCCATTGgggg                                               >  2:1188647/1‑151 (MQ=9)                                                                                                        cGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGGTTCGA‑AACCCCCGT‑AAGTTGCCGCATAAAC‑ACTTTCCCAGCGTGGTCCTTCAACCCCCCCGACA‑CCTCCCCCAAATGTTTTTCTACCAAACCTCCTgggggg                                               >  1:914088/1‑149 (MQ=9)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGG‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACg                                          <  1:1269705/151‑1 (MQ=255)                                                                                                           aCTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTCGA‑AACCCCCAT‑ACGTTTCCGTATACAC‑ACCTTTCAGGCGGGCTCCTTCAGCCCCCCCGACA‑CCCCACCAAATTGTTTTTTTACCCAAACCCATGGGGGGGaa                                            >  1:1394855/1‑149 (MQ=11)                                                                                                              ttCGTACGAGCTCTTCTTTAAATATGGCGGGGA‑GGGGGGGGGTTGA‑AACCCCCAC‑ACGTTTCCCTTTAAAC‑ACTTTTCCGGCGGGGTCCTTCCGCCAACCCGGCA‑CCCCCCCCAAATGTTTTTTTACCCAACCCCCCGGGGGGCAACcggg                                       >  2:742357/1‑151 (MQ=2)                                                                                                                 gTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGGTTGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTTCAGGCGGGCTCCTTCAGCCACCCCGGCACCCCCCCCAAATGTTTTTTTTCCCAAACCCCTTGGG‑GGGAACgggggg                                     >  1:540377/1‑148 (MQ=2)                                                                                                                  tACGAGCTCTTCTTTAAATATGGCGGGGG‑GGGGGGGGTTTGATACCCCGGAA‑GCG‑TGCCGTATAAAC‑ACTTTTCCCGAGGGCTCCCTCCGCCACCCGGGCA‑CCCCCCCAAAATGTTTTTCTACCCAAACTCATAGGTggggaccggggg                                     >  2:1083448/1‑139 (MQ=9)                                                                                                                      aGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑ACCCCCCAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCTCCCACCCGGGCA‑CCTCCCCCAATTGTTTTTCTACCCAACCTCATGGGTGGGCAAGGGgcg                                     >  1:163664/1‑145 (MQ=14)                                                                                                                       gCTCTTCTTTAAATACGGCGGTGA‑GGGGGTGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCTGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGTGCGCTACtat                              <  2:1398873/151‑1 (MQ=255)                                                                                                                          cttcttTAAATATGGCGGTGA‑GGGGGGGGTTTGA‑ACCCCCCCT‑ACCTTGCCCTATACAC‑CCTTTTCCGGCGTGCCCCCTCAGCCACCCGGCCA‑CCCCACCAAATTGTTTTTCTACCAAACCTCCTGGGGGGGAACAGGGGGCTCCTATAg                            >  2:35184/1‑150 (MQ=9)                                                                                                                            tcttTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAg                            <  2:416668/148‑1 (MQ=255)                                                                                                                               ttAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGt                       <  1:140402/150‑1 (MQ=255)                                                                                                                                tAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGtt                      <  1:532908/150‑1 (MQ=255)                                                                                                                                   atatGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTgg                    >  2:1358770/1‑149 (MQ=255)                                                                                                                                   atatGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAg                  <  2:1226553/151‑1 (MQ=255)                                                                                                                                     atGGCGGTGA‑GGGGGGGGTTCGG‑AACCCCGAT‑ACGTTGCCGTATACAC‑AACTTCCCGGGGTGTTCCCTCAGCCACCCGGACC‑CCCCACCAAAATGTTTTGCTAACAAACCTCATGGGGGGGCACGGGGGGCTACTATAGGGAGtt                      >  2:110394/1‑145 (MQ=18)                                                                                                                                      tGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                >  2:171215/1‑150 (MQ=255)                                                                                                                                       ggCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGATCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTaa               <  2:236278/150‑1 (MQ=255)                                                                                                                                       ggCGGTGA‑GGGCGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTa                <  1:691741/149‑1 (MQ=255)                                                                                                                                          ggTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACg           <  1:1061479/151‑1 (MQ=255)                                                                                                                                           gTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACgg          <  2:1390895/151‑1 (MQ=255)                                                                                                                                            tGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         >  1:39480/1‑151 (MQ=255)                                                                                                                                            tGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGt         >  2:742018/1‑151 (MQ=255)                                                                                                                                              a‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTc        >  2:608563/1‑150 (MQ=255)                                                                                                                                                 ggggggATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCa       >  2:991011/1‑149 (MQ=255)                                                                                                                                                   ggggATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaa   <  1:13691/151‑1 (MQ=255)                                                                                                                                                    gggATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCaagaag  <  1:375830/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                        GAGTTCACTGCCGTACAGGCAGCTTAGAAATTCACAGGTAACATACTCCACCCGCCCACCATGTTCACTGCCGTACAGACAGATAAAATGCGAAAAAAAAGCTCGCACTTTCGTACGAGCTCTTCTTTAAATATGGCGGTGA‑GGGGGGGATTCGA‑ACCCCCGAT‑ACGTTGCCGTATACAC‑ACTTTCCAGGCGTGCTCCTTCAGCCACTCGGACA‑CCTCACCAAATTGTTTTGCTACCAAACCTCATGGGTGGCAACGGGGCGCTACTATAGGGAGTTGGAGTAAAACGGTCAAGAAG  >  NC_000913/925751‑926047

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