Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	42.9%	101.0 / 11.2	63	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (30/6);  new base C (0/27);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 4.12e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907459‑907750                                                                                                                                                   |                                                                                                                                                     gcgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCa                                                                                                                                               <  1:576221/151‑1 (MQ=255)   ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCAcc                                                                                                                                             >  2:750079/1‑151 (MQ=255)    gTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCg                                                                                                                                            <  2:286822/151‑1 (MQ=255)      atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGcc                                                                                                                                          >  2:862922/1‑151 (MQ=255)      atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGcc                                                                                                                                          <  1:537928/151‑1 (MQ=255)          tcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATa                                                                                                                                       >  1:358573/1‑150 (MQ=255)          tcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATa                                                                                                                                       >  2:618377/1‑150 (MQ=255)           ccGCCAGAAGATCTTTTTTCCGCCCCTGAACCACCGCCCGGGCGGGGGGGGGCGGTACGGAAGTAATCACTTTGAGAAC‑CCACCCCAAAACCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CACCCAGCGCCGCTTCCCCCTCACCGCCATaa                                                                                                                                      <  2:863314/150‑1 (MQ=255)                     aTCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCaaa                                                                                                                            >  2:380805/1‑150 (MQ=255)                     aTCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAAc                                                                                                                           >  1:391874/1‑151 (MQ=255)                         tgtTTCACCCCCTGAATCACCGCCCTGGGGGGGGGGGGCGGTACGGAAGTAATCACATTAAAAAC‑CCCCCCCAAAACCCGGGGGGTGGGGG‑GGGGGAGGGGGAC‑CACCCAGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAAc                                                                                                                           <  1:79956/147‑1 (MQ=255)                           tttcccccccTGAATAACCGCCCTGGCGGTGGTGGGCGGTAGGAAAGTAACCAAATTGAGAAC‑CCACCCCAAAATCCAGGGGGGGGGGG‑GGGGGAGGGGGAA‑CAACCGGCGGCGTTTCCCGCGCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                      <  1:96602/144‑1 (MQ=255)                            ttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGCAGTAATCACATTGAGAAA‑CCCCACCGACATCCAGCGGGGGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                      >  2:487976/1‑149 (MQ=255)                              cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                   >  2:402823/1‑150 (MQ=255)                                   cgtaaaccaccccccGGGGGGGGGGGGGCGGTCGGGAAGTAACCAATTTGA‑AAA‑CCCCCCCCAAAACCCGGGGGGGGGGAGGGGGGTTGGGAAC‑CAGCCGGGGGCGCTTCCCCCGAACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                             <  2:78728/139‑1 (MQ=255)                                   ccTGAACCACCGCCCTGGCGGTGGGGGGGGGTAGGAAAGAAATCACATTAAAAAC‑CCCCCCCGAAACCCGGGGGGTGGGGG‑GGGGGTGGGGAAT‑CACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                             <  2:946063/151‑1 (MQ=255)                                    cTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                            >  1:1027813/1‑151 (MQ=255)                                      gAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTg                                                                                                            >  2:1032615/1‑149 (MQ=255)                                          cAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAcc                                                                                                      >  1:867628/1‑151 (MQ=255)                                           atggccggggggggggtggGAGGGACGAAAGAAACCGTATCGAAAAA‑CCCCACCAAAATCAAG‑GGGGGGGGG‑GGGGGGGGGGGATACACCCAGCGCGGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCAc                                                                                                       <  1:489464/135‑1 (MQ=255)                                             cGCCCGGGGGGTGGGGGGGGGAACGAAAGAAACCAC‑TTGAAAAACCCCCCCCAAAACCCAGGGGGGGGGGG‑GGGGGAGGGGGAC‑AACCCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctc                                                                                                    <  1:458188/150‑1 (MQ=255)                                               cccGGGGGGTGGGGGGCGGTACGGAAGAAATCAGATTGAAAAA‑CCACCCCAAAACCCGGGGGGTGGGGG‑GGGGGGTGGGGAA‑CAACCGGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCa                                                                                                        <  1:1032657/144‑1 (MQ=255)                                                ccTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGt                                                                                                >  2:356677/1‑151 (MQ=255)                                                  tGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                              >  1:921809/1‑151 (MQ=255)                                                  tGGCGGAGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGG‑GGGGGAGGGGGAT‑CAACCAGGGGCGCTTACCGCACACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGATCACCTCTGGTc                                                                                               >  2:469677/1‑150 (MQ=255)                                                   ggggtgggtggGGGGTGGGAAAGTAACCACATTGAAAAA‑CCACCCCAAAACCCGGGGGGGGGGGG‑GGGGGTTGGGGAT‑AAACCAGCGGCGCTCCCCCCCAACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTcc                                                                                              <  2:943624/144‑1 (MQ=255)                                                     cGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                           >  2:666081/1‑151 (MQ=255)                                                      ggtggtGGGCGGTAGGAAAGTAATCACATTGAAAAA‑CCACCCCAAACTCCAGCGGGTGGGGG‑GGGGGAGGGGTAT‑CACCCAGGGGCGCTTACCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                             <  1:629590/148‑1 (MQ=255)                                                         ggtggGCGTTGCGGAAGTATTCCCATTGAAAAC‑CCCCCCCAAAACCCGGCGGGTGGGGG‑GGGGAATGGGGAT‑CACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCt                                                                                       <  2:47585/151‑1 (MQ=255)                                                            gggCGGTACGAAAGTAATCAAATGAAAAAC‑CCCCCCCGAAACCCGGGGGGGGGGGG‑GGGGGAGGGGAAC‑CAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTc                                                                                     <  1:609140/150‑1 (MQ=255)                                                              gCGGTACGGAAGTAATCACATTGAAAAC‑CCACCCCGAAATCCGGGGGGGGGGGA‑GGGGGAGGGGGAT‑CAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGaa                                                                                  <  1:232589/151‑1 (MQ=255)                                                                  tACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATg                                                                              >  1:756927/1‑151 (MQ=255)                                                                      gAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgt                                                                          >  2:399467/1‑151 (MQ=255)                                                                       aaGTAATCACATTGAGAAA‑CCACACCGAAACCCAGGGGGTGGGGA‑GGGGGAGGGGGAT‑CAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCgg                                                                           <  1:759824/149‑1 (MQ=255)                                                                         gTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGtt                                                                       >  2:41582/1‑151 (MQ=255)                                                                          tAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                      >  2:374055/1‑151 (MQ=255)                                                                               acaTTGAAAAC‑CCACCCCAAAACCCAGGGGGGGGGGG‑GGGGGTGGGGGAT‑CACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGAGGTGTTCGCt                                                                   <  1:1024404/149‑1 (MQ=255)                                                                               acaTGGAAAAC‑CCCCCCCAAAACCCGGGGGGGGGGGA‑GGGGGATGGGGAT‑CACCCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTg                                                                  <  1:840390/150‑1 (MQ=255)                                                                                  ttGAGAAC‑CCACCCCCAAACCAAGCGGGGGGGGA‑GGGGGGTGGGGAT‑CACCCAGCGCCGCTTCCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                      <  1:35710/143‑1 (MQ=255)                                                                                   tGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                             >  2:177309/1‑151 (MQ=255)                                                                                      aaaacccccccaaaaaCCCAGGGGGGGGGGG‑GGGGGAGGGGGAT‑CAACCAGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                           <  1:1026671/138‑1 (MQ=255)                                                                                        aacccccccagAAACCCGGCGGGTGGGGG‑GGGGGATGGGGAT‑CACCCGGCGGCGCTTCCCGCTCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgatg                                                          <  2:478505/140‑1 (MQ=255)                                                                                        aacccACACCAAAACCCAGGGGGGGGGGG‑GGGGGAGGGGGAT‑CACCCGGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                           <  2:266373/146‑1 (MQ=255)                                                                                                ccGAAACCCAGGGGGTGGGGG‑GGGGGAGGGGAAT‑CACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                                                    <  1:329723/148‑1 (MQ=255)                                                                                                   aaaTCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                              >  2:788617/1‑151 (MQ=255)                                                                                                    aacccAGCGGGGGGGGG‑GGGGGAGGGGAAC‑CAACCGGCGGCGCTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                             <  2:68132/148‑1 (MQ=255)                                                                                                     acccAGGGGGGGGGGG‑GGGGGAGGGGGAT‑CACCCCGCGGCCCTCACCCCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCaa                                            <  2:58126/149‑1 (MQ=255)                                                                                                     aTCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                             >  2:929761/1‑150 (MQ=255)                                                                                                      tCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                           >  2:897118/1‑151 (MQ=255)                                                                                                      tCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAc                                           >  2:699121/1‑151 (MQ=255)                                                                                                          gCGGGGGGGGA‑GGGGGAGGGGGAT‑CAACCAGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCa                                         <  1:463979/149‑1 (MQ=255)                                                                                                           cGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAAcc                                          >  2:522671/1‑147 (MQ=255)                                                                                                               tGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTg                                  >  1:55918/1‑151 (MQ=255)                                                                                                                 gggg‑ggggGAGGGGAAT‑CAACCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgca                                <  2:342111/151‑1 (MQ=255)                                                                                                                          gatggggat‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTAGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                        >  2:479817/1‑151 (MQ=255)                                                                                                                            gggggat‑CACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTc                      <  1:927821/150‑1 (MQ=255)                                                                                                                             ggggac‑CACACAGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGtttt                       <  1:598438/149‑1 (MQ=255)                                                                                                                              gggat‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                     >  1:354693/1‑150 (MQ=255)                                                                                                                                     aaCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa               <  1:430028/150‑1 (MQ=255)                                                                                                                                       ccAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCaa              <  1:148985/149‑1 (MQ=255)                                                                                                                                       ccAGCGGCCCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATg            <  2:429222/151‑1 (MQ=255)                                                                                                                                         aGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg          >  2:113974/1‑151 (MQ=255)                                                                                                                                                 ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc  >  2:58622/1‑151 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAA‑CCACACCGAAATCCAGCGGGTGGGGA‑GGGGGATGGGGAT‑CAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGC  >  NC_000913/907459‑907750

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