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	85.2%	400.7 / 65.5	129	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand): ref base . (11/8); new base A (53/56); total (64/65)
Fisher's exact test for biased strand distribution p-value = 6.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913_3_pae_tpiA	1,544,022	0	C	A	99.2%	459.3 / NA	129	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand): ref base C (0/0); major base A (63/65); minor base T (1/0); total (64/65)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 144 total aligned reads displayed.

TGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGA  >  NC_000913_3_pae_tpiA/1543880‑1544172

|||

tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaaga                                                                                                                                                 <  1:1510324/150‑4 (MQ=255)

   gaTTAACGGTAACTGCGCTAAATAGCAGCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATaa                                                                                                                                               <  2:654088/150‑1 (MQ=255)

       aaCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGAc                                                                                                                                           >  2:1569259/1‑150 (MQ=255)

        aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAGCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AAGAAGAGAACGACt                                                                                                                                          >  2:414244/1‑150 (MQ=255)

         cGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCACAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGACAACGACTa                                                                                                                                         >  2:1461092/1‑150 (MQ=255)

              aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGCCTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTGTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTAttttt                                                                                                                                    <  3:77256/150‑1 (MQ=255)

              aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTAttttt                                                                                                                                    <  2:643245/150‑1 (MQ=255)

                 gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTtctc                                                                                                                                 <  1:721413/150‑1 (MQ=255)

                 gcgcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTtctc                                                                                                                                 >  2:1007091/1‑150 (MQ=255)

                      aaaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCAtttt                                                                                                                            <  3:370780/150‑1 (MQ=255)

                       aaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTg                                                                                                                           <  3:692405/150‑1 (MQ=255)

                         tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                         >  2:110631/1‑150 (MQ=255)

                         tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                         <  1:650691/150‑1 (MQ=255)

                         tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                         >  3:761420/1‑150 (MQ=255)

                          aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAg                                                                                                                        >  1:274473/1‑150 (MQ=255)

                              tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                    >  1:153079/1‑150 (MQ=255)

                                cAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                                  <  3:442729/150‑1 (MQ=255)

                                 aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                 <  3:484575/150‑1 (MQ=255)

                                   acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAtt                                                                                                               >  2:1518892/1‑150 (MQ=255)

                                   acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAtt                                                                                                               <  3:994826/150‑1 (MQ=255)

                                   acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAtt                                                                                                               >  3:70155/1‑150 (MQ=255)

                                   acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTACTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAtt                                                                                                               <  1:1647228/150‑1 (MQ=255)

                                     acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                             <  2:718457/150‑1 (MQ=255)

                                      aaGAAGGTAAGAAGAAGACTGGTGATTGATATCGAGCGATATCTTCTAGCAATCTGTTTTCCTTGATCATCCCAGAAGGCCGGATTTTCTGGTTTCCCGCCGCGTA‑AA‑AAGATAACGACAAACTTGCTCTGTTTGCAGCGTCCTgttgtt                                                                                                            <  4:215589/149‑1 (MQ=255)

                                       aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                         >  3:160766/1‑136 (MQ=255)

                                       aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                           <  1:1469500/150‑1 (MQ=255)

                                       aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                           >  2:1557987/1‑150 (MQ=255)

                                       aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                           >  1:1246402/1‑150 (MQ=255)

                                       aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaag                                                                                                                                                  >  2:36910/1‑108 (MQ=255)

                                            gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                      >  1:43248/1‑150 (MQ=255)

                                              aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                    <  2:1246402/150‑1 (MQ=255)

                                              aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                    >  2:596144/1‑150 (MQ=255)

                                                gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGc                                                                                                  <  4:878864/150‑1 (MQ=255)

                                                 aagTACACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                                 <  4:833715/150‑1 (MQ=255)

                                                  agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                                <  1:1361489/150‑1 (MQ=255)

                                                  agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                                <  2:1299451/150‑1 (MQ=255)

                                                  agTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcac                                                                                                <  4:527474/150‑1 (MQ=255)

                                                   gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAAc                                                                                                                                              <  4:308362/103‑1 (MQ=255)

                                                    tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTAtt                                                                                                                                       <  2:937114/109‑1 (MQ=255)

                                                      gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTGCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                            <  1:1007091/150‑1 (MQ=255)

                                                      gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                            <  1:86360/150‑1 (MQ=255)

                                                       aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                                           <  3:995632/150‑1 (MQ=255)

                                                         aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                         >  1:101568/1‑150 (MQ=255)

                                                         aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                         <  1:110631/150‑1 (MQ=255)

                                                           ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                 >  2:1634763/1‑124 (MQ=255)

                                                            tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                      >  1:1634843/1‑150 (MQ=255)

                                                            tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                      <  3:254800/150‑1 (MQ=255)

                                                                tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                                  >  1:26947/1‑150 (MQ=255)

                                                                tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                                  >  1:1228302/1‑150 (MQ=255)

                                                                  atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                                <  1:1557987/150‑1 (MQ=255)

                                                                  atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                                >  4:465305/1‑150 (MQ=255)

                                                                     tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                             <  3:584522/150‑1 (MQ=255)

                                                                     tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                             >  1:99869/1‑150 (MQ=255)

                                                                     tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATAAATAaaaaaaaaaa                                                                             <  2:117163/150‑1 (MQ=255)

                                                                      ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                            >  2:269873/1‑150 (MQ=255)

                                                                      ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                            >  3:85847/1‑150 (MQ=255)

                                                                       cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                      >  2:836845/1‑107 (MQ=255)

                                                                       cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                      <  1:836845/107‑1 (MQ=255)

                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCTAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                    >  2:74729/1‑108 (MQ=255)

                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                    <  1:74729/108‑1 (MQ=255)

                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacac                                                                                              >  2:82270/1‑130 (MQ=255)

                                                                        agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                          >  2:1333942/1‑150 (MQ=255)

                                                                              ttcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAaca                                                                    <  2:274473/150‑1 (MQ=255)

                                                                                 ttcGAGCAATATGTTTTCCTTTATTTTCCCAGTAGGCCGGATTTTCTGGTTTCCAGACGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCCGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 <  3:68537/150‑1 (MQ=255)

                                                                                 ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 <  1:1604027/150‑1 (MQ=255)

                                                                                 ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 <  1:1039763/150‑1 (MQ=255)

                                                                                 ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                 >  4:291445/1‑150 (MQ=255)

                                                                                  tcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg                                                                >  1:1169876/1‑150 (MQ=255)

                                                                                   cGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                               >  3:910452/1‑150 (MQ=255)

                                                                                       cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                           <  1:1333942/150‑1 (MQ=255)

                                                                                         atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                             <  3:419608/130‑1 (MQ=255)

                                                                                         atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                         <  3:157443/150‑1 (MQ=255)

                                                                                          tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcg                                                        <  2:675739/150‑1 (MQ=255)

                                                                                             gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                        <  1:503509/131‑1 (MQ=255)

                                                                                             gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                        >  2:503509/1‑131 (MQ=255)

                                                                                                 tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                                 <  4:85847/150‑1 (MQ=255)

                                                                                                 tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                                 <  4:910452/150‑1 (MQ=255)

                                                                                                   cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                               >  3:1009182/1‑150 (MQ=255)

                                                                                                    tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                              >  3:1016903/1‑150 (MQ=255)

                                                                                                    tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                              >  4:19674/1‑150 (MQ=255)

                                                                                                       aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTTATTTTTTAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCTGtta                                           >  4:470149/1‑150 (MQ=255)

                                                                                                              ccggaGGCCGGATTTTGTGGTTTCGAGTCGCGAA‑AA‑AAGATAACGACTCTTTTTGTCATTTTGCAGTGATATATTGTTCCTAATATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                    <  4:504969/146‑1 (MQ=255)

                                                                                                                  agcccGGAGTTTTTGGTTATCATTGGCG‑AAAA‑AAGATAACGACTATTTTTATAATTTGGCAGCGCCAGATTGTTCCTAACATGCACAAATTAGTAATGAAAAAAAGACAAAACACGGGGTAAGGCGCAATAGCCATTTATTTGAATTAAg                                <  2:980299/147‑1 (MQ=255)

                                                                                                                    gCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGTAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                              >  3:498035/1‑150 (MQ=255)

                                                                                                                     ccGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                             <  4:563158/150‑1 (MQ=255)

                                                                                                                     ccGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                             <  1:1303558/150‑1 (MQ=255)

                                                                                                                            ttCTGGTTTCCAGTCGCG‑AAAA‑AAGATAACGACTGTTTTTCTCATTTTGCAGAGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                      <  2:787235/150‑1 (MQ=255)

                                                                                                                             tCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                 >  4:536624/1‑58 (MQ=255)

                                                                                                                             tCTGGTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                 <  3:536624/58‑1 (MQ=255)

                                                                                                                                 gTTTCCAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                 >  2:155206/1‑150 (MQ=255)

                                                                                                                                     ccAGTCGCGAA‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc             >  1:202873/1‑150 (MQ=255)

                                                                                                                                      cAGTCGCG‑AAAA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt            >  3:56308/1‑150 (MQ=255)

                                                                                                                                          tcgcgaaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAATGCTGAt         <  2:1634843/143‑1 (MQ=255)

                                                                                                                                               aaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  3:573081/148‑1 (MQ=255)

                                                                                                                                               aaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  1:355579/148‑1 (MQ=255)

                                                                                                                                               aaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  1:1518892/148‑1 (MQ=255)

                                                                                                                                               aaaa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  3:971664/148‑1 (MQ=255)

                                                                                                                                                 aa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGa  <  3:291445/150‑1 (MQ=255)

                                                                                                                                                 aa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGa  <  4:603022/150‑1 (MQ=255)

                                                                                                                                                 aa‑aaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGa  >  1:367495/1‑150 (MQ=255)

|||

TGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AA‑AAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGA  >  NC_000913_3_pae_tpiA/1543880‑1544172

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