Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,578,143	0	A	C	40.5%	136.8 / 21.7	79	L79V (TTA→GTA)	maeB	malic enzyme: putative oxidoreductase/phosphotransacetylase
Reads supporting (aligned to +/- strand):  ref base A (39/8);  new base C (0/32);  total (39/40)
Fisher's exact test for biased strand distribution p-value = 7.02e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGT  >  NC_000913/2577998‑2578285                                                                                                                                                   |                                                                                                                                                  tCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑cccc                                                                                                                                              <  2:693081/150‑1 (MQ=255)     gcgcGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑c                                                                                                                                          <  1:1179597/149‑1 (MQ=255)        cggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑ccgc                                                                                                                                    <  2:502704/151‑1 (MQ=255)         ggcgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑ccgcc                                                                                                                                   <  1:465733/151‑1 (MQ=255)           cgACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGt                                                                                                                                 >  2:905177/1‑151 (MQ=255)             aCAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTAc                                                                                                                               <  2:31833/151‑1 (MQ=255)              cAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTAcc                                                                                                                              >  2:1006406/1‑151 (MQ=255)                     aaTAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTag                                                                                                                         >  1:1222565/1‑149 (MQ=255)                       tAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTag                                                                                                                         >  2:181474/1‑147 (MQ=255)                           tttGTCCGGGTCGAGTTCGTCAACTTCAATGTAAAATACATCAACCCCGGCGAATTTCTAAAAC‑AGACCGCCCTTGCCTCCCATCCCCGGTTTCCCGGCCAGCGCGCCAATGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCAc                                                                                                                  <  2:106447/150‑1 (MQ=255)                             tGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCg                                                                                                                >  1:528753/1‑150 (MQ=255)                              gTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGcc                                                                                                              >  1:185909/1‑151 (MQ=255)                                 cGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGccacc                                                                                                           >  1:431224/1‑151 (MQ=255)                                  gggTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTTCCGTTAGAGATCACCGccac                                                                                                            >  1:558762/1‑149 (MQ=255)                                  gggTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGccacc                                                                                                           >  1:666705/1‑150 (MQ=255)                                    gTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAgg                                                                                                        >  1:279798/1‑151 (MQ=255)                                       gAGTTCTTCAACTTCAATGTAAAATACATCAACCCCGGCGATTTTCTAAAAC‑AGACCCCCCTTCCCTCCCACCCCCGGTTTGCCTGCCAGCGCGCAATTGTTCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGtt                                                                                                      <  2:597716/150‑1 (MQ=255)                                         gTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTAcc                                                                                                   >  1:453945/1‑151 (MQ=255)                                          ttcttcAACTCCAATGTAAAATACATCAACCCCGGCAAATTTCTTAAAC‑AGACCCCCCTTCCCTTCCACCCCCGGTTGCCCTGCCGGCGCGCAAATGTTCCCTAC‑CCCCCG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTa                                                                                                     <  2:896368/148‑1 (MQ=255)                                            cGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCg                                                                                                >  2:1236471/1‑151 (MQ=255)                                              tCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCggg                                                                                              >  1:1259033/1‑151 (MQ=255)                                                accTTCATTGTAAAATTCATCAACCCCGCCGATTTTCTTAAAC‑AGAACCCCCTTCCTTCCCATCCCCGTTTCGCTGGCCGCGGCACCACTGTTGCCCAC‑CCCGAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCt                                                                                                  <  2:744334/143‑1 (MQ=255)                                                aaCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCg                                                                                            >  2:258168/1‑151 (MQ=255)                                                aaCTTCAATGTAAAATACATAACCCCCGGCAAATTTCTAAAAC‑AGACCCCCTTTCCCTCCCATCCCCGTTTTCCCTCCCAGCGCGCCAATGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCgg                                                                                               <  2:201622/148‑1 (MQ=255)                                                 aCTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCgg                                                                                           >  1:807206/1‑151 (MQ=255)                                                   ttCAATGTAAAATCCACCAACCCCGGAGATTTTCTAAAAC‑AGACCGCCCTT‑CCCTCCCCCCCCCGTTTTCCCGCCCGGCCCCCAATTGTGCCCAACCCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGc                                                                                             <  1:111809/147‑1 (MQ=255)                                                      aatttaaaaTACACCAACCCCGGCGAATTTCTTAAAC‑AGACCCCCCTTCCCTCCCACCCCCGGTTTGCCTCCCAGCGCGCCAATGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGtat                                                                                        <  2:806692/143‑1 (MQ=255)                                                      aaTGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTatat                                                                                      >  2:24666/1‑151 (MQ=255)                                                      aaTGTCAAATACATCAACCCCGGCGAATTTTTAAAAC‑AGACCCCCCTTCCCTTCCACCCCCGTTTTCCCCCCCAGCGCCCCAATGTTGCCTAA‑CCCCAG‑CG‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGta                                                                                         <  1:124591/148‑1 (MQ=255)                                                      aaTGTAAAATACATCAACCCGGGCAATTTTCTAAAAC‑AGACCCCCTTTCCCTCCCATCCCGGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGta                                                                                         <  2:499669/148‑1 (MQ=255)                                                       aTGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATAtt                                                                                     >  1:1261778/1‑151 (MQ=255)                                                            aaaTACATCACCCCCGGCGATTTTCTAAAAC‑AAACCCCCCTCCCCTCCCACCCCCGGTTTGCCTGCCAGCGCGCAATTGTTGCCAAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAg                                                                                <  2:1185822/151‑1 (MQ=255)                                                            aaaTACATCACCCCCGGCAATTTTCTAAAAC‑AGACCCCCCTTCCCTCCCACCCCCGTTTTCCCTGCCAGCGCCCAATTGTTCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTg                                                                                   <  1:1165022/148‑1 (MQ=255)                                                            aaaTACATAAACCCCGGCAATTTTCTTAAAC‑AGACCCCCCTTCCCTTCCACCCCCGGTTTGCCTGCCAGCGCCCCAATGTTGCCAAC‑CCCCAG‑CC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTa                                                                                 <  1:434975/150‑1 (MQ=255)                                                               tACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCt                                                                             >  2:505305/1‑151 (MQ=255)                                                                 cATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCttt                                                                           >  1:882921/1‑151 (MQ=255)                                                                  aTCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                          >  1:839411/1‑151 (MQ=255)                                                                  aTCAACCCCGCCAATTTCCTAAAAC‑AGAACCCCCTTCCCTCCCACCCCCTGTTCCCCTCCCGGCGCCCCAATGTTCCCTAC‑CCCCAC‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                            <  2:764463/149‑1 (MQ=255)                                                                   tCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtttt                                                                          >  2:1245785/1‑150 (MQ=255)                                                                     aaccccGGCGATTTTCTAAAAC‑AGACCCCCCTTGCCTCCCATCCCCGGTTGGCCTCCCAGCGCCCAAATGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCtt                                                                            <  2:1111929/143‑1 (MQ=255)                                                                     aaTCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAAc                                                                       >  1:559067/1‑151 (MQ=255)                                                                      aTCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                      >  2:760354/1‑151 (MQ=255)                                                                       tCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACgg                                                                     >  2:210760/1‑151 (MQ=255)                                                                        ccgggCGATTTTCTAAAAC‑AAACCCCCTTTGCCTCCCATCCCCGTTTTGCCTGCCGGCGCGCAATTGTTCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                      <  2:235882/146‑1 (MQ=255)                                                                        cccGGCAATTTTCTAAAAC‑AGACCGCCTTTGCCTTCCATCCCCGGTTTGCCTGCCAGCGCGCAATTGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACg                                                                      <  2:263981/149‑1 (MQ=255)                                                                        cccGGCAATTTTCTAAAAC‑AGACCGCCCTTGCCTCCCACCCCCGGTTTGCCTGCCAGCGCGCCAATGTTCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGc                                                                              <  2:874249/141‑1 (MQ=255)                                                                           ggCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                                 >  1:575183/1‑151 (MQ=255)                                                                            gaaatttccaaaaac‑aaacCGCCTTTGCCTCCCACCCCGGTTTTCCCGGCCGGCGCGCAATTTTTGCCAAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                <  1:804200/140‑1 (MQ=255)                                                                            gCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCtt                                                                >  2:373042/1‑151 (MQ=255)                                                                               aaTTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCttttt                                                             >  1:436195/1‑151 (MQ=255)                                                                                aTTTCCTAAACGAGACCCCCCTTCCCTTCCACCCCCGTTTTGCCGCCCGGCGCCCCACTTTTGCCCAC‑CCCCGG‑CC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCt                                                                 <  1:949247/147‑1 (MQ=255)                                                                                 tttCTTAAAC‑AGACCCCCCTTCCCTTCCATCCCCGGTTTCCCTGCCAGCGCCCAATTGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTc                                                            <  2:1020859/150‑1 (MQ=255)                                                                                  ttCTTAAAA‑AGACCGCCCTTCCCTCCCACCCCCGGTTTGCCTGCCAGCGCGCAAATGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCg                                                           <  2:1099183/150‑1 (MQ=255)                                                                                       aaaC‑AGACCGCCCTTGCTTCCCACCCCCGTTTTGCCTGCCAGCGCGCCATTGTTCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAtt                                                        <  2:234705/148‑1 (MQ=255)                                                                                       aaaC‑AGAACCCCCTTGCCTCCCATCACCGGTTTGCCTGCCAGCGCGCCATTTTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGAttt                                                       <  2:2349/149‑1 (MQ=255)                                                                                               aCGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                >  1:778862/1‑149 (MQ=255)                                                                                                cccctgtgcttcccATCCCGGGTTTGCCTCCCGGCGCGCCAATGTTCCCTAC‑CCCCAG‑CC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACa                                                <  1:765602/136‑1 (MQ=255)                                                                                                 ccccTTCCCTCCCACCCCCGTTTGGCCGGCCAGCCCGCATATGTGCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                            <  1:1233405/150‑1 (MQ=255)                                                                                                  cctttGCCTCCCATCCCCGGTTTGCCTGCCGGCGCGCAAATTTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTg                                           <  2:853556/148‑1 (MQ=255)                                                                                                  cccTTCCCTCCCATCCCCGTTTTGCCCGCCAGCGCCCCATTTTTCCCTAC‑CCCCAC‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACaa                                               <  1:229646/147‑1 (MQ=255)                                                                                                  cccTTCCCTCCCACCCCCGGTTTGCCTGCCACCGCGCCAATGTTGCCTAC‑CCCCAC‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGt                                            <  2:1219323/150‑1 (MQ=255)                                                                                                   ccTTGCTTTCCACCCCCGGTTTCCCTGCCGGCCCCCCAATGTTCCTTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                          <  2:968290/151‑1 (MQ=255)                                                                                                   ccTTCCCTCCCACCCCCGTTTTGCCTGCCGGCGCGCAAATTTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                          <  2:591671/151‑1 (MQ=255)                                                                                                   ccTTCCCTCCCACCCCCGTTTTGCCTGCCAGCGCCCAATTGTTGCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGc                                          <  1:245443/151‑1 (MQ=255)                                                                                                     ttGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCgg                                        >  2:421008/1‑151 (MQ=255)                                                                                                        cttcccacccccGTTTTGCCTCCCCGCGCCCCATTTTTCCAAAC‑CCCCAC‑CC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                      <  2:1207729/140‑1 (MQ=255)                                                                                                        ccTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCa                                      >  1:615792/1‑150 (MQ=255)                                                                                                        ccTCCCACCCCCGGTTTGCCTGCGAGCGCCCAATTGTTCCCT‑A‑CCCCAA‑CATCCGCCGTACCGTTAGAGATCGCCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCaa                                     <  1:1121293/151‑1 (MQ=255)                                                                                                         cTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                    >  2:797918/1‑151 (MQ=255)                                                                                                          ttCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACg                                   >  1:1238834/1‑151 (MQ=255)                                                                                                           tCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGc                                  >  1:244973/1‑151 (MQ=255)                                                                                                            ccATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGcc                                 >  1:548607/1‑151 (MQ=255)                                                                                                               tCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTgg                              >  2:570905/1‑151 (MQ=255)                                                                                                               tCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTgg                              >  2:189673/1‑151 (MQ=255)                                                                                                                ccccggttggctgcccAGCGCGCAATTTTTCCCTAC‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAAc                                    <  1:1188772/130‑1 (MQ=255)                                                                                                                 cccGTTTTGCCTGCCGCCGCCCCATGGTTGCCTAC‑CCCCAG‑CC‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTg                            <  1:947505/150‑1 (MQ=255)                                                                                                                               cAGCGCGCCGATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc              >  1:418308/1‑151 (MQ=255)                                                                                                                                   gcgcCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGcc              >  2:367821/1‑147 (MQ=255)                                                                                                                                       cAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATcgcg      >  1:788451/1‑151 (MQ=255)                                                                                                                                             ttcccaa‑CCCCAACCA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCtgt  <  1:799308/145‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                  TCGAGCGCGGCGACAACTTCAATAAATTTGTCCGGGTCGAGTTCGTCAACTTCAATGTCAAATACATCAATCCCGGCGAATTTCTTAAAC‑AGAACGCCCTTGCCTTCCATCACCGGTTTGCCTGCCAGCGCGCCAATGTTGCCTAA‑CCCCAG‑CA‑CCGCCGTACCGTTAGAGATCACCGCCACCAGGTTACCTCGGGCGGTATATTTGTAGGCTTTTAACGGGTCTTTTTCGATTTCAAGACAAGGTGCGGCAACGCCTGGTGAGTAGGCCAGCGCCAGATCGCGCTGT  >  NC_000913/2577998‑2578285

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