BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913_3_pae_tpiA	1,544,022	1 bp→AA	coding (39/1389 nt)	narU ←	nitrate/nitrite transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913_3_pae_tpiA	1,544,021	1	.	A	82.8%	299.9 / 60.6	99	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand): ref base . (11/6); new base A (45/37); total (56/43)
Fisher's exact test for biased strand distribution p-value = 5.93e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01
*	NC_000913_3_pae_tpiA	1,544,022	0	C	A	95.0%	329.8 / 6.1	99	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand): ref base C (3/2); new base A (53/41); total (56/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.05e-01

Only 100 of 110 total aligned reads displayed.

AACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543887‑1544170

|||

aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑a‑aaaagat                                                                                                                                              <  4:215439/144‑5 (MQ=255)

 aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACt                                                                                                                                        >  1:541998/1‑150 (MQ=255)

    gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGCGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTAtt                                                                                                                                     >  4:36660/1‑150 (MQ=255)

              tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCAttt                                                                                                                           >  4:689318/1‑150 (MQ=255)

              tAAATAGCATCCAGACACAGAAAGCAAGAAGAAGACAACTGACTGATAACCAGAAATTTCTACGAGCAATATGTTTTCCTTTATTTTCCCAGAAGACCGGATTTTCTGAATTCCAGTCGCAA‑A‑A‑AAAAATAACAACAATTTTTCTCAttt                                                                                                                           >  2:379120/1‑150 (MQ=255)

                 aTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCTAGTCGC‑G‑CAA‑AAAGATAAAAACTATTTTTTTCATTTTgt                                                                                                                        >  2:64083/1‑149 (MQ=255)

                 aTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGc                                                                                                                        <  1:169930/150‑1 (MQ=255)

                    gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTAttt                                                                                                                                    >  2:580197/1‑135 (MQ=255)

                    gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTAttt                                                                                                                                    <  1:580197/135‑1 (MQ=255)

                     cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                    <  1:843699/150‑1 (MQ=255)

                        ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTGCTCATTTTGCAGTGCCa                                                                                                                 >  3:562555/1‑150 (MQ=255)

                          aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               <  3:453169/150‑1 (MQ=255)

                              acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCGCAGGGGGCCGTGTTGTCTGTTTTCCAGTGTC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                           >  1:1343347/1‑150 (MQ=255)

                              acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                           >  3:919664/1‑150 (MQ=255)

                                    ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     >  4:66992/1‑150 (MQ=255)

                                     gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCATTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                    >  2:1095965/1‑150 (MQ=255)

                                        agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                 >  2:1118988/1‑150 (MQ=255)

                                          aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                               <  4:562555/150‑1 (MQ=255)

                                             tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTTTTTTTTTGCAGTGCCATATAGTTCCTCACATacacac                                                                                            >  4:580477/1‑150 (MQ=255)

                                               gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          >  2:441725/1‑150 (MQ=255)

                                               gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          <  4:65766/150‑1 (MQ=255)

                                               gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          >  3:65766/1‑150 (MQ=255)

                                                    ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     >  4:431391/1‑150 (MQ=255)

                                                    ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  3:259704/150‑1 (MQ=255)

                                                    ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTACAGAAGA‑AAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  4:627761/150‑1 (MQ=255)

                                                     tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    >  4:204604/1‑150 (MQ=255)

                                                         tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTAttttt                                                                                                                                  <  4:616652/100‑1 (MQ=255)

                                                           atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaa                                                                             <  3:801815/150‑1 (MQ=255)

                                                              tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           >  3:972094/1‑150 (MQ=255)

                                                              tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTACAGTAGA‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  4:710650/150‑1 (MQ=255)

                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAAc                                                                                                                                            <  3:968066/84‑1 (MQ=255)

                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       >  2:652541/1‑137 (MQ=255)

                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       <  1:652541/137‑1 (MQ=255)

                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAAc                                                                                                                                            >  4:968066/1‑84 (MQ=255)

                                                               ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          >  3:392720/1‑150 (MQ=255)

                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTtc                                                                                                                                 <  2:391674/94‑1 (MQ=255)

                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTtc                                                                                                                                 >  1:391674/1‑94 (MQ=255)

                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                         >  4:80808/1‑150 (MQ=255)

                                                                cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACACGCACACATTGGTAATGAAAAAAAGa                                                                         >  4:195540/1‑150 (MQ=255)

                                                                 agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAGGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                             <  2:341568/129‑1 (MQ=255)

                                                                 agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAGGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                             >  1:341568/1‑129 (MQ=255)

                                                                 agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTTTTTTTTTTGTTCGC‑GAA‑A‑AAAGATAACGACTATTTTCCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                        >  4:735310/1‑150 (MQ=255)

                                                                   agaTTTCTTCAAACAATATGTTTTCCTTTATTTTCCAAGAAGTCCGGATTTACTGCTACCTACTTTC‑GAA‑A‑AAAGAGAACGACTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      <  1:64083/150‑1 (MQ=255)

                                                                    gaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaa                                                                     >  1:616152/1‑150 (MQ=255)

                                                                     aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    <  1:561502/150‑1 (MQ=255)

                                                                      tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAac                                                                   >  1:1316567/1‑150 (MQ=255)

                                                                      tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAac                                                                   >  3:572446/1‑150 (MQ=255)

                                                                        tcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                 >  4:703320/1‑150 (MQ=255)

                                                                          ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTtctc                                                                                                                               <  2:924015/86‑1 (MQ=255)

                                                                          ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTtctc                                                                                                                               >  1:924015/1‑86 (MQ=255)

                                                                          ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               <  3:66992/150‑1 (MQ=255)

                                                                          ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               <  3:204604/150‑1 (MQ=255)

                                                                          ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               >  1:477714/1‑150 (MQ=255)

                                                                                     tGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCa                                                    <  2:477714/150‑1 (MQ=255)

                                                                                      gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                   >  3:171898/1‑150 (MQ=255)

                                                                                       ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                                  >  3:368508/1‑150 (MQ=255)

                                                                                            cTTTATTTTCCCAGAAGGCCGGCTTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGAGAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  3:850054/1‑150 (MQ=255)

                                                                                            cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACACGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  3:195540/150‑1 (MQ=255)

                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  3:463952/149‑1 (MQ=255)

                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  4:463952/1‑149 (MQ=255)

                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            >  4:704993/1‑150 (MQ=255)

                                                                                             tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑AGAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            >  2:71598/1‑150 (MQ=255)

                                                                                                 ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat                                        >  3:68294/1‑150 (MQ=255)

                                                                                                 ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat                                        >  3:69031/1‑150 (MQ=255)

                                                                                                   ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                      <  4:235680/150‑1 (MQ=255)

                                                                                                    tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAg                                     <  3:97154/150‑1 (MQ=255)

                                                                                                    tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAg                                     <  3:385023/150‑1 (MQ=255)

                                                                                                      ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                   >  3:806141/1‑150 (MQ=255)

                                                                                                       cAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                  >  1:769674/1‑150 (MQ=255)

                                                                                                       cAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                  >  4:682986/1‑150 (MQ=255)

                                                                                                       cAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTATCATTTTGCAGTGCCATATAGTACCTCACATGCAGACATTGGTAATGTAAAAAAGGCAAAACACGAGTTAAGGCGCAATAGCCAGTTATTAGAAt                                  >  3:338460/1‑150 (MQ=255)

                                                                                                          aaGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  4:590998/150‑1 (MQ=255)

                                                                                                          aaGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  2:229998/150‑1 (MQ=255)

                                                                                                               cGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  2:1414625/146‑1 (MQ=255)

                                                                                                               cGGATTTTCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              >  1:1414625/1‑146 (MQ=255)

                                                                                                               cGGATTTTCTGGTTTCCAGTCGC‑G‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGa                         >  1:285195/1‑150 (MQ=255)

                                                                                                                     ttCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:1118988/150‑1 (MQ=255)

                                                                                                                      tCTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  4:368508/150‑1 (MQ=255)

                                                                                                                       cTGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                  <  4:68294/150‑1 (MQ=255)

                                                                                                                       cTGGTTTCCAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                  <  4:69031/150‑1 (MQ=255)

                                                                                                                        tGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTg                 >  1:1431799/1‑150 (MQ=255)

                                                                                                                        tGGTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTg                 <  3:687382/150‑1 (MQ=255)

                                                                                                                          gTTTCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTAATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa               <  2:190109/150‑1 (MQ=255)

                                                                                                                           tttCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                            >  2:626053/1‑88 (MQ=255)

                                                                                                                           tttCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg              >  1:538875/1‑150 (MQ=255)

                                                                                                                           tttCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg              >  1:810755/1‑150 (MQ=255)

                                                                                                                             tCCAGTCGC‑GAA‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTg            <  1:441725/150‑1 (MQ=255)

                                                                                                                               cAGTCGC‑G‑AAA‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  1:1153085/1‑150 (MQ=255)

                                                                                                                                aGTCGC‑G‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                       <  1:1128585/103‑1 (MQ=255)

                                                                                                                                aGTCGC‑G‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                       >  2:1128585/1‑103 (MQ=255)

                                                                                                                                    tcgcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  4:572446/143‑1 (MQ=255)

                                                                                                                                     cgcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt      <  2:1316567/144‑1 (MQ=255)

                                                                                                                                     cgcgaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  1:1209169/143‑1 (MQ=255)

                                                                                                                                     cgcgaaa‑aaaGATAACGAATTTTTTTCTAATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAtt      <  4:338460/144‑1 (MQ=255)

                                                                                                                                       cgaaa‑aaaGATAACGACTATTTTTCTTATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGa    >  2:157963/5‑150 (MQ=255)

                                                                                                                                        gaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa   <  3:636514/147‑1 (MQ=255)

                                                                                                                                        gaaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa   >  2:1093200/4‑150 (MQ=255)

                                                                                                                                         aaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  3:36660/148‑1 (MQ=255)

                                                                                                                                         aaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  1:1231623/3‑150 (MQ=255)

                                                                                                                                         aaa‑aaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  2:769674/148‑1 (MQ=255)

|||

AACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑A‑AAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543887‑1544170

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 37 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`