Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	802,111	0	G	T	39.1%	74.9 / 24.9	64	A75A (GCG→GCT)	ybhI	putative DASS family tricarboxylate or dicarboxylate transporter
Reads supporting (aligned to +/- strand):  ref base G (23/16);  new base T (0/25);  total (23/41)
Fisher's exact test for biased strand distribution p-value = 3.60e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCAGGATTAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCA  >  NC_000913/801962‑802256                                                                                                                                                       |                                                                                                                                                  gCAGGATTAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGt                                                                                                                                                  <  2:959224/151‑1 (MQ=255)  cAGGATTAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGtt                                                                                                                                                 <  1:645482/151‑1 (MQ=255)    ggATTAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGttt                                                                                                                                                >  1:778550/1‑150 (MQ=255)      aTTAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTaaa                                                                                                                                             <  1:229526/151‑1 (MQ=255)                cTGGCGTGGGGGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCACGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCAccgcc                                                                                                                                    >  2:758266/1‑150 (MQ=255)                   gCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTAt                                                                                                                                >  1:135790/1‑151 (MQ=255)                      tGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTa                                                                                                                              >  1:881928/1‑150 (MQ=255)                          tGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCgg                                                                                                                         >  2:307805/1‑151 (MQ=255)                             tttttGTTATTTCCCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTCCCCGGAC‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAAGGGGGGTGGTCGGTAATTTACCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGtt                                                                                                                       <  2:54790/150‑1 (MQ=255)                              tttgGGTATTTCCCTGGCGGCCTTTGTGGGGCGGGTTATCAA‑CCTTTCCCCGAACCCGGTCTTATTGTTAATGGCCTTTGCTGCCTCAATGGTGGGGGTCGGTAATTTACCCGCCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGtt                                                                                                                       <  1:234232/149‑1 (MQ=255)                                 tggtttttccgggggggcAATTGTGGGCGTGTTAACCAAGCTTTTCCCGGAA‑CTTTTCGTCTTTTTAATTCCCGTTGCTCCCTCAATGGTGTTGGTGGTTAACTTAGCCGACGGGGCGTTTAAACCCACCGCCGTATTAAGCGGTTActc                                                                                                                   <  2:917679/135‑1 (MQ=255)                                 ggttttttCCCGGGGGGCCTTTGGGGGGCTTGTTATAAGGCCTTCCCCGGAC‑CTTTTCGTACTTTTAATTGCCGTTCCTCCCTAATGGGTGGTGGTGGGTACTTTTCCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTActc                                                                                                                   <  2:1171811/145‑1 (MQ=255)                                    tATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTc                                                                                                                >  1:115864/1‑150 (MQ=255)                                       ttcccgggcggcCATTGTGGGGCTGTTTAACAAGCTTTCCCCGGAA‑CCTGTCGTACTGTTAATGGCCGTTGCTGCCTAAATGGGGGGGGTGGGTAACTTACCCGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                             <  1:155165/144‑1 (MQ=255)                                       ttACCTGGCGGCGTTTGGGGGGTTGTTTATCAAGCTTTTCCCGGAC‑CCTTTGGTATTGTAAATTGCTGTTGCTGCCTAAAGGGTGGTGGTCGTTACCTTACCGGACGTTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAgg                                                                                                             <  1:582837/150‑1 (MQ=255)                                          ccgggcggcCTTTGTGGGGCTGGTTATCAAGCTTTCCCCGGAC‑CCTGTCGTACTGTTAATTGCCGTTGCTCCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGGGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTa                                                                                                           <  1:1116233/146‑1 (MQ=255)                                          ccTGGCGCCCTTTTGGGGGCTGTTTAACAAGCTTTTCCGGAAC‑CCTGTCGTCTTGTTAATTGCGTTTGCTGCCCCAATGGTGGTGTTCGGTAATTTACCCGCCGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGt                                                                                                            <  1:574692/148‑1 (MQ=255)                                            tGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                       >  2:882617/1‑151 (MQ=255)                                             ggcggcCTTTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTTTCGTACTTTTAATTGCCGTTGCTGCCTAAAGGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTaccac                                                                                                       <  1:1048852/150‑1 (MQ=255)                                                       gtgGGGCTGGTTATCAAGCTTTCCCCGGAC‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTTACCCGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgt                                                                                              <  1:256230/149‑1 (MQ=255)                                                       gtgGGGCTGGTTATAAACCTTTTCCCGGAC‑CCTGTCTTATTTTTATTTGCCTTTGCTGCCTCAATGGTGGTGGTCGGTACCTTACCCGACGGTGCTTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGgtgt                                                                                            <  1:996349/151‑1 (MQ=255)                                                        tgGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTCCCACCTGGCTGGTGtt                                                                                           >  2:242876/1‑151 (MQ=255)                                                         ggggTTGGTTATCAAGCTTTCCCCGGAA‑CTTGTCTTATTGTTAATTGCCGTTGCTGCCTAAATGGTGGTGGTCGGTAATTTATCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGtt                                                                                           <  1:135459/150‑1 (MQ=255)                                                              tGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGACGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTtctc                                                                                        >  2:795645/1‑148 (MQ=255)                                                                   aTCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                 >  2:399578/1‑150 (MQ=255)                                                                      atgcttttCCGAGAT‑CCTGTCGTATTGTTAATTACCTTCGCTGCCACAATGGTGGTGGTGGAAATCTTAGCCGTGGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTAcc                                                                              <  2:268800/145‑1 (MQ=255)                                                                      aatccTTTCCCGAAA‑CCGGCCGTTTTTTAAAGCGCCGTTGCGGCCTCATTGGGGTGGGTCGGTAACTTATCCGACGGTGCTTTAAAACCCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGttt                                                                                 <  1:117902/144‑1 (MQ=255)                                                                         ccTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTa                                                                           >  2:133098/1‑150 (MQ=255)                                                                           ttccccGGAA‑CCTTTCGTTCTTTAAATTGCCTTTGCTCCCTCAATGGTGGTGGTGGGAAACTTACCCGCCGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                            <  1:576432/144‑1 (MQ=255)                                                                             tCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgcgc                                                                      >  1:385241/1‑151 (MQ=255)                                                                              cccGGAC‑CTTGTCGTACTTTTAATTCCCTTTGCTCCCTAAAGGGTGGTGGTCGGTAACTTATCCGACGGGGTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCtt                                                                            <  2:219045/144‑1 (MQ=255)                                                                              cccGGAC‑CCTGTGGTACTTTTATTTGCCTTTGCCGCCTCAAGGGTGGTGGTGGGTATTTTATCGGCCGGTGTTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAgc                                                                        <  2:1058433/148‑1 (MQ=255)                                                                               ccGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGcgccg                                                                    >  2:584112/1‑151 (MQ=255)                                                                                    acccTTTCGTACTGTTATTGGCCTTTGCTGCCACAATGGGGGTGGTCGGTATCTTATCCGAGGGTCTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCa                                                                  <  1:539995/147‑1 (MQ=255)                                                                                       cTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTtgt                                                             >  1:849749/1‑151 (MQ=255)                                                                                            gTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacc                                                         >  1:501278/1‑150 (MQ=255)                                                                                             tACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                       >  1:726645/1‑151 (MQ=255)                                                                                              actttaatttGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAATTAACCCGCCGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGaccac                                                       <  1:509861/142‑1 (MQ=255)                                                                                                tGTAAATGGCCGTTGCGCCCTAAAGGGGGGGGTGCGGTAACTTTCCCCACGTTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGacca                                                        <  2:236949/147‑1 (MQ=255)                                                                                                         ccGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTa                                            >  2:968003/1‑150 (MQ=255)                                                                                                         ccGTTGCTGCCTCAAGGGTGGGGGTCGTTAACTTACCCGACGGTGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGt                                             <  1:290163/149‑1 (MQ=255)                                                                                                            ttGCTCCCTCAAGGGGGTTGGTCGGTAATTTATCCGGGGGTGCTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAAc                                         <  2:583930/150‑1 (MQ=255)                                                                                                             tgctgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                       >  1:241880/1‑151 (MQ=255)                                                                                                                tgcCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                    >  1:474396/1‑151 (MQ=255)                                                                                                                 gcCTAAAGGGGGGTGTTGGTTAACTTATCCGACGGTGCGTTAAAAACCTCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGTCCACCGGTTTAGGTAAAc                                         <  1:627812/145‑1 (MQ=255)                                                                                                                  ccTAAATGGTGGTGGTGGGTACTTTCTCCGAGGGGGCTTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTAtt                                    <  2:387393/149‑1 (MQ=255)                                                                                                                   ctaattggtgggggTGGTTAACTTACCCGAGGGTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGt                                       <  1:897521/140‑1 (MQ=255)                                                                                                                      aaTGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTAt                              >  2:627420/1‑151 (MQ=255)                                                                                                                      aaTGGTGGTGGTCGGTAACTTACCCGACGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTa                               <  1:201588/150‑1 (MQ=255)                                                                                                                         ggtggtggtCGGTATCTTACCCGACGGTGCGTTAAAACCCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATct                            <  1:882677/150‑1 (MQ=255)                                                                                                                         gggggtggtCGTTATCTTACCCGACGGTGGTTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctg                           <  2:1154830/148‑1 (MQ=255)                                                                                                                               ggttggtttcTTACCCGAGGGTGGTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATctgc                          <  2:1105178/137‑1 (MQ=255)                                                                                                                                gtgggtattttACCCGGCGGTGCTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTg                    <  1:584140/142‑1 (MQ=255)                                                                                                                                  cGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                  >  2:862028/1‑151 (MQ=255)                                                                                                                                      aaCTTTCCCGTCGGTGCTTAAAAACCCCCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGAt                      <  2:1177818/143‑1 (MQ=255)                                                                                                                                         ttACCCGACGGTGTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGt                  <  1:339756/144‑1 (MQ=255)                                                                                                                                           acccGACGGTGCGTTAAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg          <  1:806405/148‑1 (MQ=255)                                                                                                                                           aTCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCgg          >  1:114505/1‑150 (MQ=255)                                                                                                                                             ccgacgtgtcGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa        <  2:137733/141‑1 (MQ=255)                                                                                                                                             ccGAGGGTGTTTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTaa       <  2:1169137/151‑1 (MQ=255)                                                                                                                                             ccGACGGTGCGTTTAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGt         <  2:455227/149‑1 (MQ=255)                                                                                                                                             ccGACGGTGCGTTTAAAACCGCCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTa        <  2:778510/150‑1 (MQ=255)                                                                                                                                                aCGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAacac    >  2:1139694/1‑151 (MQ=255)                                                                                                                                                  ggTGCGTTAAAACCCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAAcacca  <  2:447237/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                  GCAGGATTAAGCGAACTGGCGTGGGTGCTTTTTGGTATTTACCTGGCGGCCATTGTGGGGCTGGTTATCAAGCCTTTCCCGGAA‑CCTGTCGTACTGTTAATTGCCGTTGCTGCCTCAATGGTGGTGGTCGGTAACTTATCCGACGGTGCGTTTAAAACCACCGCCGTATTAAGCGGTTACTCTTCAGGTACCACCTGGCTGGTGTTCTCGGCGTTTACCTTAAGCGCCGCATTTGTGACCACCGGTTTAGGTAAACGTATTGCCTATCTGCTGATTGGTAAAATCGGTAACACCA  >  NC_000913/801962‑802256

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