Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,573,500	0	A	G	51.0%	57.0 / 77.6	51	pseudogene (1087/1503 nt)	yjiT	pseudogene
Reads supporting (aligned to +/- strand):  ref base A (4/21);  new base G (26/0);  total (30/21)
Fisher's exact test for biased strand distribution p-value = 1.11e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-09

AGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGG  >  NC_000913/4573355‑4573650                                                                                                                                                   |                                                                                                                                                                 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:346905/151‑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‑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)                                                                                                                                                 ggATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCg    <  1:246087/151‑1 (MQ=255)                                                                                                                                                   atga‑tTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCggg  >  1:150197/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                                 AGATGGTGAGGAAGTCTCTGGCCTTGGGCCAGCCTATGCTTCTCTGGAAAACAGACAGGCAACAGTTCGATTACGTAAAAGCGAAGTGAGATTTGGCAGGCAAAATCCATCGGCAGGTTTGTCGTTAGTGGCTCGT‑GCTGGAGGGATGA‑TTG‑TT‑GGGAGCA‑TTAAACTTGATGACA‑GT‑G‑AAA‑TTG‑CCATTGGT‑GAGGTGCCGTTAACCTTCATCGTTGATGCAGATCAATGGCTGTTACAGGGACAGGCTTCTTGCAGTGTGCGAAGCAGCGATGTTCTGATTGTGCTCCCTCGGG  >  NC_000913/4573355‑4573650

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: ‑