Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,097,581	0	A	G	46.0%	84.6 / 11.4	63	noncoding (72/88 nt)	serX	tRNA‑Ser
Reads supporting (aligned to +/- strand):  ref base A (7/27);  new base G (29/0);  total (36/27)
Fisher's exact test for biased strand distribution p-value = 1.12e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAA  >  NC_000913/1097436‑1097728                                                                                                                                                    |                                                                                                                                                     cACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGa                                                                                                                                                 <  1:823923/151‑1 (MQ=31) cACTGCCTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCg                                                                                                                                                  <  1:721724/150‑1 (MQ=14)       gTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGAGTCGAAacccc                                                                                                                                           >  1:506467/1‑151 (MQ=18)          cACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGa                                                                                                                                         <  1:766882/150‑1 (MQ=38)            cGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATAc                                                                                                                                      <  2:172119/151‑1 (MQ=39)                   tCCAACTGGCTGGCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGGCTGG‑GGGGGGGAATTAAAccccccgacgatggc                                                                                                                                 >  1:381227/1‑140 (MQ=11)                    ccAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGGTTGGACCCCCCACACGTTTCCg                                                                                                                               >  1:167618/1‑150 (MQ=11)                    ccAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCg                                                                                                                               <  2:982125/150‑1 (MQ=255)                        cTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAACCCCCCCTACGTTTCCGcata                                                                                                                           >  1:318316/1‑150 (MQ=11)                                 gCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGAACCCCCGCTACGTTGCCGGATACACACttt                                                                                                                   >  2:821129/1‑149 (MQ=255)                                  ccAGTTGTCGACCCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTcc                                                                                                                 <  1:198690/150‑1 (MQ=255)                                    aGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAgg                                                                                                              <  2:1058132/151‑1 (MQ=255)                                       tGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGAACCCCCCATACGTTGCCGGATACACACTTTCCCgg                                                                                                              >  2:364476/1‑148 (MQ=255)                                              ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGCACCCCCCATACGCTGCCGCATAAACACTTTTCCAGCGGGCT‑cc                                                                                                     >  1:244996/1‑149 (MQ=14)                                              ccccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑AGA‑GGGGGGGAGTCGGAACCCCCCTACCTTTCCGCCTACACACTTTCCCGGCGTGcc‑ccct                                                                                                   >  1:546156/1‑149 (MQ=14)                                               cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAACCCCCGCTACGTTGCCGTATACACACTTTTCCGGCGGGGT‑CCtt                                                                                                   >  1:1082468/1‑150 (MQ=17)                                               cccGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGAACCCCCCATACGTTTCCGGATAAACACCTTTCAGGGGGGCT‑cc                                                                                                     >  2:932654/1‑148 (MQ=14)                                                     cGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCa                                                                                             <  1:609744/150‑1 (MQ=255)                                                      ggggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCAc                                                                                            <  2:489576/150‑1 (MQ=255)                                                        ggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TTA‑GGGGGGGGTTTGCAGCCCCGCGCCGGTTCCGTATAAACAATTTCCAAGGGGGGT‑CCTTCTGCCACCCg                                                                                         >  2:408441/1‑151 (MQ=11)                                                        ggCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCg                                                                                         <  1:611611/151‑1 (MQ=255)                                                         gCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCgg                                                                                        <  2:175666/151‑1 (MQ=255)                                                             ctcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGAACCCCCCAAACGTTGCCGGATAAAAACTTTCCAGGCGGGCT‑CCTTCAGCCACCCgg                                                                                        >  1:583029/1‑147 (MQ=17)                                                              tcATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAACCCCCCATACGGTGCCGGATACACACTTTCCACGCGTGCT‑CCTTCAGCCCCCCGGACAcc                                                                                   >  2:949716/1‑151 (MQ=17)                                                               cATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGaca                                                                                     <  2:830535/148‑1 (MQ=255)                                                                aTCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGGTTGCAACCCCCAAACGTTTCCGTATAACCACATTTCAAGCGGGCTCCCTTCACCCAccgggcacac                                                                                   >  1:923277/1‑141 (MQ=2)                                                                      ccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCaaa                                                                           <  2:37641/151‑1 (MQ=255)                                                                      ccGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑ATA‑GGGGGGGGTTCGAACCCCCCATAAGTTGCCCGATAAAAAATTTTCAGGCGTGGT‑CCTTCCGCCACTCGGACACCTCAcc                                                                              >  2:1101231/1‑148 (MQ=17)                                                                           gtgtgCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAAACCCCGATACGTTGCCGTATACACACCTTTCAGGCGTGCT‑CCCTCAGCCCCCCGGGCACCCCCCCAAATTGt                                                                       >  1:703514/1‑150 (MQ=14)                                                                               gCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGtt                                                                   <  1:741506/150‑1 (MQ=255)                                                                                  atatACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGAACCAAACCCCCATACCATTCCGTATAAACACTTTTCCAGCGCGCT‑CCCTCAACCACTTGGACAGCTCTCCAAATTTTTGTTcc                                                                 >  2:260940/1‑149 (MQ=9)                                                                                   tataCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑AGA‑GGGGGGGGGTTGAACCCCCGATACGTTGCCGTATACACAACTTCCAGGCGGGCT‑CCTTCAGCCCCTCGGACACCCCACCAAAATTTCggttct                                                                >  1:195656/1‑144 (MQ=11)                                                                                     taCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGAGTCGAACCCCCCAAAAGTTGCCGTGAAAAAACTTTCCCGGGGGGGG‑CCCTTAGCCCCCCGGGCACCCCACCaaaatttggt                                                                    >  2:336839/1‑136 (MQ=2)                                                                                      aCGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTATTAACCCCGCAATCTTGTCGTATAAAAACTTTACTGGGGTGCT‑CCTTCAGCCCCCCGGACACCTCCCCCACTTTTTGGTCCTGt                                                              >  1:983276/1‑148 (MQ=2)                                                                                       cGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTg                                                          <  2:541311/151‑1 (MQ=255)                                                                                        gAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGGACCCCCCAGAACTTTCCGTATAAACACCTTCCCGGGGTGCT‑CCTTCAGCCACCCGGGAACCCCCCCCAAATGTCGTTCCTGTCtt                                                           >  1:221539/1‑149 (MQ=14)                                                                                         aaaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGAACCCCCGCTACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCCGGCACCCCCCCAAATTTTCGTTCCTGTCtt                                                           >  2:905144/1‑148 (MQ=14)                                                                                         aaaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGGTCGAACCCCCCCTACGTTGCCGCCTACACACTTTCCAGGCGGGCT‑CCCTTCGCCACTCGGGCACCCCACCAAATTGTCGTTCCTGTCtt                                                           >  2:664228/1‑148 (MQ=14)                                                                                         aaaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTg                                                          <  1:1003995/149‑1 (MQ=255)                                                                                          aaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAACCCCCGCTACGTTGCCGGATACACACTTTTCAGGCGGGGT‑CCTTCAGCCCCTCGGGCACCCCACCAAATTGTCGTTTCTGTCTTGCt                                                        >  2:1015486/1‑150 (MQ=14)                                                                                          aaaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTg                                                       <  1:1009477/151‑1 (MQ=255)                                                                                           aaaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTgg                                                      <  2:549437/151‑1 (MQ=255)                                                                                            aaaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑AAA‑GGGGGGGGTTCGAAACCCCCACAACTTTGCGTATAAACACTTTTCACGGGGGGG‑CCCTCCACCACCCGGAGACCTCCCCAAATTTTCGGTCCCGTCtttcctg                                                      >  1:681102/1‑145 (MQ=2)                                                                                             aaaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑AGA‑GGGGGGGGTTTGAAACCCCCCTAGCTTTCCGTATAAAAACTTTCCAGGCGGGCG‑CCTTCACCCACCCGGACCCCCCCCCAAATTGTTGTTCCTTTCTTTCTgg                                                      >  1:1086528/1‑149 (MQ=9)                                                                                              aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGAGGGGGGGGGTTCG‑ACCCCCCCTACGCTGCCGCATTCACACTTTCCAGGCGTGGT‑CCCTCCGTCACTCGGCCACCTCACCCAATTGTCGGTCCTGTCTTgttcg                                                      >  1:180561/1‑144 (MQ=14)                                                                                              aaaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAAc                                                   <  2:497560/151‑1 (MQ=255)                                                                                               aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAACCCCCCATACGTTGCCGGATAAACACCTTCCCGGCGTGCT‑CCCTCCGCCACTCGGGCACCCCACCAAATTTTCGTTTCCGTCttt                                                          >  2:1002035/1‑142 (MQ=14)                                                                                               aaGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACg                                                  <  1:514757/151‑1 (MQ=255)                                                                                                   ccGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGGACCCCCCCTAAGTTGCCGGATAAACACTTTTCCGGCGTGCT‑CCTTCAGCCACCCGGGCAACCCACCCAATTGTCGTTTCTGTCTTTCTGGGACggg                                                >  2:1040841/1‑149 (MQ=14)                                                                                                     gTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGG‑GGGGGGGGTTTGGAACCCCGCTACCTTTCCCTATACACACTTTCCCGGCGGGCC‑CCCTTTCCCACCCGGACACCCCACCAAATTTTCGTTCCTGTCTTGCTGGGACgg                                                 >  1:4368/1‑146 (MQ=9)                                                                                                        cTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTa                                           <  1:757223/149‑1 (MQ=255)                                                                                                             gTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGCACCCCCCGTACGTTTCCGCCTACACACTTTCCACGCGTGCT‑CCTTCTGCCACTCGGGCACCCCTCCAAAATTTCGGTCCTGTCTTTCTGGAAACGGCGCTCAttt                                       >  1:749818/1‑148 (MQ=14)                                                                                                              tACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTTGAACCCCCCCTACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTTAGCCACTCCGACACCCCCCCAAATTTTCGTTCCTGTCTTTCTGGAACGGGCGCTAAattttgg                                   >  2:1015023/1‑147 (MQ=17)                                                                                                               aCGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGTTCGAACCCCCGATACGTTTCCGTATAAACACTTTTCAGGCGGGCT‑CCCTTAGCCACTCGGGCACCTCACCAAAATTTCGTTTCTGTCTTGCTGGGAAGGGCGCTAATTTTggg                                   >  2:232867/1‑150 (MQ=14)                                                                                                                  aGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGGATCCAACCCCCCCTACCTTTCCGTATAAACACCTTCCCGGGGGGGT‑CCTTCAGCCACCCGGGCACCCCACCCAAATTTCGTTCCTGTCTTTCTGGAACGGGCGCCAATTTTGGGaaa                                >  1:184124/1‑150 (MQ=12)                                                                                                                           ttAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCt                       <  1:735847/150‑1 (MQ=255)                                                                                                                            tAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTg                      <  2:562188/150‑1 (MQ=255)                                                                                                                              aaTATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCCGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCCAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAgg                   >  2:60799/1‑151 (MQ=255)                                                                                                                                 atGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTca                <  2:835840/151‑1 (MQ=255)                                                                                                                                 atGGCGG‑TGA‑GGGGGGGATACGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTc                 <  2:1053158/150‑1 (MQ=255)                                                                                                                                      gg‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACaaa           <  1:821112/151‑1 (MQ=255)                                                                                                                                          gA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACtt        <  1:877515/151‑1 (MQ=255)                                                                                                                                                  ggATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGaa  >  1:857882/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                     CACTGCGTTCACGGTCTGTCCAACTGGCTGTCGCCAGTTGTCGAACCCCGGTCGGGGCTTCTCATCCCCCCGGTGTGTGCAATATACGAAAAAAAAGCCCGTACTTTCGTACGAGCTCTTCTTTAAATATGGCGG‑TGA‑GGGGGGGATTCGAACCCCCGATACGTTGCCGTATACACACTTTCCAGGCGTGCT‑CCTTCAGCCACTCGGACACCTCACCAAATTGTCGTTCCTGTCTTGCTGGAACGGGCGCTAATTTAGGGAAATCATGACCTGAGGTCAACAAACTTTTTGAA  >  NC_000913/1097436‑1097728

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