Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	907,602	0	A	C	41.0%	162.0 / 10.8	83	G227G (GGT→GGG)	ybjT	putative NAD‑dependent oxidoreductase
Reads supporting (aligned to +/- strand):  ref base A (40/9);  new base C (0/34);  total (40/43)
Fisher's exact test for biased strand distribution p-value = 2.58e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCG  >  NC_000913/907454‑907751                                                                                                                                                        |                                                                                                                                                      gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTAcc                                                                                                                                                     <  4:349341/151‑1 (MQ=255) gTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTAcc                                                                                                                                                     <  3:671047/151‑1 (MQ=255)   aGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTAcca                                                                                                                                                    <  4:81132/150‑2 (MQ=255)       cgGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCAc                                                                                                                                               >  4:46479/1‑151 (MQ=255)        ggTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCAcc                                                                                                                                              >  3:698003/1‑151 (MQ=255)           atcatcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGcc                                                                                                                                           >  1:43511/1‑151 (MQ=255)             catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCa                                                                                                                                          <  3:460146/150‑1 (MQ=255)             catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCAt                                                                                                                                         >  3:759990/1‑151 (MQ=255)             catcCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCAt                                                                                                                                         >  3:1245247/1‑151 (MQ=255)                   ccAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGGGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATg                                                                                                                                   >  3:237266/1‑151 (MQ=255)                     agcagTTCGTTTTTCACCCCCTAAATAACCGCCCTGGGGGGGGGGGGCGGTACGGAAGTAATCACATTGAAAACCCCCCCCAAAA‑CCCG‑GCGGGGGGGGAGGGGGA‑GGGGGATCAACCAGCGGCGCTTCCCGCTCACCGCCATAAAATGtt                                                                                                                                 <  4:144977/151‑1 (MQ=255)                        agATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCa                                                                                                                               >  1:152583/1‑150 (MQ=255)                         gATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCaa                                                                                                                              >  2:536135/1‑150 (MQ=255)                           tCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACt                                                                                                                           >  2:245251/1‑151 (MQ=255)                               gtTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                       >  2:556066/1‑151 (MQ=255)                                tttCCGCCCCTAAATCACCGCCCTGGGGGGGGGGGGCGGTAGGGAAGTAATCACAGTGAAAACCCACCCCGAAA‑TCCA‑GCGGGTGGGGAGGGGGT‑GGGGGATCAACCGGCGCCGCTTCCCCCACCCCGCCATAAAATGTTCAAACTGTTg                                                                                                                       <  4:78384/150‑1 (MQ=255)                                tttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                                      >  2:142450/1‑151 (MQ=255)                                tttCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGGGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTg                                                                                                                       >  3:559058/1‑150 (MQ=255)                                 ttCAGCCCCTGAATCACCGCCCGGGGGGTGGGGGGCGGTACGGAAGTAATCACATTGAAAACCCACCCCAAAA‑TCCG‑GGGGGGGGGGGGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCCCTCACCGCCATAAAATGTTCAAACTGTTGc                                                                                                                      <  2:516134/150‑1 (MQ=255)                                   cAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTg                                                                                                                    >  1:400225/1‑150 (MQ=255)                                        ccgaaaccaccGCCGTGGCGGTGGTGGGCGGGACGAAGGAAATAACATGAAGAACCCCCACGGAAC‑CCCG‑GGGGGGGGGGGGGGGGT‑GGGGAATCAGCCGGCGGCGCTCCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCt                                                                                                                     <  4:283992/134‑1 (MQ=255)                                            aaTCACCGCCCTGGGGGTGGTGGGGGGTAGGAAACTAAACCCAAGGAAAACCCACCCCAAAA‑CCGA‑GGGGGTGGGGAGGGGGA‑GGGGAATCACCCGGCGGCTCTTCCCTCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACt                                                                                                              <  1:731558/147‑1 (MQ=255)                                            aaTCAACGCCCTGGCGGTGGCGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGc                                                                                                          >  4:168960/1‑151 (MQ=255)                                              acaacgcagagggggtggaGGGTGGTAGAGAATTAATAACATCAAACACCCAC‑CCGAAATTCCA‑GCGGGT‑GGGAGGTGGATTGGGGACAGGACAGCGGCATTGCCCCCTCACCGCCATAAAATGTTCAAACAGTTGCTGATa                                                                                                                 <  4:329014/129‑1 (MQ=255)                                                aaCGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACct                                                                                                      >  3:537219/1‑151 (MQ=255)                                                  cGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                    >  1:133179/1‑151 (MQ=255)                                                    ccggggggggggtggGGGGTACGGAATTAACCTCATAAAAAACCCCCCCCGAAA‑TCCT‑GGGGGGGGGGAGGGGGA‑GGGGGTGCAACAGGCGGGCCTAACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGa                                                                                                            <  4:376520/131‑1 (MQ=255)                                                    cccGGGCGTGGGGGGGGGGTACGAAAGAATCCACATAGAAAACCCCCACCAAAC‑CCCA‑GGGGGGGGGGAGGGGTT‑GGGAAACCACCCGGCGCCCCTTCCCTCCCTCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                    <  1:744498/149‑1 (MQ=255)                                                        ggCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcg                                                                                              >  3:1020979/1‑151 (MQ=255)                                                         gCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcgg                                                                                             >  2:506004/1‑151 (MQ=255)                                                           ggtggtGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                           >  1:923378/1‑151 (MQ=255)                                                           ggtggtGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                           >  1:1017973/1‑151 (MQ=255)                                                           ggggggggggggTAGGAAATTAACCACATGGAAAACCCCCCCCAAAA‑CCCG‑GGGGGGGGGGGGGGGGA‑GGGGGACCAACCGGCGGCTCTCCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                           <  1:924066/145‑1 (MQ=255)                                                           gggggggggCGGTACGGAAGCAAACACA‑GGAAGAACCCCCCCGAAA‑ACCAGGCGGGGGGGGAGGGGGA‑GGGGGAAAAACCCGCGGCGTTCACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcg                                                                                           <  3:862718/145‑1 (MQ=255)                                                            gtggtgGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                          >  1:323801/1‑151 (MQ=255)                                                            gggggggggggTAGGGAATTAATCAAATGGAAAACCCCCACCAAAA‑CCCG‑GGGGGGGGGGGGGGGGG‑GGGGGATCACCCCGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggcgg                                                                                          <  1:557523/146‑1 (MQ=255)                                                              gggggggggTAGGAAATTAATCACATTGAAAACCCCCCCCAAAA‑CCCA‑GGGGGGGGGGGGGGGGG‑GGGGGACCACCCAGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCcggc                                                                                            <  3:410000/144‑1 (MQ=255)                                                                 ggggggTACGAAAGTAATCACATGGAAAACCCCCCCCGAAA‑CCCG‑GCGGGGGGGGGGGGGGA‑TGGGGATCACCCGGCGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                                                     <  2:94695/151‑1 (MQ=255)                                                                    cGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCg                                                                                     >  3:383814/1‑148 (MQ=255)                                                                    cGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGa                                                                                    >  2:564204/1‑149 (MQ=255)                                                                     ggTAGGAAAGAAATCAAATTAAAAACCCACCCAAAAA‑CCCG‑GGGGGGGGGGGGGGGGA‑GGGGAATCACCCGGCGGCCCTTCCCGCCCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGa                                                                                 <  3:471063/151‑1 (MQ=255)                                                                           gAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACctct                                                                                                    >  3:247110/1‑126 (MQ=255)                                                                            aaGTACCCACATTGAGAACCCACACAGAAA‑CCCG‑GGGGGGGGGGGGGGGGA‑GGGGAATCCCCCGGCGGCCCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                          <  1:356943/151‑1 (MQ=255)                                                                                aaTCACATTGAAAAACCCCACCAAAA‑CCCA‑GGGGGGGGGGGGGGGGA‑TGGGGATCAACCAGGGGCGCTTCCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTc                                                                       <  3:569065/150‑1 (MQ=255)                                                                                    acaTTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTgg                                                                  >  2:158204/1‑151 (MQ=255)                                                                                       ttGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                               >  1:1014052/1‑151 (MQ=255)                                                                                       tggagaAACCCCCCCGAAA‑CCCA‑GGGGGTGGGGAGGGGGA‑GGGGGACCACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGcc                                                                <  2:548433/148‑1 (MQ=255)                                                                                        tGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTTGCCgg                                                              >  3:806463/1‑151 (MQ=255)                                                                                        tGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCgg                                                              >  1:370401/1‑151 (MQ=255)                                                                                          agaAACCCCACCGAAA‑CCCG‑GGGGGGGGGGGGGGGGA‑GGGGGAACACCCAGCGGCGCTTCCCGCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGgat                                                            <  4:122608/151‑1 (MQ=255)                                                                                          agaAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGgtg                                                                          >  1:844684/1‑137 (MQ=255)                                                                                            aacccccacCAAAA‑CCCG‑GGGGGGGGGGGGGGGGG‑GGGGAACCACCCGGGGGCGCTCCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCg                                                               <  1:804057/140‑1 (MQ=255)                                                                                            aaaCCACCCCAAAA‑CCCG‑GGGGGGGGGGAGGGGGA‑GGGGGACCACCCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGga                                                             <  3:77056/148‑1 (MQ=255)                                                                                               ccACACCGAAA‑CCCA‑GGGGGGGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                       <  2:307513/151‑1 (MQ=255)                                                                                               ccACACCAAAA‑CCCA‑GGGGGTGGGGGGGGGGA‑GGGGAATCAACCAGCGGCCCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCt                                                       <  4:273194/151‑1 (MQ=255)                                                                                                 acacCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTa                                                      >  3:1151895/1‑150 (MQ=255)                                                                                                  cacCAAAA‑CCCA‑GGGGGGGGGGGGGGGGA‑GGGGGACCACCCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaa                                                     <  3:1320103/150‑1 (MQ=255)                                                                                                     cGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaac                                                 >  1:1113272/1‑151 (MQ=255)                                                                                                      ccaaaacccgggggggggggggggggT‑GGGGGACCAACCGGCGGCCCTCCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTaacaa                                                  <  4:382367/139‑1 (MQ=255)                                                                                                      aaaa‑CCCA‑GGGGGGGGGGGGGGGGA‑TGGGGATCACCCGGCGGCCCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACg                                                <  1:618496/150‑1 (MQ=255)                                                                                                         aaacccgggggggggggggggggA‑GGGGGATCACCCGGCGGCGCTTCCCCCTCCCCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGc                                               <  1:121036/142‑1 (MQ=255)                                                                                                         a‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCa                                              >  1:1234219/1‑150 (MQ=255)                                                                                                                    ggggggggggggA‑GGGGGACCACCCAGCGGCGCTCCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGt                                    <  1:797495/149‑1 (MQ=255)                                                                                                                         ggggggga‑gggggaTCACCCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                               <  3:522376/149‑1 (MQ=255)                                                                                                                          aGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagc                               >  2:147083/1‑150 (MQ=255)                                                                                                                           gggggg‑ggggAATCACCCGGCGGCGTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTgcagca                              <  3:803132/150‑1 (MQ=255)                                                                                                                            gggga‑tggggatACCCCGGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAg                            <  1:1066802/151‑1 (MQ=255)                                                                                                                            gggga‑gggggatCACCCAGCGGCGCTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAg                            <  2:256039/151‑1 (MQ=255)                                                                                                                             ggga‑gggggacCAACCGGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGt                           <  1:208840/151‑1 (MQ=255)                                                                                                                               ga‑tggggatCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGttt                         >  1:750991/1‑151 (MQ=255)                                                                                                                                   ggggatCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTcc                      >  1:294301/1‑151 (MQ=255)                                                                                                                                   ggggacCACCCAGCGGCGCTTCCCGCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTc                       <  3:69180/150‑1 (MQ=255)                                                                                                                                       atCACCCGGCGGCGCTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgc                   <  1:710790/150‑1 (MQ=255)                                                                                                                                       atCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgcg                  >  3:208629/1‑151 (MQ=255)                                                                                                                                          acccggcggccTTTCCCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                <  3:792095/145‑1 (MQ=255)                                                                                                                                          aaCCAGCGGCGTTTCCCCCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCa                <  2:35527/150‑1 (MQ=255)                                                                                                                                            ccggcgccgCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCaa               <  1:584435/146‑1 (MQ=255)                                                                                                                                            ccAGCGGCGCTTCCCCCCCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAgc                   <  1:993204/145‑1 (MQ=255)                                                                                                                                             cAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATgg            >  3:1179813/1‑151 (MQ=255)                                                                                                                                              aGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATggg           >  2:311043/1‑151 (MQ=255)                                                                                                                                               gcggcgCTTCCCGCACACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGc          <  3:748352/151‑1 (MQ=255)                                                                                                                                                      ttACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGc   >  1:31158/1‑151 (MQ=255)                                                                                                                                                       tACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCg  >  4:68324/1‑151 (MQ=255)                                                                                                                                                        |                                                                                                                                                      GTAGCGCGGTATCATCCGCCAGCAGATCGTGTTTCAGCCCCTGAATCAACGCCCTGGCGGTGGTGGGCGGTACGGAAGTAATCACATTGAGAAACCACACCGAAA‑TCCA‑GCGGGTGGGGAGGGGGA‑TGGGGATCAACCAGCGGCGCTTACCGCTCACCGCCATAAAATGTTCAAACTGTTGCTGATAACTGAGCACCTCTGGTCCGGCGGCTTCGAAGATGCGGTGTTCGCTGGCCGGATGATCTAACAACGCAACCAGATAGTGCAGCAAGTTTTCCAGCGCAATGGGCGTGGTGCG  >  NC_000913/907454‑907751

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