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.A91.0% 222.5 / 15.0 67coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (6/0);  new base A (33/28);  total (39/28)
Fisher's exact test for biased strand distribution p-value = 3.65e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 240.5 / NA 67L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (39/28);  total (39/28)

ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGA  >  NC_000913_3_pae_tpiA/1543878‑1544172
                                                                                                                                                 |||                                                                                                                                                      
aTTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑aaaaa                                                                                                                                                    >  1:185182/1‑149 (MQ=255)
     gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATaa                                                                                                                                               <  2:136083/150‑1 (MQ=255)
     gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATaa                                                                                                                                               <  2:769594/150‑1 (MQ=255)
     gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATaa                                                                                                                                               >  4:234365/1‑150 (MQ=255)
              tAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGa                                                                                                                                            <  2:724302/144‑1 (MQ=255)
              tAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGa                                                                                                                                            >  1:724302/1‑144 (MQ=255)
                 cTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTtc                                                                                                                                   <  1:564723/150‑1 (MQ=255)
                  tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAA‑GATAACGACTATTTTtct                                                                                                                                  >  2:235414/1‑150 (MQ=255)
                           tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCTGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCa                                                                                                                         >  4:480094/1‑150 (MQ=255)
                             gCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGt                                                                                                                       <  3:89818/150‑1 (MQ=255)
                              cATCCAGACACAGAGGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCAAGTAGAAAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                      <  2:61992/150‑1 (MQ=255)
                              cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                      >  1:557287/1‑150 (MQ=255)
                               aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                     <  1:752808/150‑1 (MQ=255)
                                tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCTTTTTGCAGTGcc                                                                                                                    >  4:424215/1‑150 (MQ=255)
                                 ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                                                   >  3:239123/1‑150 (MQ=255)
                                   aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                 >  3:55447/1‑150 (MQ=255)
                                      cacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                              <  2:549064/150‑1 (MQ=255)
                                       acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                             >  2:665179/1‑150 (MQ=255)
                                         aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                           >  1:113684/1‑150 (MQ=255)
                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                         <  3:73470/150‑1 (MQ=255)
                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                         <  1:244225/150‑1 (MQ=255)
                                             ggCAAGAAGTAGACTACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                       <  2:446420/150‑1 (MQ=255)
                                             ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                       >  2:354186/1‑150 (MQ=255)
                                               cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGCTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcaca                                                                                                     >  2:589033/1‑150 (MQ=255)
                                                 agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                   >  1:760745/1‑150 (MQ=255)
                                                 agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                   <  2:760745/150‑1 (MQ=255)
                                                 agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑A‑AAAACGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                   >  1:38237/1‑150 (MQ=255)
                                                   aagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGca                                                                                                 >  1:358892/1‑150 (MQ=255)
                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                               <  2:281707/150‑1 (MQ=255)
                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                            >  1:134600/1‑150 (MQ=255)
                                                           aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                         <  2:329371/150‑1 (MQ=255)
                                                           aaCTGACTGATATACAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAAAGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACCTTgg                                                                                         >  3:26726/1‑150 (MQ=255)
                                                                  tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                                  >  2:497939/1‑150 (MQ=255)
                                                                  tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                                  >  4:483530/1‑150 (MQ=255)
                                                                  tgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaa                                                                                  >  2:383850/1‑150 (MQ=255)
                                                                    atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                                >  2:510101/1‑150 (MQ=255)
                                                                    atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                                >  2:510139/1‑150 (MQ=255)
                                                                        ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                            >  1:577266/1‑150 (MQ=255)
                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAAc                                                                                                                                              <  3:197344/83‑1 (MQ=255)
                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAAc                                                                                                                                              >  4:197344/1‑83 (MQ=255)
                                                                         cAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGa                                                                           >  3:326491/1‑150 (MQ=255)
                                                                            agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                      <  2:487045/136‑1 (MQ=255)
                                                                            agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                      >  1:487045/1‑136 (MQ=255)
                                                                            agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                        <  1:657022/150‑1 (MQ=255)
                                                                              aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                      <  1:665179/150‑1 (MQ=255)
                                                                              aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                      >  2:572656/1‑150 (MQ=255)
                                                                               tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAac                                                                     <  1:235414/150‑1 (MQ=255)
                                                                               tttCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAac                                                                     >  1:323412/1‑150 (MQ=255)
                                                                                 tcttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAacac                                                                   <  3:492253/150‑1 (MQ=255)
                                                                                         cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                           <  1:354186/150‑1 (MQ=255)
                                                                                         cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                           <  1:572656/150‑1 (MQ=255)
                                                                                            tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑aa                                                                                                                                                       >  4:78535/1‑55 (MQ=255)
                                                                                            tatGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑a‑aa                                                                                                                                                       <  3:78535/55‑1 (MQ=255)
                                                                                                ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                     <  2:786924/117‑1 (MQ=255)
                                                                                                ttttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                     >  1:786924/1‑117 (MQ=255)
                                                                                                     cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAAAAAAAAAAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                               <  2:358892/150‑1 (MQ=255)
                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                              >  3:244107/1‑150 (MQ=255)
                                                                                                        tATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                                            >  4:258593/1‑150 (MQ=255)
                                                                                                                cAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAAt                                    >  2:427072/1‑150 (MQ=255)
                                                                                                                  gAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTa                                  >  3:473387/1‑150 (MQ=255)
                                                                                                                   aaGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                                 <  3:297891/150‑1 (MQ=255)
                                                                                                                      gCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGa                              <  1:427072/150‑1 (MQ=255)
                                                                                                                        cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATg                            >  1:789854/1‑150 (MQ=255)
                                                                                                                           aTTTTCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAAt                         >  2:334167/1‑150 (MQ=255)
                                                                                                                             tttCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTg                       >  3:117145/1‑150 (MQ=255)
                                                                                                                              ttCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                      <  4:108056/150‑1 (MQ=255)
                                                                                                                              ttCTGGTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                      <  2:260174/150‑1 (MQ=255)
                                                                                                                                   gTTTCCAGTCGCGAA‑AAAA‑GATAACGACTATTTTTCTTTTTTTTCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                 >  4:185551/1‑150 (MQ=255)
                                                                                                                                        cAGTCGCG‑AAAAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt            >  2:483624/1‑150 (MQ=255)
                                                                                                                                                 aaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa    <  1:449339/148‑1 (MQ=255)
                                                                                                                                                 aaaaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTCGGTGAAGTGCTGATTGaaa    >  2:821301/3‑150 (MQ=255)
                                                                                                                                                   aaaa‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGa  >  2:184410/1‑150 (MQ=255)
                                                                                                                                                 |||                                                                                                                                                      
ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAA‑GATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAGA  >  NC_000913_3_pae_tpiA/1543878‑1544172

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