Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,386,862	0	G	A	42.9%	80.2 / 62.5	77	intergenic (+31/‑72)	yfbL/yfbM	putative M28A family peptidase/DUF1877 family protein
Reads supporting (aligned to +/- strand):  ref base G (7/37);  new base A (33/0);  total (40/37)
Fisher's exact test for biased strand distribution p-value = 5.75e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  NC_000913/2386713‑2387011                                                                                                                                                      |                                                                                                                                                          ccGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAg                                                                                                                                                           <  2:815623/150‑1 (MQ=255) ccGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGc                                                                                                                                                          <  2:248175/151‑1 (MQ=255)       tttATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACAt                                                                                                                                                     <  1:1171052/150‑1 (MQ=255)               aataaACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACATggccaattt                                                                                                                                            >  1:989895/1‑143 (MQ=255)                       aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGCAAAAAAGACATGGCC‑ctttttttttt                                                                                                                                     >  1:339146/1‑148 (MQ=255)                       aTACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGCCCCAAAAAAAAGACACGACA‑attttttttt                                                                                                                                      >  1:159/1‑148 (MQ=255)                          ccACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTTGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCa                                                                                                                                <  2:693415/151‑1 (MQ=255)                          ccACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCt                                                                                                                                  <  2:940032/149‑1 (MQ=255)                          ccACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCa                                                                                                                                <  1:578611/151‑1 (MQ=255)                            aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑t                                                                                                                              <  1:134014/150‑1 (MQ=255)                            aCTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAGGAC‑ATTTTTTTTTtcccc                                                                                                                                 >  2:396266/1‑146 (MQ=255)                                gCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGt                                                                                                                          <  1:277227/150‑1 (MQ=255)                                 cccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGTAC‑ATTTTATTTTTCCTCA‑TCTGttt                                                                                                                        >  1:551032/1‑151 (MQ=255)                                  ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCAAAAAAAAAGAAATGCACAATTTTTTTTTCCCTC‑‑TCTGtttt                                                                                                                       >  1:3406/1‑151 (MQ=255)                                  ccGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCTAAAAAAAGAAAAGCAC‑AATTTTTTTTTTCTCA‑TCTGtttt                                                                                                                       >  2:150596/1‑151 (MQ=255)                                    ggTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACAAGCAC‑ATTTTTTTTTTCCTCA‑TCTGtttt                                                                                                                       >  1:60178/1‑149 (MQ=255)                                    ggTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGGCAGAAAAAAAGCAAAGGAC‑AATTTTTTTTTCCCCA‑TCTGgtt                                                                                                                        >  2:871204/1‑148 (MQ=255)                                        aTACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATg                                                                                                                 <  1:42225/151‑1 (MQ=255)                                         tACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAACGACAAGCAC‑ATTTTTTTTTTCCTCA‑TCtttttt                                                                                                                       >  1:853392/1‑144 (MQ=255)                                           ccGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑Ga                                                                                                             <  1:888126/151‑1 (MQ=255)                                                    gaTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑aaa                                                                                                   <  1:363784/151‑1 (MQ=255)                                                     aTTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAAt                                                                                                  <  2:937399/151‑1 (MQ=255)                                                      ttGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCAC‑ATTTTTTTTTTCCTCA‑TCTTTTTTAAAAAGC‑GGAACATG‑aaa                                                                                                   >  1:622444/1‑149 (MQ=255)                                                       tGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGAAATGCAC‑ATTTTTTTTTTCCTCC‑TCTTTTTTTAAATGG‑GAAAAAAG‑aaaat                                                                                                 >  1:839463/1‑148 (MQ=255)                                                        gAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAAAGACAAGCAC‑AATTTTTTTTTCCTCC‑TCTGTTTTTAGATGG‑GAAAAAAG‑aaa                                                                                                   >  1:993947/1‑147 (MQ=255)                                                          aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAAc                                                                                             <  1:1094550/151‑1 (MQ=255)                                                          aTTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAAc                                                                                             <  2:386279/151‑1 (MQ=255)                                                           ttATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCACAAAAAACGACATGAAC‑ATTTTTTTTTTCCTCATTCTTTTTTAA‑ATGC‑GGAAAATG‑aaaattaa                                                                                              >  1:557522/1‑144 (MQ=255)                                                               cAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACggagg                                                                                        <  2:146378/151‑1 (MQ=255)                                                                aGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACggagga                                                                                       <  1:780681/151‑1 (MQ=255)                                                                  aaaaTGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGaaa                                                                                     <  1:309588/151‑1 (MQ=255)                                                                         cTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTTCCGAAAAAAAGACAATCAC‑AATTTTTTTTTTCCCA‑TCTGTTTTTAGATGG‑GAAAAAAG‑AAAATTAACAGAGGGAaaataaa                                                                               >  2:1047311/1‑150 (MQ=255)                                                                         cTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGAAAATCAA‑AATTTTTTTTTCCTCC‑TCTGTTTTTAAATTG‑GAAAAAAGAAAATTCAAAGG‑GGAaaaaaaaa                                                                               >  2:483857/1‑150 (MQ=255)                                                                          tCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCTGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATg                                                                             >  2:1144201/1‑151 (MQ=255)                                                                           cAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAaaataaat                                                                              <  1:923013/149‑1 (MQ=255)                                                                               tAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCAC‑AATTTTTTTTTCCTCC‑TCTGTTTTTAGAAGCGGAAAAATA‑AAATTAAAAGAGGGAAAAAAAAATg                                                                             >  2:696544/1‑147 (MQ=255)                                                                                aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTGCAGAAAAAACGAAAAGAAC‑AATTTTTTTTTCCTCC‑TCTGTTTTTATATAG‑GAAAAACA‑AAAATTAACGGAGGGAaaata                                                                                 >  2:796221/1‑141 (MQ=255)                                                                                aGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAAt                                                                        <  2:417104/150‑1 (MQ=255)                                                                                 gTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGa                                                                      <  2:1063669/151‑1 (MQ=255)                                                                                  tGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGACATACAC‑ATTTTTTTTTTCCTCA‑TCTGTTTTTAAATGCGGAAAAATA‑AATTTTAAGGGGGGAAAAAAAAAAGGGa                                                                          >  1:23763/1‑147 (MQ=255)                                                                                   ggATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGAt                                                                     <  1:584466/150‑1 (MQ=255)                                                                                     atcgtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGTCAGAAAAAAAGACAAACAC‑ATTTTTTTTTTCCTTC‑TTTGTTTTTAAAAGG‑GAAACATA‑AAATTTAAAAGAGGAAAAAAAAATagg                                                                           >  2:836617/4‑142 (MQ=255)                                                                                       ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                                <  1:118741/151‑1 (MQ=255)                                                                                       ggTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACAAGCAC‑AATTTTTTTTTCCTCC‑TCTGTTTTTAAAAGC‑GAAACATG‑AAATTTAACGGGGGGAAAAAAAATAGGGAAGAtttg                                                                  >  1:243272/1‑147 (MQ=255)                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCt                                                                <  1:755079/150‑1 (MQ=255)                                                                                        gtgtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTa                                                               <  2:8309/151‑1 (MQ=255)                                                                                          gtTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGACAAGCAC‑AATTTTTTTTTCCTCA‑TCTGTTTTTAAAAGA‑GAAAAATG‑AAATTCAACGGAGGGAAAAAAAAATGGAATGATTTGGTAtt                                                             >  1:962646/1‑151 (MQ=255)                                                                                             ataACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTg                                                          <  2:1043713/151‑1 (MQ=255)                                                                                              taACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCTGAAAAAACGAAAAGCAC‑AATTTTTTTTTTCTCC‑TCTGTTTTTAGATGC‑GGAACATG‑AAATTTAAAGGAGGAAAAATAAATGGGAATGATTGGCTAtttt                                                           >  2:45389/1‑149 (MQ=255)                                                                                                   ttGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGaa                                                     <  1:171742/150‑1 (MQ=255)                                                                                                    tGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCCGAAAAAAAGAAAAGAAC‑AATTTTTTTTTTCTCC‑TCTTTTTTTAAATGC‑GAAAAAAG‑AAAATTAAAAGAGGAAAAAAAAAAAGGGATGAATTGCTATTTTTTAGaaa                                                    >  2:33052/1‑150 (MQ=255)                                                                                                     gTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAAt                                                   <  2:729904/150‑1 (MQ=255)                                                                                                       tataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAAAGAAAAACAC‑ATTTTTTTTTTCCTCA‑TCTGTTTTTAGATGA‑GGAACATG‑AAAATTAAAAGAGGGAAAAAAAAAGGGGAAGATTGGCTATTTTTTAGaaa                                                    >  2:35122/1‑147 (MQ=255)                                                                                                        ataCAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGa                                                <  2:661314/150‑1 (MQ=255)                                                                                                              caGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGaa                                         <  2:408243/151‑1 (MQ=255)                                                                                                               aGTAAATAATAAAAATTTCTCGAGGCGGGTTCAAAAAAAAAGACATGCAA‑ATTTTTTTTTTCCTCA‑TCTTTTTTTAAAAGC‑GAAAAATG‑AAATTTAAACGAGGAAAAAAAAAATGGGATGATTGGCTATTTTGGAGAAAATGAATCAGaa                                         >  1:47039/1‑150 (MQ=255)                                                                                                                 tAAATAATAAAAATTTCTCGAGGCGGGGCCCGAAAAAACGGAAAGCAA‑AATTTTTTTTTTCTCC‑TCTGTTTTTAGATGG‑GAAACAAT‑AAAATTAACAGAGGGAAAAAAAATGGGGATGATTGGCTATTTTTCAGAAAAATATTCAGaaaa                                       >  2:249734/1‑150 (MQ=255)                                                                                                                 tAAATAATAAAAATTTCTCGAGGCGGAGCCTAAAAAAACGGCATGCAC‑AATTTTTTTTGCCTTC‑TCTGTTTTTAAATGG‑GGAACAAG‑AAATTTTACGGAGGAGAAAAAAATGTGGAAGATTTGCGCTTTTTCTGAAAATGGATCAGaaaa                                       >  1:709657/1‑150 (MQ=255)                                                                                                                   aataatAAAAATTTCTCGAGGCGGGGCCCGAAAAAAAGACAAGAAC‑AATTTTTTTTTTCTCC‑TCTGTTTTTAGAAGC‑GAAAAAAGAAAATTAAAAGGGGGAAAAA‑AAATAGGGATGATTGGCTATTTTTGTGAAAATTAATCAGaaaaaa                                     >  2:1018407/1‑148 (MQ=255)                                                                                                                     taataaAAATTTCTCGAGGCGGGGCCAGAAAAAAAGACATGCAA‑AATTTTTTTTTTCTCC‑TCTGTTTTTAGAAGG‑GAAAAAAA‑AAATTTAACGGGGGAAAAAAAAATTGGGAAGATTTGCTATTTTGTAGAAAATGACACAGAAAAAAtt                                   >  2:1099143/1‑150 (MQ=255)                                                                                                                      aataaAAATTTCTCGAGGCGGGTCCAGAAAAAACGACATGGAC‑AATTTTTTTTTCCTCA‑TCTTTTTTTAGATGC‑GAAACACG‑AAAATCCACGGTGGAAAAAAAAAATGGGATGATTGGCTATTTTTTAGAAATTTAATCAGgaaa                                       >  2:530235/1‑145 (MQ=255)                                                                                                                         aaaaaTTTCTCGAGGCGGGTCCAGAAAAAAAGAAATGCAC‑ATTTTTTTTTTCCTCC‑TCTTTTTTTAGAAGG‑GAAAAAAA‑AAATTTAAAGGAGGAAAAAAAAATTGGAAAGATTGGGTATTTTTTAGAAAATGAATCaaaaaaaaat                                   >  2:536109/1‑142 (MQ=255)                                                                                                                          aaaaTTTCTCGAGGCGGGTCCCCAAAAAACGGCAAGCAC‑ATTTTTTTTTTCCTCC‑TCTTTTTTTAAATGC‑GAAAAATGAAAATTAAAAGAGGGAAAAAAAAATGGG‑ATGATTGGCTATTTTGTAGAAAATGAATTAGAAAAGAATAAACa                              >  2:538993/1‑150 (MQ=255)                                                                                                                                   cGAGGCTGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCtt                       <  1:1164364/148‑1 (MQ=255)                                                                                                                                   cGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTg                      <  2:727394/149‑1 (MQ=255)                                                                                                                                      ggCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCt                 <  2:1026067/151‑1 (MQ=255)                                                                                                                                      ggCGGGGCCAGAAAAAGCGACACGCAC‑AATTTATTTTCCCTCA‑TCTGTTTTAAGATGG‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTTACTCAGAAAAGATTAACCAACTGCTTGAGTc                  >  2:648423/1‑150 (MQ=255)                                                                                                                                       gCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCt                 >  1:317310/1‑150 (MQ=255)                                                                                                                                        cGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAc               <  2:460590/151‑1 (MQ=255)                                                                                                                                        cGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTAc               >  1:889411/1‑151 (MQ=255)                                                                                                                                          ggTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg             <  2:658749/151‑1 (MQ=255)                                                                                                                                          ggTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg             <  1:45392/151‑1 (MQ=255)                                                                                                                                          ggTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg             >  1:435509/1‑151 (MQ=255)                                                                                                                                          ggTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTg             >  2:878730/1‑151 (MQ=255)                                                                                                                                                gAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAAc       <  2:996246/151‑1 (MQ=255)                                                                                                                                                   aaaGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCtt    <  1:879978/151‑1 (MQ=255)                                                                                                                                                    aaGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCttt   >  2:1020541/1‑151 (MQ=255)                                                                                                                                                     aGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTg  <  1:120670/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                          CCGCTTTTTATCGTAATAAACAATACCACTTGCCCGGTGATACCGCAGACAGATTGAATTATCAGAAAATGGCTCAGGTAGTGGATGGTGTTATAACTTTGTTATACAACAGTAAATAATAAAAATTTCTCGAGGCGGGTCCAGAAAAAGCGACATGCAC‑ATTTTATTTTCCCTCA‑TCTGTTTTAAGATGC‑GAAACATG‑AAATTCAACGGAGGAAAAATAAATGGGAATGATTGGCTATTTTGCAGAAATTGACTCAGAAAAGATTAACCAACTGCTTGAGTCTACTGAAAAACCTTTG  >  NC_000913/2386713‑2387011

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