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.5%	155.5 / 22.3	55	coding (40/1389 nt)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (5/3);  new base A (22/25);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.69e-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	98.2%	174.6 / NA	55	L13F (TTG→TTT)	narU	nitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (26/28);  total (27/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543882‑1544170                                                                                                                                               |||                                                                                                                                                    agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑A‑aaaaagat                                                                                                                                              >  1:979943/1‑145 (MQ=255)  gaTTAACGATAACTGCGCTTAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTCCCTCGAGCAATATGTATGCCTTTATTCTCCCAGAAGGCCGGATTCTCAGATTTCGAGTC‑GCG‑A‑aaaaagata                                                                                                                                             >  2:855248/1‑144 (MQ=255)               tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTGTtc                                                                                                                                 <  2:10973/149‑1 (MQ=255)                         aGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCCGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:287616/149‑1 (MQ=255)                           cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGACCGGATTTTATTGTGTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                     <  2:167480/149‑1 (MQ=255)                            aTCCAGACACCGAGGGCAAGCAGTAGACAACTGACTGATAGCCAGAGACTTCTTCGAACACTGTGATTTCTTTTGTTTTCCCAGAAGGCCGGATTTTCTGTTTTCTTTTC‑GCGAC‑TCAGGATAACCACTATTTTTCTCGTTTTCCACTg                                                                                                                    >  2:886916/1‑149 (MQ=255)                             tCCACAAACAGAAGGCAAGCAGTAAACAACAGACTGTTGTCGGGAGGGTTCGTGGAGAAATATGGTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTCCAGTGc                                                                                                                   >  2:93420/1‑149 (MQ=255)                                    acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                            >  1:114579/1‑149 (MQ=255)                                     caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                           <  2:979943/149‑1 (MQ=255)                                        aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCGAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                        <  2:871015/149‑1 (MQ=255)                                            cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                    >  2:909723/1‑149 (MQ=255)                                                aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACGATGTTTCTCATTTTGCAGTGCCATATTGTTCGTCACATGc                                                                                                >  2:133418/1‑149 (MQ=255)                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTCTATTTGCCTTTTTGCATTGCCATATTGTTCCTCACATGCACACAt                                                                                          <  1:496777/149‑1 (MQ=255)                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacacct                                                                                          >  2:248730/1‑147 (MQ=255)                                                      aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTGTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                          <  1:187991/149‑1 (MQ=255)                                                        aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        >  1:520462/1‑149 (MQ=255)                                                           tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     >  1:965989/1‑149 (MQ=255)                                                            gactgaTATCCAGAGATTTCTTCGAGCAATGTGTTTTCCTGTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    <  1:588463/149‑1 (MQ=255)                                                               tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGa                                                                                 >  2:934208/1‑149 (MQ=255)                                                                   atCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaa                                                                             >  2:1004230/1‑149 (MQ=255)                                                                     ccAGAGAGATCTCCGATGAATAATTTGTTCTTTATTTTACCAGGAGCCCGGATTTTTTGGTTTCCAGTCGGCA‑A‑GAAAGATAACGACTGTTTTCCTCATTTTGCAGCGCCATATTGTTCCTCACATTCAAACATTGGTAATGaaaaaaa                                                                           <  2:28030/149‑1 (MQ=255)                                                                      cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           >  1:705902/1‑148 (MQ=255)                                                                      cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           <  2:705902/148‑1 (MQ=255)                                                                      cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGACAAAAAg                                                                          >  1:164364/1‑149 (MQ=255)                                                                            tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTGTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    <  1:839774/149‑1 (MQ=255)                                                                               cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                 >  1:31934/1‑149 (MQ=255)                                                                                ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                                >  2:715756/1‑149 (MQ=255)                                                                                  ggagaAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTGTCCAGTC‑GCGAA‑AAAAGATAACGACTAGTTCTCCCTTTTTGCGGTGCCATATTGTTCCGCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAg                                                              <  1:399085/148‑1 (MQ=255)                                                                                   gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAgg                                                             <  1:969558/149‑1 (MQ=255)                                                                                       aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCAAAGGCGCACTGGGTTAATGAAAAAAAGACAAAACACGAGGTAAg                                                         >  1:533227/1‑149 (MQ=255)                                                                                             ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGCCAAAACACGAGGTAAGgcg                                                      >  1:886870/1‑146 (MQ=255)                                                                                             ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgcg                                                      <  2:886870/146‑1 (MQ=255)                                                                                             ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCaa                                                   >  1:547352/1‑149 (MQ=255)                                                                                              tttCATTTACTTTCACAGAAGGCCGCATTTTCCGGTTTCCAGTT‑GCGAA‑AAAAGATAACGACTATGTTTCCCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAAt                                                  <  2:644098/149‑1 (MQ=255)                                                                                               tttcTTTATTTTCCCATAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                 <  2:141119/146‑1 (MQ=255)                                                                                               ttCGTTTATTGTCCCAGAAGTCCGTATTTTCTGGTTTCTAGTC‑GCGAA‑AAAAGCTAGCGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                 <  2:510212/149‑1 (MQ=255)                                                                                               ttCCTTCATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAACACAAAACACGAGGTAAGGCGCAATa                                                 <  2:578245/149‑1 (MQ=255)                                                                                                   tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             <  1:261203/149‑1 (MQ=255)                                                                                                   tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                             >  1:748730/1‑149 (MQ=255)                                                                                                   tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAGTGGTAATGAAAACAAGACAAAACACGAGGTGAGGAGCAATAACCa                                             >  1:171886/1‑149 (MQ=255)                                                                                                          tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGGCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGGCGAATAGACAGttatta                                      >  1:687531/1‑149 (MQ=255)                                                                                                            ccAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGa                                    <  2:965989/149‑1 (MQ=255)                                                                                                             cAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                   <  1:461016/149‑1 (MQ=255)                                                                                                                 agtccGGATTTCCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  1:636964/146‑1 (MQ=255)                                                                                                                   gCCGTATTTTCTGGTTTCCAGTA‑GCGAA‑AAAAGATAACGACTAATTTTCTCATTTTGCAGGGCCATATTGTTCCTCACATGCACACAATGGTAATAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAgg                             <  2:261829/149‑1 (MQ=255)                                                                                                                   gCCGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAACACCCAGGTAACACCCAATAGCCAGGTATTACAATTaaag                             >  2:53314/1‑147 (MQ=255)                                                                                                                   gCCGAATTTTCTGGTTTCCAGTC‑TCGAA‑AAAAGATAACGACTATTCTTCTCATTTTGCAGTGACATATTGTTCCTCACATGCACACAGTGGTAATGAACAAAAGACAAAATACGAGGTATGGCGCAATAGCCAGTTATTAGAATTAAgg                             <  1:855248/149‑1 (MQ=255)                                                                                                                     cGGATTTTCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                           >  1:378814/1‑149 (MQ=255)                                                                                                                          tttCTCGTTACCAGTT‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAtt                      <  1:385005/149‑1 (MQ=255)                                                                                                                           ttatgttttCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCAATATTGTTCCGCACATGCACACATAAGAAAAAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                     <  2:774772/143‑1 (MQ=255)                                                                                                                            tCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  2:748730/149‑1 (MQ=255)                                                                                                                            tCTGGTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:319664/149‑1 (MQ=255)                                                                                                                               ggTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTg                 >  2:592577/1‑149 (MQ=255)                                                                                                                                gTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAACAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGACTTAAGGATGAATTGGGTGa                >  1:165102/1‑149 (MQ=255)                                                                                                                                gTTTCCAGTC‑GCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGCGGTAAGGCGCAATAGCCAGCTATAAGAATTAAGGAAGGATTGGGGGa                >  2:87791/1‑149 (MQ=255)                                                                                                                                          tcgcgaaaaaaGATAACGAACATTTTTCGCATTTTGCAGTGCCATATTGTTCCTCACATGCAAACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGa        <  2:51034/142‑1 (MQ=255)                                                                                                                                           cgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt       <  2:114579/143‑1 (MQ=255)                                                                                                                                                aaaaaGATAACTACTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:886916/148‑1 (MQ=255)                                                                                                                                                aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:592577/148‑1 (MQ=255)                                                                                                                                               |||                                                                                                                                                    AGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTC‑GCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543882‑1544170

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