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_tpiA1,544,0211.A83.7% 253.0 / 47.2 86coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (9/5);  new base A (35/37);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 3.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_tpiA1,544,0220CA98.8% 295.2 / NA 86L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base A (43/42);  minor base T (1/0);  total (44/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_pae_tpiA/1543877‑1544171
                                                                                                                                                  |||                                                                                                                                                     
tATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACGGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAGGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑a                                                                                                                                                   <  4:416654/150‑1 (MQ=255)
 aTTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aa                                                                                                                                                  <  3:221298/150‑2 (MQ=255)
  ttGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aag                                                                                                                                                 >  3:46768/1‑148 (MQ=255)
         tAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCTGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGa                                                                                                                                           >  2:427622/1‑150 (MQ=255)
         tAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGa                                                                                                                                           <  1:216351/150‑1 (MQ=255)
           aCGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACt                                                                                                                                         >  1:430714/1‑150 (MQ=255)
                   tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtct                                                                                                                                 >  1:443188/1‑150 (MQ=255)
                       cTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCAtt                                                                                                                             >  3:318382/1‑150 (MQ=255)
                         aaaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCAtttt                                                                                                                           <  2:59576/150‑1 (MQ=255)
                                aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                    <  4:310183/150‑1 (MQ=255)
                                aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aa                                                                                                                                                  >  4:444159/1‑118 (MQ=255)
                                aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aa                                                                                                                                                  <  3:444159/119‑2 (MQ=255)
                                    aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                >  2:335041/1‑150 (MQ=255)
                                          aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCTTTTTGCAGTGCCATATTGTTc                                                                                                          >  3:401094/1‑150 (MQ=255)
                                          aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          >  1:606475/1‑150 (MQ=255)
                                           gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCCTTATTTACCCAGAAGGCCGGATTTTCTGGTGAAAAAAAAAAAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                         <  4:200045/150‑1 (MQ=255)
                                            aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                        >  1:484466/1‑150 (MQ=255)
                                            aaGGCAAGAAGTAGACAACAACATGAGGACGTGAGAACTCGCCGAGCGATAGGTTTGCCTTGATATTCCCAGAAGGCCGGATTTTCTAGTTTCCAGACGAGAA‑AAA‑AAAAAACGACTATTTTACGCAGGTTGCATTGCCAGATGATTCCt                                                                                                        >  4:117113/1‑150 (MQ=255)
                                                      gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                              <  4:246383/150‑1 (MQ=255)
                                                       tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aaga                                                                                                                                                <  4:93683/98‑4 (MQ=255)
                                                       tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaa‑aaga                                                                                                                                                >  3:93683/1‑95 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTTTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATGTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  1:718178/1‑150 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           <  2:85339/150‑1 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  4:46808/1‑150 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           <  4:46768/150‑1 (MQ=255)
                                                         gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  2:316464/1‑150 (MQ=255)
                                                             actgacttATATCAAAAGATTACTTCGATAAATATTTTTTCCTTTATTTTATTAGAAGGCCGGATTTTCTGGTTAAAAATAAAAAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGt                                                                                       <  2:545724/150‑1 (MQ=255)
                                                                gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                    >  2:295110/1‑150 (MQ=255)
                                                                   tgaTATCCAGAGATTTCTTCGAGCCATATGTTTTCCTTTATTTTCCCAGAATGCCGGATTTTCTGGTTTCCCGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCTATAGTGTTCCTCACATTCACAAATTGGTAATGaa                                                                                 >  2:475566/1‑150 (MQ=255)
                                                                   tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAAc                                                                                                                                             >  3:365205/1‑90 (MQ=255)
                                                                   tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAAc                                                                                                                                             <  4:365205/90‑1 (MQ=255)
                                                                     atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                               >  1:502921/1‑150 (MQ=255)
                                                                        tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                            <  4:278704/150‑1 (MQ=255)
                                                                         ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCTTAGTTGGTTTTATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTTCAGTGCCATATTTTTCCTCACATGCACACATTGGTAATTTAAATAAg                                                                           >  2:645415/1‑150 (MQ=255)
                                                                         ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAATGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaaa                                                                           >  2:645333/1‑149 (MQ=255)
                                                                         ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑T‑AAACAGCTAACGTCTATTTGTCTCATTTTGCGGTGCCATATTGGTCCTCACATGCACTCATTTGTAATGAAAAAAAg                                                                           >  1:641655/1‑150 (MQ=255)
                                                                          cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          <  2:295959/118‑1 (MQ=255)
                                                                          cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          >  1:295959/1‑118 (MQ=255)
                                                                          cAGAGATTTCTTCGAGCAATATGTTTCCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          <  3:117113/150‑1 (MQ=255)
                                                                             agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                       <  1:584894/150‑1 (MQ=255)
                                                                             agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                       <  1:585027/150‑1 (MQ=255)
                                                                                   tttcGAGTGATATGTTTGCCTTTATTTTCCCAGAAGGCCGGATTTTATGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCCCATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                                 <  2:641655/149‑1 (MQ=255)
                                                                                   cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                                 <  2:430714/150‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        <  4:105994/94‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        >  3:105994/1‑94 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                <  2:502921/150‑1 (MQ=255)
                                                                                    ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                <  3:198863/150‑1 (MQ=255)
                                                                                       gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt                                                             <  2:484466/150‑1 (MQ=255)
                                                                                        aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacac                                                                                             <  1:306473/117‑1 (MQ=255)
                                                                                        aGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacac                                                                                             >  2:306473/1‑117 (MQ=255)
                                                                                         gacatatGTGCTACTCTATGTTCCCAGGAGGCTGGTTTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAAGAAAAAAAGACAAAACACGAGGTaa                                                           <  2:71586/147‑1 (MQ=255)
                                                                                              atGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTTCATTACCATATTATTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGCGGTAAGgcgc                                                      >  2:720098/1‑150 (MQ=255)
                                                                                                gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                    >  1:232674/1‑86 (MQ=255)
                                                                                                gTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                    <  2:232674/86‑1 (MQ=255)
                                                                                                     aataaaatttCACAGAAGGCCGGATTTTCTGGTTTCAAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCAATATTGTTCCTCACATGCACACAAAAAAAAAAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                               <  4:284780/143‑1 (MQ=255)
                                                                                                      cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          >  1:119976/1‑122 (MQ=255)
                                                                                                      cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          <  2:119976/122‑1 (MQ=255)
                                                                                                      cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  1:60881/1‑150 (MQ=255)
                                                                                                      cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTAGTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  1:441644/1‑150 (MQ=255)
                                                                                                       tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  1:161327/1‑104 (MQ=255)
                                                                                                       tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           <  2:161327/104‑1 (MQ=255)
                                                                                                       tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                             >  4:240822/1‑150 (MQ=255)
                                                                                                         tATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                                           >  1:435485/1‑150 (MQ=255)
                                                                                                           ttttCCCAGAAGGCCGGATTTTCTGGGTTCCAGTCGCG‑AAAAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAAAAAACAAGAGGTAAGGCGCAATAGCCAGttat                                         >  3:450954/1‑150 (MQ=255)
                                                                                                            tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtc                                                                                                                                  >  3:336897/1‑60 (MQ=255)
                                                                                                            tttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTtc                                                                                                                                  <  4:336897/60‑1 (MQ=255)
                                                                                                             ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                       <  3:46808/150‑1 (MQ=255)
                                                                                                              tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAg                                      <  2:443188/150‑1 (MQ=255)
                                                                                                              tCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACAAGAGGTAAGGCGCAAAAGCCAGTTATTAg                                      >  1:233893/1‑150 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCCGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          <  3:92196/112‑1 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                          >  4:92196/1‑112 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                               >  2:371228/1‑139 (MQ=255)
                                                                                                                ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGcc                                               <  1:371228/139‑1 (MQ=255)
                                                                                                                              tttCTGGTTTCCAGCAGCGAA‑AAA‑AGATAACGCCTATTTTCCTCATTTTGCAGTGCCCTATTGTTCCTCACATGCAGCCATTGGTCATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                      <  3:467370/150‑1 (MQ=255)
                                                                                                                               atCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAAAACGAGGTAAGGCGCAAAAGACAGTTATTAGAATTAAGGATGAATTgg                     >  2:26155/2‑150 (MQ=255)
                                                                                                                                tCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                    >  1:368941/1‑150 (MQ=255)
                                                                                                                                tCTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                    >  1:368959/1‑150 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                     >  2:515014/1‑68 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                     >  2:515113/1‑68 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGt                   <  1:26155/150‑1 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATGGTTCCTcac                                                                                                     <  1:515014/68‑1 (MQ=255)
                                                                                                                                 cTGGTTTCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATGGTTCCTcac                                                                                                     <  1:515113/68‑1 (MQ=255)
                                                                                                                                     tttCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                            <  3:148688/121‑1 (MQ=255)
                                                                                                                                     tttCCAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                            >  4:148688/1‑121 (MQ=255)
                                                                                                                                        ccAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc            <  2:368959/150‑1 (MQ=255)
                                                                                                                                        ccAGTCGCGAA‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc            <  2:368941/150‑1 (MQ=255)
                                                                                                                                                   aaaa‑aGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  >  1:561238/2‑150 (MQ=255)
                                                                                                                                                  |||                                                                                                                                                     
TATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAA‑AGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_pae_tpiA/1543877‑1544171

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