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.2%	364.3 / 74.8	118	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand): ref base . (15/6); new base A (46/51); total (61/57)
Fisher's exact test for biased strand distribution p-value = 5.58e-02
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	100.0%	431.6 / NA	118	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand): ref base C (0/0); new base A (61/57); total (61/57)

Only 100 of 131 total aligned reads displayed.

TTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543876‑1544170

|||

ttATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑aa‑a                                                                                                                                                    <  1:386963/150‑1 (MQ=255)

 tATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aa                                                                                                                                                   >  2:410188/1‑150 (MQ=255)

   ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaag                                                                                                                                                 <  1:1026859/150‑3 (MQ=255)

   ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaag                                                                                                                                                 >  1:628995/1‑148 (MQ=255)

    tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCATTTATATTCCCAATAGTCCGGATTTTCTGATTTCCAGTCGAG‑A‑aa‑aaaga                                                                                                                                                <  2:331103/150‑4 (MQ=255)

      agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aaagata                                                                                                                                              >  4:620306/1‑145 (MQ=255)

        aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AA‑GATAAc                                                                                                                                            >  2:344718/1‑150 (MQ=255)

                               gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                     >  2:722737/1‑150 (MQ=255)

                                aaaccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTAAAAAAAAAAAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                    <  4:248953/147‑1 (MQ=255)

                                  tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  >  1:154077/1‑150 (MQ=255)

                                           aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTAtttt                                                                                                                                   <  1:1017347/124‑1 (MQ=255)

                                           aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTAtttt                                                                                                                                   >  2:1017347/1‑124 (MQ=255)

                                           aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  3:563661/150‑1 (MQ=255)

                                           aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  1:410188/150‑1 (MQ=255)

                                             aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               >  2:815563/1‑142 (MQ=255)

                                              aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      <  3:199474/150‑1 (MQ=255)

                                               ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  2:49032/150‑1 (MQ=255)

                                                  aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  >  1:1006605/1‑150 (MQ=255)

                                                  aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  >  1:34509/1‑150 (MQ=255)

                                                     aagTAGACAACTGACTGATATCGAGAGTTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTTGTTTCCAGTCGCG‑A‑AACAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                               >  4:163623/1‑150 (MQ=255)

                                                          gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATATTCTGGTTTCCAGACGCG‑AAAA‑AA‑AATAACGACTATTTTTAACATTTTGCATTGCTATATTATACCTTACATACACACAt                                                                                          >  2:581652/1‑150 (MQ=255)

                                                             aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATaa                                                                                                                                             <  3:552264/96‑1 (MQ=255)

                                                             aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATaa                                                                                                                                             >  4:552264/1‑96 (MQ=255)

                                                             aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTAttttt                                                                                                                                  >  2:1065149/1‑107 (MQ=255)

                                                               ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GAAAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   <  1:1160002/120‑1 (MQ=255)

                                                               ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AA‑GAAAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   >  2:1160002/1‑120 (MQ=255)

                                                                tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    >  1:441013/1‑150 (MQ=255)

                                                                tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AA‑AAGGATAACGCCTATTTTTATCATTTTGCAGTGCCATCTTGTTCCTCACATGCACACATTGGTaa                                                                                    >  3:669854/1‑150 (MQ=255)

                                                                      atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑aa‑a                                                                                                                                                    <  4:615104/80‑1 (MQ=255)

                                                                      atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑aa‑a                                                                                                                                                    >  3:615104/1‑80 (MQ=255)

                                                                        atCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                            <  2:628995/150‑1 (MQ=255)

                                                                          ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                 >  2:305285/1‑127 (MQ=255)

                                                                          ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                 <  1:305285/127‑1 (MQ=255)

                                                                          ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                          >  2:1222674/1‑150 (MQ=255)

                                                                           cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATaa                                                                                                                                             <  1:146830/82‑1 (MQ=255)

                                                                            agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       <  4:55109/135‑1 (MQ=255)

                                                                            agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                        >  1:1077323/1‑150 (MQ=255)

                                                                            agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aa‑aa                                                                                                                                                   >  3:266107/1‑75 (MQ=255)

                                                                              agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATaa                                                                                                                                             <  2:595311/79‑1 (MQ=255)

                                                                              agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATaa                                                                                                                                             >  1:595311/1‑79 (MQ=255)

                                                                              agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         >  1:1121434/1‑115 (MQ=255)

                                                                              agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  2:1121434/115‑1 (MQ=255)

                                                                                aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  1:1194876/104‑1 (MQ=255)

                                                                                     ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               >  2:238854/1‑150 (MQ=255)

                                                                                     ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               >  2:290657/1‑150 (MQ=255)

                                                                                          gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                          <  1:238854/150‑1 (MQ=255)

                                                                                             atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTg                                                                                                                         <  1:808915/84‑1 (MQ=255)

                                                                                             atatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTg                                                                                                                         >  2:808915/1‑84 (MQ=255)

                                                                                              tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   <  2:994440/89‑1 (MQ=255)

                                                                                              tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcg                                                      >  1:479194/1‑150 (MQ=255)

                                                                                              tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   >  1:994440/1‑89 (MQ=255)

                                                                                                 gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  <  1:1257579/103‑1 (MQ=255)

                                                                                                 gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  >  2:1257579/1‑103 (MQ=255)

                                                                                                     tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                       <  2:673509/126‑1 (MQ=255)

                                                                                                     tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                       >  1:673509/1‑126 (MQ=255)

                                                                                                        tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCGGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           >  2:2353/1‑119 (MQ=255)

                                                                                                        tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCGGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  1:2353/119‑1 (MQ=255)

                                                                                                            ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             >  2:1146898/1‑129 (MQ=255)

                                                                                                            ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                   <  3:564923/91‑1 (MQ=255)

                                                                                                             tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                       <  2:479194/150‑1 (MQ=255)

                                                                                                             tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatt                                       <  1:31046/150‑1 (MQ=255)

                                                                                                              ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                      >  2:226173/1‑150 (MQ=255)

                                                                                                              ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                      <  1:713518/150‑1 (MQ=255)

                                                                                                              ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                      <  1:344718/150‑1 (MQ=255)

                                                                                                                cccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                    <  4:24905/150‑1 (MQ=255)

                                                                                                                 ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               >  4:337465/1‑74 (MQ=255)

                                                                                                                 ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               <  3:337465/74‑1 (MQ=255)

                                                                                                                   aGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      >  4:240368/1‑113 (MQ=255)

                                                                                                                   aGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      <  3:240368/113‑1 (MQ=255)

                                                                                                                    gAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTa                                <  2:441013/150‑1 (MQ=255)

                                                                                                                         ccGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  >  1:855936/1‑63 (MQ=255)

                                                                                                                         ccGGATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  2:855936/63‑1 (MQ=255)

                                                                                                                         ccGGATTTTCTGGTTTCCAATGACG‑AAAA‑AA‑GATAACACCTATTTATCTCATTTTGCAGTGCAATATTGTGCCTCACATGCAAACAAAGGTAAAGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                           <  4:506781/150‑1 (MQ=255)

                                                                                                                           ggATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                           >  2:33295/1‑68 (MQ=255)

                                                                                                                            gATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                    >  1:1006629/1‑58 (MQ=255)

                                                                                                                            gATTTTCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                    <  2:1006629/58‑1 (MQ=255)

                                                                                                                                ttCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:736287/150‑1 (MQ=255)

                                                                                                                                ttCTGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    >  1:1245784/1‑150 (MQ=255)

                                                                                                                                   tGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTtct                                                                                                                                <  3:467261/39‑1 (MQ=255)

                                                                                                                                   tGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTTTTTTTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGTTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTAGGTg                 >  4:741270/1‑150 (MQ=255)

                                                                                                                                   tGGTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                                                 <  3:77689/54‑1 (MQ=255)

                                                                                                                                    ggTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      <  4:756965/64‑1 (MQ=255)

                                                                                                                                    ggTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                      >  3:756965/1‑64 (MQ=255)

                                                                                                                                     gTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa               >  4:124126/1‑150 (MQ=255)

                                                                                                                                     gTTTCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa               <  1:1129518/150‑1 (MQ=255)

                                                                                                                                      tttCCAGTCGCGAA‑AA‑AA‑TATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAg              <  4:669854/150‑1 (MQ=255)

                                                                                                                                      tttCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   <  2:1197931/49‑1 (MQ=255)

                                                                                                                                       ttCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                      >  2:271066/1‑77 (MQ=255)

                                                                                                                                       ttCCAGTCGCGAA‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                      <  1:271066/77‑1 (MQ=255)

                                                                                                                                          cAGTCGCG‑AAAA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  1:1148470/1‑150 (MQ=255)

                                                                                                                                          cAGTCGCG‑AAAA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  2:706714/1‑150 (MQ=255)

                                                                                                                                                gcgaaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  3:30426/105‑1 (MQ=255)

                                                                                                                                                gcgaaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  4:30426/6‑110 (MQ=255)

                                                                                                                                                gcgaaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTg     <  4:728286/145‑1 (MQ=255)

                                                                                                                                                 cgaaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGa    >  1:785691/5‑150 (MQ=255)

                                                                                                                                                   aaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  4:32603/3‑150 (MQ=255)

                                                                                                                                                   aaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:988392/148‑1 (MQ=255)

                                                                                                                                                   aaaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  3:666833/3‑150 (MQ=255)

                                                                                                                                                    aaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGTGCAATAGc                                               >  2:1203073/2‑104 (MQ=255)

                                                                                                                                                    aaa‑aa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAGGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGTGCAATAGc                                               <  1:1203073/103‑1 (MQ=255)

|||

TTATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AA‑AA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543876‑1544170

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