Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,153,573	0	T	G	43.5%	99.1 / 24.4	62	Y92D (TAC→GAC)	pabC	4‑amino‑4‑deoxychorismate lyase component of para‑aminobenzoate synthase multienzyme complex
Reads supporting (aligned to +/- strand):  ref base T (1/34);  new base G (27/0);  total (28/34)
Fisher's exact test for biased strand distribution p-value = 1.00e-16
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAA  >  NC_000913/1153426‑1153718                                                                                                                                                    |                                                                                                                                                    atCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑Ca                                                                                                                                                 <  2:1001736/151‑1 (MQ=255)        gACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGACGTGGGGA‑TAGcacact                                                                                                                                          >  1:812621/1‑149 (MQ=255)         aCTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGG‑GAGcacaatt                                                                                                                                         >  2:397630/1‑148 (MQ=255)            aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGaa                                                                                                                                       <  1:309815/150‑1 (MQ=255)            aCAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGaa                                                                                                                                       <  1:451332/150‑1 (MQ=255)              aGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGaa                                                                                                                                       <  1:776645/148‑1 (MQ=255)                gATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCg                                                                                                                                  <  1:382947/151‑1 (MQ=255)                 aTGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGG‑AAGCACCTTGAAAAGc                                                                                                                                   >  1:1121265/1‑149 (MQ=255)                  tGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑CAGCACATTTAAAAGCgg                                                                                                                                 >  1:321187/1‑150 (MQ=255)                   gCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGG‑GAGCACAATGaaaa                                                                                                                                     >  2:425765/1‑143 (MQ=255)                      tGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGAcc                                                                                                                              <  2:197034/149‑1 (MQ=255)                         cAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCa                                                                                                                          <  1:170375/150‑1 (MQ=255)                           gCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑CAGCACATTTAACAACGGGCCGGcca                                                                                                                         >  1:460276/1‑147 (MQ=255)                               ttGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑AAGCACATTGAACAACGGACCgg                                                                                                                            >  1:901513/1‑142 (MQ=255)                                        tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTcc                                                                                                          <  1:1014257/150‑1 (MQ=255)                                        tCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTcc                                                                                                          <  2:169007/150‑1 (MQ=255)                                         ccTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGt                                                                                                        <  2:499922/151‑1 (MQ=255)                                           tgtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGA‑AAGCACATTTTAAAGCGGAACGGGCACGC‑GGGGTCTCTcc                                                                                                          >  1:818930/1‑147 (MQ=255)                                            gtgACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGG‑CAGCACCTTTAAAAACGGGCCGGCAACCC‑GGGGTCTCTcccgt                                                                                                       >  1:131844/1‑146 (MQ=255)                                              gACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACgg                                                                                                   <  2:686583/151‑1 (MQ=255)                                                  tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCt                                                                                                 <  2:613535/149‑1 (MQ=255)                                                  tCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCt                                                                                                 <  1:825630/149‑1 (MQ=255)                                                     ggCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTAt                                                                                              <  2:815681/149‑1 (MQ=255)                                                        cTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGG‑AAGCACAATGAAAAGCGGAACGGCAACGC‑GGGGTCTCTCCCTTACGGccttttc                                                                                            >  1:857014/1‑142 (MQ=255)                                                        cTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGG‑CAGCACCATTAACAACGGGGCGGGAACGC‑GGGATCTCTCCCTTACgggctt                                                                                               >  1:777268/1‑141 (MQ=255)                                                          cAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGca                                                                                        <  1:837985/150‑1 (MQ=255)                                                            gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGG‑CAGCACATTTAAAAACGGACCGGCGACAC‑GGGTGCTTTCCGGTTCGGccttttc                                                                                            >  2:17039/1‑138 (MQ=255)                                                            gCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑GAGCAAACTTAACAACGGGGCCGCGACGC‑GGGGTTTCTCCCGTTCGGCt                                                                                                 >  2:804330/1‑139 (MQ=255)                                                             cTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑CAAAACAATGAACAACGGAACCGCCACCC‑GGGTTCTCTCCCCTACGGCTTAttc                                                                                            >  1:914829/1‑141 (MQ=255)                                                               tGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTa                                                                                   <  2:1124866/150‑1 (MQ=255)                                                                   cAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGAcc                                                                              <  1:1519/151‑1 (MQ=255)                                                                    aaGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑CAGCACATTTAAAAGCGGACCGGCCACGC‑GGGGTCTCTCCCTTTCGGCTTATTCTTCACCATACGAcc                                                                              >  1:411603/1‑150 (MQ=255)                                                                     agagATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGACCACAAAATTGAAAAG‑GGGCCGGCGAAGC‑GGGTTATCCCCGTTACGGCTTTTGCTTcacccttcta                                                                                >  2:324069/1‑140 (MQ=255)                                                                       agaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGttt                                                                          <  2:370879/151‑1 (MQ=255)                                                                        gaTGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGG‑CAGCACATTTTAAAACGGGACGGCCACGC‑GGGGTTTCTCCGTTACGGCCTTTTCTGCCCCTTACGACCgct                                                                           >  1:197786/1‑147 (MQ=255)                                                                                 gCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑AAGGACATTTAAAAACGGACCGGCAACGC‑GGGGTTTTTCCCGTACGGCCTCTTCTTCACATTACGACCCGTTGCGTaa                                                                    >  2:500137/1‑147 (MQ=255)                                                                                       agcagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGAt                                                           <  1:683901/150‑1 (MQ=255)                                                                                         cagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTAc                                                        <  1:489203/151‑1 (MQ=255)                                                                                         cagAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTAc                                                        >  2:1122498/1‑151 (MQ=255)                                                                                           gAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGt                                                      <  1:397640/151‑1 (MQ=255)                                                                                             aCAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTg                                                    <  2:684975/151‑1 (MQ=255)                                                                                              cAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGGGGGGA‑CAGCACATTGAAAAGCGGAGCCGCAACGC‑GGGGTCTCTCCCTTTCCGCCTATTCTGCACATTACGACCGCTTTCGTAACGAggggggtatgtt                                                     >  1:540656/1‑141 (MQ=255)                                                                                                gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGGGAGGGGG‑CAGCACAATGAACAACGGGCCCGGCACGC‑GGGGTCTCTCCGTTTCGGCGTTTTCTGCACATTACGACCGTTTGCGTAAAGAGGGGGGTACGtttggg                                                 >  2:527140/1‑147 (MQ=255)                                                                                                gCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGGG‑GAGCACAATGAACAACGGACCGGCAACGC‑GGGGTTTTTCCCTTTCGGGTTATTCTGCACCATACGACCGTTTGCGGTAAGAGGGGGTTTCGtt                                                     >  2:946454/1‑147 (MQ=255)                                                                                                    aaTGGTGTGCTGAAGGTCGTGAGCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑TGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTaa                                             <  2:644437/151‑1 (MQ=255)                                                                                                     aTGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑CAGCACATTTAACAGCGGGCCGGCAAAGC‑GGGGTCTCTCCCTTTCGGCTTATTCTGCACATTAAGACCCGTTTCGTAACAAGGGGGGTACGTTGGGGCTaa                                             >  1:163413/1‑150 (MQ=255)                                                                                                          gtgCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCg                                        <  2:818910/150‑1 (MQ=255)                                                                                                              tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑AAGGACATTTGAAAGCGGAACGGCCAAGC‑GGGTTCTCTCCGTTTCGGCTTATTCTGCCCAATACGACCGGTTTGGGAACGAGGGGGGTACGGTGGCGCTAAGccc                                         >  2:386191/1‑145 (MQ=255)                                                                                                              tGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGAGCGGGGA‑GAGCACAATTAAAACCGGAACGGGAACGCGGGATTTTTCCCGTTCGGGCTT‑TACTGCTCATTAAGACCGTTTGGGGGAAGAGGGGGGTTTGGTGGCGCCAAGccccgcg                                     >  1:1118309/1‑145 (MQ=255)                                                                                                                aaaGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGc                                 <  1:709762/151‑1 (MQ=255)                                                                                                                 aaGTCGTGATCAGTCGCGGTAGTGGCGGGCGGGGGGA‑CAGCACAATGAAAAACGGACCGGCAACGC‑GGGTTCTCTCCCTTTCGGCTTATCCTGCCCATTACGACCGTTTTCGTAACGAGGGGGTTACGTTTGCGCCAAAccccg                                       >  1:1150188/1‑142 (MQ=255)                                                                                                                     cGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTgggg                            <  1:166885/151‑1 (MQ=255)                                                                                                                       tGATCAGTCGCGGTAGTGGCGGGCGAGGGGA‑AAACACATTTAAAAGAGGGCCGGGAAAGC‑GGGGTTTCTCCGTTACGGCTTATTCTTCACATTACGACCCCTTGCGTAACGAGGGGGGTTCGTTTGGGCTAAGCCCggggggggggggggg                          >  1:1140969/1‑138 (MQ=255)                                                                                                                         aTCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGgcgc                         <  2:285612/150‑1 (MQ=255)                                                                                                                             gTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATc                     <  2:767701/150‑1 (MQ=255)                                                                                                                               cgcgGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCt                   <  2:369875/150‑1 (MQ=255)                                                                                                                                 cgGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCAt                <  2:762047/151‑1 (MQ=255)                                                                                                                                   gTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCt              <  2:463135/151‑1 (MQ=255)                                                                                                                                      gTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTg            <  2:769684/150‑1 (MQ=255)                                                                                                                                        ggcgggcgAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAg         <  2:321183/151‑1 (MQ=255)                                                                                                                                               gAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTaa  <  2:1020003/151‑1 (MQ=255)                                                                                                                                               gAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTa   >  1:974977/1‑150 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATCCAGCGACTACAGGATGCTTGTCAGCGGTTGATGATTTCCTGTGACTTCTGGCCTCAGCTTGAACAAGAGATGAAAACGCTGGCAGCAGAACAGCAAAATGGTGTGCTGAAAGTCGTGATCAGTCGCGGTAGTGGCGGGCGAGGGTA‑CAGCACATTGAACAGCGGACCGGCAACGC‑GGATTCTCTCCGTTACGGCTTATCCTGCACATTACGACCGTTTGCGTAACGAGGGGATTACGTTGGCGCTAAGCCCGGTGCGGCTGGGGCGCAATCCTCATCTTGCAGGTATTAA  >  NC_000913/1153426‑1153718

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