Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,847,971	0	A	C	46.0%	111.2 / 12.6	74	G218G (GGT→GGG)	uhpC	membrane protein regulates uhpT expression
Reads supporting (aligned to +/- strand):  ref base A (38/2);  new base C (0/34);  total (38/36)
Fisher's exact test for biased strand distribution p-value = 4.59e-19
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCGCCCGGACCACATAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGAATG  >  NC_000913/3847822‑3848118                                                                                                                                                        |                                                                                                                                                    gccgccCGGACCACATAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCac                                                                                                                                                    >  2:114928/1‑151 (MQ=255) gccgccCGGACCACATAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCac                                                                                                                                                    >  2:889149/1‑151 (MQ=255)      ccGGACCACATAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCaccgacc                                                                                                                                               >  2:1135756/1‑151 (MQ=255)           ccACATAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCgg                                                                                                                                          >  1:454068/1‑151 (MQ=255)              caTAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGATGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTa                                                                                                                                        >  1:154368/1‑150 (MQ=255)                 ataCCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAAcc                                                                                                                                    >  1:1201994/3‑151 (MQ=255)                 aGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAAcc                                                                                                                                    >  2:299994/1‑151 (MQ=255)                    ccAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTa                                                                                                                                  >  2:817281/1‑150 (MQ=255)                     cAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAAc                                                                                                                                >  1:70861/1‑151 (MQ=255)                      aGCACATAGCAAAAC‑GAAAGCAGCCAGCTATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTATACCTAACg                                                                                                                               >  2:389750/1‑151 (MQ=255)                        cacaTAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACg                                                                                                                               >  2:450985/1‑149 (MQ=255)                              agaaaaCGGAAAGCAGCCAGATATCGGGTTTCACCAACACATATTTGGGGAGGATCTTTTTCCGCGCCAC‑CCCCCCCCCTTCTTGTTTTTGAGAATTTCCCGGCGGGTCGTGTCGCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCgg                                                                                                                        <  1:559007/149‑1 (MQ=255)                                aaaaC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAC‑CCCTGCCCCTTCTTGTTGTGGAGCATTTTCCAGCGCGTGGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCgg                                                                                                                        <  2:96009/148‑1 (MQ=255)                                aaaaC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCg                                                                                                                     >  2:244203/1‑151 (MQ=255)                                aaaaC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAAGATCTCTTAACGAGTCAA‑CCCTGCCCCTTCTTGTTGTTAAGCAATTTCCAGCGCGTCGTATCGACATTCACCGACCGCCGGTAAACCTAACGCCTCCggac                                                                                                                      >  1:1006265/1‑148 (MQ=255)                                aaaaC‑GAAAGCAGCCACATATACGGATTCAGCAACACATTTTTGTTGAGGATCTCTTTCCGCGTCAC‑CCCTCCCCTTTTTTTTTGTTGAGCACTTTCCAGCGCGTCGTGTCACCATTCCCCGACCGCCGGTAAACCTAACGCCTGCg                                                                                                                         <  2:1013945/147‑1 (MQ=255)                                  aaC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAC‑CCCTGCCCCTTCTTGTTGTTGAGCAATTCCCAGCGCGTCGTGTCGCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGAt                                                                                                                   <  2:611895/151‑1 (MQ=255)                                   aagaaaaGCAGCCAGATATACGGATTCAGCAACACATATTGGGTGAGGTTCTCTTTCCGCGTCAA‑CCCTGCCCCTTTTTGTTTTTGAGCAATTCCCAGCGCGTCGGGTCGCCTTTCCCCAACCGCCGGTAAACCTAACGCCTGCGGGCGAt                                                                                                                   <  2:634257/147‑1 (MQ=255)                                     caaaaGCACCCGGATAAAGGAATTCAGCAACATTTATTGGGGGAGCTTCTCTTTCCGCGTCAC‑CCCTGCCCCTTCTTGTTGTTGAGCATTTCCGGGGGCTGGGTGCCGCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATc                                                                                                                  <  1:1219660/148‑1 (MQ=255)                                         agcagcCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGt                                                                                                              >  2:895626/1‑150 (MQ=255)                                           cagcCAGATATACGGATTCAGCAACACTTATTTGGTGAGGATCTCTTTACGCGTCAC‑CCCGGCCCCTTCTTGTTGTGGAGCAATTCCCAGCGCGTGGTGCCGCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAg                                                                                                            <  1:225876/150‑1 (MQ=255)                                           cagcCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTCCGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTCCCAGCGCGTTGGGTCCCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAg                                                                                                            <  2:76822/150‑1 (MQ=255)                                              ccAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTTTTTCCGCGTCAA‑CCCGCCCCCTTCTTGTTGTTGAGCAATTCCCAGCGCGTCGTGTCGCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAgcc                                                                                                          <  2:968378/149‑1 (MQ=255)                                              ccAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAgccgc                                                                                                        >  2:895774/1‑151 (MQ=255)                                              ccAAATAAACGGATTCGGCAACACATATTGGGTGGGTTTCTCTTT‑CCCCTCAATCCCGCCCCCTTCTTGTTGTGGAGCATTTCCCGCCGCGGCGTTCCCCCTTTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAgcc                                                                                                          <  2:770809/149‑1 (MQ=255)                                              caaGATATACGGTTCAAGCAACACATATTTGGTGAGGATCTTTTTCCGCGTCAC‑CCCCCCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCTTTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAgccg                                                                                                         <  2:909148/148‑1 (MQ=255)                                                aGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAgccgcc                                                                                                       >  2:361634/1‑150 (MQ=255)                                                 gATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAgccgcc                                                                                                       >  2:240866/1‑149 (MQ=255)                                                   tataCGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCcagc                                                                                                    >  2:269803/1‑150 (MQ=255)                                                       cGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAC‑CCCGCCCCCTTCTTTTTGTTGAGCAATTCCCAGCGCGTCGTGTCGCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCagag                                                                                                <  2:156396/150‑1 (MQ=255)                                                       cGGATTAAGCAACAATTTTTTGGTGAGGTTCTTTTTACGCGTCAC‑CCCTGCCCTTTCTTGTTGTTGAAAATTTCCCGGCGCGTGGTGTCCCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCcagcag                                                                                                  <  1:1116972/148‑1 (MQ=255)                                                        ggATTCAACAACAAATTTTGGGTGGGGTTTTTTTTACGCGTAAC‑CCCTCCCCCTTCTTGTTGTTGAGCAATTCCCAGCGGGTCGTGTCGCATTTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCaga                                                                                                 <  1:114927/148‑1 (MQ=255)                                                           ttCAGCAACACATATTTGGTGGGGATCTCTTTACGCGTCAC‑CCCTGCCCCTTCTTTTTGTTGAGCAATTCCAGGCGCGTCGTGTCCCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGaaa                                                                                             <  2:601542/149‑1 (MQ=255)                                                           ttCACCACCACATATTTGGTGGGGTTTTCTTTACGCGTCAC‑CCCGGCCCCTTCTTGTTGTTGAGAAATTTCCAGCGCGTCGGGTCGCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGaaa                                                                                             <  2:460537/149‑1 (MQ=255)                                                           tccggCAACACATATTTGGTGAGGATTTCTTTCCGCGTCAA‑CCCGGCCCCTTTTTGTTGTGGAGCAATTCCCAGCGCGTCGTGTCGCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCagag                                                                                                <  1:264332/142‑1 (MQ=255)                                                             caccaacaaATATTTGGGGAGGTTCTCTTTACGCGTCAC‑CCCTCCCCCTTCTTGTTTTTGAGAATTTCCCGGCGGGTCGTGCCCCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGaaaaa                                                                                           <  1:514672/146‑1 (MQ=255)                                                              aGCAACACATTTTTGGTGAGGATCTCTTTACGCGTCAC‑CCCTCCCCCTTCTTGTTGTGGAGAATTTTCCAGCGCGTCGTGTCGCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATc                                                                                         <  2:757492/150‑1 (MQ=255)                                                               gCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGCAAAATccc                                                                                       >  2:464830/1‑151 (MQ=255)                                                                   catattttgggggAGGTTTTCTTGACGCGCCAC‑CCCGCCCCCTTGTTTTTGTTGGGCATTTCCCAGCGGGTCGTGCCGCCATCCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACg                                                                                   <  1:86949/139‑1 (MQ=255)                                                                   cacaTATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACg                                                                                   >  2:868650/1‑151 (MQ=255)                                                                     caTATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGAc                                                                                 >  1:400649/1‑151 (MQ=255)                                                                      atatTTGGGAGGGGTCTCTTCCCGCCACAA‑CCCGCCCCCTCCTTGTTGTGGGCCATGCCCCGGCGCGTCGTACCGACATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGAc                                                                                 <  2:1111984/150‑1 (MQ=255)                                                                       tatTTGGGGAGGATCTTTTTCCGCGTCAC‑CCCTCCCCCTTTTTTTTGTTGAACAATTCCCAGCGCGTCGTGCCGCCTTTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATcccc                                                                                      <  1:1206760/144‑1 (MQ=255)                                                                         tttGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTa                                                                               >  1:242825/1‑149 (MQ=255)                                                                            ggTGAGGTTCTCTTTACGCGTCAC‑CCCGCCCCCTTTTTGTTGTTGAGCATTTCCCGGCGCGTCGTGCCCCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGcc                                                                          <  1:497166/151‑1 (MQ=255)                                                                                aGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATAc                                                                      >  2:201323/1‑151 (MQ=255)                                                                                      tctTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAAcc                                                                  >  2:1008344/1‑149 (MQ=255)                                                                                        tttGCGCGCCAC‑CCCTCCCCTTTTTTTTTTTGGAGAAATTCCCGCCGCGTCTGGTCCCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATa                                                                       <  1:392420/142‑1 (MQ=255)                                                                                        tttACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGc                                                               >  1:646390/1‑150 (MQ=255)                                                                                           tccgtctcacccttcccccTTCTTTTTGTTAAGCAATTCCCAGCGCGCCGTGCCCCAATTCCCCGCCCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCaa                                                             <  1:450964/135‑1 (MQ=255)                                                                                                tCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAatcatcat                                                      >  2:924534/1‑151 (MQ=255)                                                                                                tCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAatcatcat                                                      >  1:829248/1‑151 (MQ=255)                                                                                                 cAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACACCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAatcatcatc                                                     >  2:1214148/1‑151 (MQ=255)                                                                                                    ccccTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCCCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATccc                                                   <  2:1084924/150‑1 (MQ=255)                                                                                                     cccTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGCAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCg                                                  >  2:825504/1‑151 (MQ=255)                                                                                                          ccccTTTTTGTTGTTGAGCAATTCCCAGCGCTGCGGGTCCCCATTCCCCGGCCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGc                                                <  2:1176012/148‑1 (MQ=255)                                                                                                          ccccTTTTGTTTGTTGAGCAATCACCAGCGTGGCGTCTCCCCTCCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACg                                             <  1:1013377/151‑1 (MQ=255)                                                                                                          ccccTTCTTGTTTTTGAGCAATTCCCGGCGGTTCGTGCCGCCATCCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCa                                               <  1:554365/149‑1 (MQ=255)                                                                                                            ccTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGcc                                           >  2:49167/1‑151 (MQ=255)                                                                                                             cttcttGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGcc                                           >  2:878057/1‑150 (MQ=255)                                                                                                                  tgttgttgAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGt                                     >  2:752919/1‑151 (MQ=255)                                                                                                                       ttaaccAATTCCCAGCGCTTCGTGTCCCAATCCCCCGCCCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATg                                 <  2:421140/145‑1 (MQ=255)                                                                                                                         gAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATgc                                >  2:562974/1‑149 (MQ=255)                                                                                                                             atttcccggcggGTCGTGCCGCCATCCCCCAACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATgc                                <  1:483112/137‑1 (MQ=255)                                                                                                                               ttcccggcgcGTCGTGCCCCCATCCCCCGCCCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATgcagcgc                           <  1:361616/142‑1 (MQ=255)                                                                                                                                  caggcgcttggcgccGCAATTCCCGGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGc                     <  2:301787/137‑1 (MQ=255)                                                                                                                                     gcgcGTCGTGTCCCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGc                  <  2:248390/151‑1 (MQ=255)                                                                                                                                            gtgtCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATcac            >  1:288758/1‑150 (MQ=255)                                                                                                                                             tgtCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATcaca           >  1:108569/1‑150 (MQ=255)                                                                                                                                             tgtCCCCATTCCCCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATca             <  2:245988/148‑1 (MQ=255)                                                                                                                                                cGCCGTTCCCCGCCCGCCGGTAACCCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATc              <  1:428850/144‑1 (MQ=255)                                                                                                                                                 gCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATggg      >  1:750814/1‑151 (MQ=255)                                                                                                                                                   cATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGa     >  2:1200034/1‑150 (MQ=255)                                                                                                                                                   cATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGa     >  2:247691/1‑150 (MQ=255)                                                                                                                                                   cATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGa     >  2:526854/1‑150 (MQ=255)                                                                                                                                                       cACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGAATg  >  2:828235/1‑149 (MQ=255)                                                                                                                                                       cACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGAATg  >  1:297346/1‑149 (MQ=255)                                                                                                                                                        |                                                                                                                                                    GCCGCCCGGACCACATAGACCAGCACATAGCAAAAC‑GAAAGCAGCCAGATATACGGATTCAGCAACACATATTTGGTGAGGATCTCTTTACGCGTCAA‑CCCTGCCCCTTCTTGTTGTTGAGCAATTTCCAGCGCGTCGTGTCGCCATTCACCGACCGCCGGTAAACCTAACGCCTGCGGGCGATCGCGTAGCCGCCAGCAGAGAAAAATCCCCACGACTATCGCCATACAACCAGCAATCATCATCCCGGCACGCCAGCCGTAATGCAGCGCAGCCGCTGCCATCACAATGGGAATG  >  NC_000913/3847822‑3848118

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