Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,496	0	A	G	51.0%	32.0 / 75.8	50	pseudogene (1083/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand):  ref base A (4/20);  new base G (25/0);  total (30/20)
Fisher's exact test for biased strand distribution p-value = 3.76e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.16e-09

CATTTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCC  >  NC_000913/4573346‑4573644                                                                                                                                                        |                                                                                                                                                               cATTTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGa                                                                                                                                                                <  1:281450/151‑1 (MQ=255)          aGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GGTGGGGGGGGGA‑tt                                                                                                                                                      >  1:364744/1‑151 (MQ=255)           gATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑ttg                                                                                                                                                     <  2:356643/151‑1 (MQ=255)                     aaGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGGGGGGTGG‑TTT‑TT‑tgggggat                                                                                                                                         >  1:151110/1‑146 (MQ=255)                         ctctGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GCTGGGGGGGTGA‑TTG‑TT‑GGGGGCA‑tt                                                                                                                                       >  1:175725/1‑148 (MQ=255)                              gCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGGGGGGTGG‑TTG‑TT‑GGGGGGA‑TTAAAAtt                                                                                                                                 >  1:165068/1‑149 (MQ=255)                                   gggCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACa                                                                                                                          <  2:493252/151‑1 (MQ=255)                                        aGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGGGGGGTGA‑TTG‑TT‑GGGGGGA‑TTAAACTTGATGAAA‑ggg                                                                                                                      >  1:691050/1‑147 (MQ=255)                                        aGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GGTGGGGGGGTGA‑TTT‑TT‑GGGGGCA‑TTAAAAtt                                                                                                                                 >  1:239546/1‑139 (MQ=255)                                               gCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCAT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑cc                                                                                                          <  2:599754/150‑1 (MQ=255)                                                     ctGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGGGGGGGGG‑TTG‑TT‑GGGGGGA‑TTAAAATTGATGACA‑GT‑G‑AAA‑TTG‑GCATTGGT‑g                                                                                                  >  1:145273/1‑151 (MQ=255)                                                       ggAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGt                                                                                                    <  2:343677/148‑1 (MQ=255)                                                         aaaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGGGGGGAGA‑TTT‑TT‑GGGGGGATTTAAATTTGTTG‑CA‑GTGG‑AAA‑TTG‑CATTTggggggg                                                                                                >  2:74719/1‑146 (MQ=255)                                                          aaaCAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGACTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAgg                                                                                               <  2:37042/149‑1 (MQ=255)                                                              agacagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCG‑GGTGGGGGGGGGG‑TTT‑TT‑GGGGGGA‑TTAAAATTGTTGGCAGGG‑G‑AAA‑TGG‑CCTTTgggggggg                                                                                               >  2:6226/1‑141 (MQ=255)                                                                acagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGGGGGGGGA‑TTG‑TT‑GGGGGGA‑TTTAAATTTATGACA‑GTGG‑AAA‑TTT‑CCTTTGG‑‑GAGGGGCCGtt                                                                                        >  1:553682/1‑150 (MQ=255)                                                                acagGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCCG‑GGTGGGGGGGGGA‑TTT‑TT‑GGGGGGA‑TTAAAATTTTTGACA‑AG‑G‑AAA‑ATG‑GCATTTGG‑GAgggggcggt                                                                                        >  2:691385/1‑143 (MQ=255)                                                                  agGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGAGGGGTGA‑TTG‑TT‑GGGGGGA‑TTAAAATTGATGACA‑GG‑G‑AAA‑ATG‑CCATTGGT‑gggggggcggt                                                                                        >  1:247499/1‑141 (MQ=255)                                                                  agGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GGTGGGGGGGAGATTTT‑TT‑GGGGGAT‑TTAAACTTGTTGA‑A‑GG‑G‑AAATTTG‑CCTTTgggggggggcggt                                                                                          >  2:743209/1‑137 (MQ=255)                                                                      aaCAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GCTGGGGGGGTGG‑TTT‑TT‑TGGGGCA‑TTAAAATTTTTGACA‑GT‑G‑AAA‑TTT‑CCCTTGGG‑GGGGGGCCGTTTACCtt                                                                                  >  1:418995/1‑150 (MQ=255)                                                                             cGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGCGGGGGGGTGG‑TTG‑TT‑TGGGGCA‑TTAAACTTTCTGGCA‑GT‑GAAAA‑TGG‑CCTTTGGG‑GGGGGGCGTTTAACTTTctcctt                                                                            >  1:392994/1‑145 (MQ=255)                                                                             cGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTg                                                                          <  1:134986/151‑1 (MQ=255)                                                                              gATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GGTTGGGGGGAGA‑TTG‑TT‑GGGGGGA‑TTATAATTGTTGGAA‑CT‑G‑AAA‑ATG‑TCCTTGGG‑GGGGGGGCGTTAAACTTccttctt                                                                           >  2:552689/1‑143 (MQ=255)                                                                              gATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GCTGGGGGGGTGA‑TTG‑TT‑GGGGGGA‑TTAAACTTGATGACA‑GG‑G‑AAA‑TTG‑CCCTTGGG‑GAGGGGGCGTTTAACTTCAttgttg                                                                          >  1:283532/1‑150 (MQ=255)                                                                                 tACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATg                                                                       <  1:113736/150‑1 (MQ=255)                                                                                  aCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCa                                                                     <  1:229372/151‑1 (MQ=255)                                                                                  aCGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GGTGGGGGGGTGA‑TTG‑TT‑GGGGGGA‑TTAAAATTTGTTACA‑GT‑G‑AAA‑ATG‑TCCTTTGG‑GGGGGGCCGTTTACCCTCCTCGtttgt                                                                        >  1:239988/1‑145 (MQ=255)                                                                                   cGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGGGGGAGAA‑ATA‑TT‑GGGGGGA‑TTAAAAATTATGAAC‑AT‑G‑AAA‑TTT‑GCATTGGG‑GGGGGGGCGGTAACCTTTATAGTTGAt                                                                        >  2:200560/1‑147 (MQ=255)                                                                                        aaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGTGGCTGG‑GGGGGGA‑TTT‑TTGGGGGGCT‑TAAAACTTTTGGAAG‑GG‑G‑AAA‑TTG‑‑CATTGGT‑GGGGTGGCGTTTAACTTCCTCGTTTATGCAggt                                                                  >  2:415518/1‑146 (MQ=255)                                                                                        aaGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GCTGGGGGGGTGG‑TTG‑TT‑TGGGGGA‑TTTAAATTTGTGACA‑GT‑G‑AAA‑ATT‑CCCTTGGTGGGGGGGCCGTT‑ACCTTCATCGTTTATGCa                                                                     >  1:253438/1‑145 (MQ=255)                                                                                          gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGGGGGGTGG‑TTG‑TT‑TGGAGCA‑TTAAAACTTATGACA‑GT‑G‑GAA‑ATG‑CCCTTGGG‑GGGGGGCCGTTAACCCTTATCGTTGGTGGAGAACaa                                                               >  1:583535/1‑149 (MQ=255)                                                                                          gCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATg                                                             <  2:24388/151‑1 (MQ=255)                                                                                            gAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGGGGGGTGG‑TTG‑TT‑TGGAGCA‑TTAAAATTTATGACA‑GT‑G‑GAA‑TTT‑CCATTTGT‑GGGGGGCCGTTTAACTTCATCGTTGATGCAGATCAATgg                                                            >  1:292303/1‑150 (MQ=255)                                                                                                 gagaTTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGtt                                                       <  2:648181/150‑1 (MQ=255)                                                                                                      ttGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GGTGGGGGGGTGA‑TTGTTT‑GGGGGCT‑TTAAA‑TTTATGAAAGGG‑G‑AAA‑TTG‑CCTTTGGG‑GGGGGGCGGTTAA‑CTTTATCTTTTATGGAGATCAATTGCTGTTTCAggg                                                 >  1:608648/1‑151 (MQ=255)                                                                                                        ggcaggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGAc                                               <  1:66107/151‑1 (MQ=255)                                                                                                           aggcaAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GGTGGGGGGGTGA‑TTT‑TT‑GGGGGGA‑TTTAAATTTATGAAA‑GG‑G‑GAA‑ATG‑CCCTTGGG‑GAGGGGCCGTTAAACTTCATCGTTGATGGAGATCAATGGGTGTTTCAGGGGCAgg                                            >  1:651173/1‑151 (MQ=255)                                                                                                              caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGctt                                         <  1:6226/151‑1 (MQ=255)                                                                                                              caAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGct                                          <  1:46211/150‑1 (MQ=255)                                                                                                               aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCGGGGGGGGGGG‑TTT‑TT‑GGGGGCA‑TTTAAATTTTTGACA‑GT‑G‑AAA‑ATGTCCTTTGTG‑GGGGGGCCGTT‑ACCTTTCTTGTTTATGCCGATAAAAAGCTGTTTTAGGGGCAGGctt                                         >  2:466306/1‑150 (MQ=255)                                                                                                               aaaaTCCATCGGCAGGTTTGTCGTTAGTGGCTCGG‑GCTGGGGGGGTGA‑TTG‑TT‑GGGGGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑ATT‑CCATTGGG‑GAGGGGCCGGTAAACTTCATCGTTGTTGCAGAGAAAAGGCTGTTTAAGGGACAGGctt                                         >  1:228554/1‑150 (MQ=255)                                                                                                                            aGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGa                           <  1:440609/151‑1 (MQ=255)                                                                                                                                 tGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagca                       <  1:693033/150‑1 (MQ=255)                                                                                                                                  gTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagcagc                     <  1:291005/151‑1 (MQ=255)                                                                                                                                  gTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagcagc                     <  1:346905/151‑1 (MQ=255)                                                                                                                                  gTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAagcag                      <  2:591745/150‑1 (MQ=255)                                                                                                                                         gTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGGGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGGGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGtt               >  2:730496/1‑150 (MQ=255)                                                                                                                                             cTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGAt          >  1:264030/1‑151 (MQ=255)                                                                                                                                               cGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtg        <  2:36630/151‑1 (MQ=255)                                                                                                                                                gT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATtgt       >  1:498111/1‑151 (MQ=255)                                                                                                                                                      ggAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTcc  <  1:658277/151‑1 (MQ=255)                                                                                                                                                        |                                                                                                                                                               CATTTGTGAAGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCC  >  NC_000913/4573346‑4573644

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