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.A88.3% 252.6 / 27.6 77coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (4/5);  new base A (35/33);  total (39/38)
Fisher's exact test for biased strand distribution p-value = 7.37e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 279.7 / NA 77L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (39/38);  total (39/38)

ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_pae_tpiA/1543878‑1544171
                                                                                                                                                  |||                                                                                                                                                    
aTTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑aaaaa                                                                                                                                                  <  3:56057/150‑2 (MQ=255)
  tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑aaaaaga                                                                                                                                                <  2:199259/150‑4 (MQ=255)
  tGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑aaaaaga                                                                                                                                                >  3:74505/1‑147 (MQ=255)
    agaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑A‑aaaaagata                                                                                                                                              <  1:18590/150‑6 (MQ=255)
     gaTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATaa                                                                                                                                              >  3:152647/1‑150 (MQ=255)
      aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAAc                                                                                                                                             >  2:351736/1‑150 (MQ=255)
             gTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAtt                                                                                                                                      >  1:912174/1‑150 (MQ=255)
                  tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTtct                                                                                                                                 >  4:437882/1‑150 (MQ=255)
                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCAt                                                                                                                              >  4:86575/1‑150 (MQ=255)
                     gcTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCAt                                                                                                                              >  3:68237/1‑150 (MQ=255)
                      cTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAttttt                                                                                                                                   <  4:381318/144‑1 (MQ=255)
                      cTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAttttt                                                                                                                                   >  3:381318/1‑144 (MQ=255)
                       tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCAttt                                                                                                                            <  1:487162/150‑1 (MQ=255)
                           tAGCATCCAGACACCGTGGCCGAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        <  1:271195/150‑1 (MQ=255)
                           tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        <  2:160616/150‑1 (MQ=255)
                           tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        <  2:912174/150‑1 (MQ=255)
                           tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        >  1:942991/1‑150 (MQ=255)
                                 ccAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCa                                                                                                                  <  1:149666/150‑1 (MQ=255)
                                   aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                                >  3:32578/1‑150 (MQ=255)
                                     acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAtt                                                                                                              >  1:188576/1‑150 (MQ=255)
                                         aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                          >  1:219387/1‑150 (MQ=255)
                                          gAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                         >  2:387464/1‑150 (MQ=255)
                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                        <  3:86575/150‑1 (MQ=255)
                                           aaGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCt                                                                                                        <  4:39108/150‑1 (MQ=255)
                                            aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTc                                                                                                       <  2:188576/150‑1 (MQ=255)
                                              gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATCTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                     <  1:632757/150‑1 (MQ=255)
                                                aagacgTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCTATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                   >  2:34675/1‑150 (MQ=255)
                                                  gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        >  3:262444/1‑127 (MQ=255)
                                                  gaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                        <  4:262444/127‑1 (MQ=255)
                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAt                                                                                                                                       >  4:593557/1‑109 (MQ=255)
                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                              >  3:117679/1‑150 (MQ=255)
                                                     gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                              >  1:186234/1‑150 (MQ=255)
                                                     gTAGACAACTGACTGATATCCAGAGATTTCTGCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTAt                                                                                                                                       <  3:593557/109‑1 (MQ=255)
                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  2:2212/1‑150 (MQ=255)
                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  1:1029353/1‑150 (MQ=255)
                                                        gACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAt                                                                                           >  1:101405/1‑150 (MQ=255)
                                                            actgactgATATTCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTTTTCCTCACATGCACACATTGGt                                                                                       <  4:524419/150‑1 (MQ=255)
                                                               gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                    <  2:301789/150‑1 (MQ=255)
                                                               gactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAAt                                                                                    <  1:330087/150‑1 (MQ=255)
                                                                    atatCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaa                                                                               >  2:470850/1‑150 (MQ=255)
                                                                      atCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaa                                                                             >  1:342194/1‑150 (MQ=255)
                                                                        ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                           <  1:351736/150‑1 (MQ=255)
                                                                        ccAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAg                                                                           >  3:356718/1‑150 (MQ=255)
                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                         >  3:85376/1‑150 (MQ=255)
                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                         >  1:502733/1‑150 (MQ=255)
                                                                          agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                         >  1:638406/1‑150 (MQ=255)
                                                                            agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                       >  2:194172/1‑150 (MQ=255)
                                                                              aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                     <  2:143350/150‑1 (MQ=255)
                                                                                  cttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACg                                                                 <  1:40834/150‑1 (MQ=255)
                                                                                   ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                                <  2:608665/150‑1 (MQ=255)
                                                                                      gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt                                                             >  1:257277/1‑150 (MQ=255)
                                                                                      gAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGt                                                             >  4:376286/1‑150 (MQ=255)
                                                                                        gCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTaa                                                           >  3:455126/1‑150 (MQ=255)
                                                                                                 tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATa                                                  <  4:318957/150‑1 (MQ=255)
                                                                                                   tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAAAAAAAAAGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                                <  2:942991/150‑1 (MQ=255)
                                                                                                     cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              <  1:387464/150‑1 (MQ=255)
                                                                                                     cTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCa                                              >  4:236326/1‑150 (MQ=255)
                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                             >  2:898806/1‑150 (MQ=255)
                                                                                                      tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                             >  4:448353/1‑150 (MQ=255)
                                                                                                        tATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                                           <  1:34675/150‑1 (MQ=255)
                                                                                                         tttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtta                                          <  1:386655/149‑1 (MQ=255)
                                                                                                         aTTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtta                                          <  2:1025147/150‑1 (MQ=255)
                                                                                                            ttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttatta                                       <  2:219387/150‑1 (MQ=255)
                                                                                                             tCCCATAAGGCCGTATTTTATGGTTACCAGTTGCAA‑A‑AAAAAATACCGCCTATTTTTCACAATTTTCAGTGTTATATTGTTCCTATAATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTTAGGTGCAATAGCCAGTTATTAg                                      <  4:286838/150‑1 (MQ=255)
                                                                                                               ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGaa                                    >  1:893638/1‑150 (MQ=255)
                                                                                                                cAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                        <  2:108859/113‑1 (MQ=255)
                                                                                                                cAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACa                                                                        >  1:108859/1‑113 (MQ=255)
                                                                                                                       ccGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                            <  2:502733/150‑1 (MQ=255)
                                                                                                                       ccGGATTTTCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGAt                            >  4:225128/1‑150 (MQ=255)
                                                                                                                              ttCTGGTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                     <  2:1029353/150‑1 (MQ=255)
                                                                                                                                   gTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCTTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                >  2:416833/1‑150 (MQ=255)
                                                                                                                                   gTTTCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGaa                <  1:734796/150‑1 (MQ=255)
                                                                                                                                    tttCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                         <  1:119413/60‑1 (MQ=255)
                                                                                                                                    tttCCAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTcc                                                                                                         >  2:119413/1‑60 (MQ=255)
                                                                                                                                       ccAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc            <  1:275158/150‑1 (MQ=255)
                                                                                                                                       ccAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGc            <  1:194172/150‑1 (MQ=255)
                                                                                                                                        cAGTCGC‑GAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt           <  2:342194/150‑1 (MQ=255)
                                                                                                                                                 gaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt        <  3:504329/143‑1 (MQ=255)
                                                                                                                                                 gaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGAt        >  4:504329/4‑146 (MQ=255)
                                                                                                                                                 gaaaaaaGATAACGACTAACTTTCTCATTTTGCAGTGCCATATTGTTCTTCACATGAACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaa    <  2:409876/147‑1 (MQ=255)
                                                                                                                                                   aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  >  2:345988/2‑150 (MQ=255)
                                                                                                                                                   aaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAg  >  3:524433/2‑150 (MQ=255)
                                                                                                                                                  |||                                                                                                                                                    
ATTGAGATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGC‑G‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAAG  >  NC_000913_3_pae_tpiA/1543878‑1544171

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: ‑