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.A89.7% 320.7 / 31.1 97coding (40/1389 nt)narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base . (9/1);  new base A (41/46);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 1.60e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01
*NC_000913_3_pae_tpiA1,544,0220CA100.0% 352.5 / NA 97L13F (TTG→TTT) narUnitrate/nitrite transporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (50/47);  total (50/47)

Only 100 of 103 total aligned reads displayed.

ATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543884‑1544170
                                                                                                                                           |||                                                                                                                                                   
aTTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAAc                                                                                                                                            >  4:478336/1‑150 (MQ=255)
          aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTAttttt                                                                                                                                  <  1:663181/150‑1 (MQ=255)
          aCTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTAttttt                                                                                                                                  <  1:51717/150‑1 (MQ=255)
           gtGCGCTAAATAGCATCCAGACAAAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTtc                                                                                                                                 <  2:582979/149‑1 (MQ=255)
            tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTtct                                                                                                                                >  2:872618/1‑150 (MQ=255)
            tGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTtct                                                                                                                                >  3:140483/1‑150 (MQ=255)
            tGCGCTAAATAGCATCCAGACACAGAACGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCGCATAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTtct                                                                                                                                >  1:800727/1‑150 (MQ=255)
                 tAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCAttt                                                                                                                           >  1:390705/1‑150 (MQ=255)
                  aaaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCAtttt                                                                                                                          <  1:985993/150‑1 (MQ=255)
                   aaTAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTg                                                                                                                         >  4:120727/1‑150 (MQ=255)
                     tAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  2:674144/150‑1 (MQ=255)
                        cATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTg                                                                                                                    >  2:42716/1‑150 (MQ=255)
                         aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   >  2:375143/1‑150 (MQ=255)
                         aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   >  2:845428/1‑150 (MQ=255)
                         aTCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGc                                                                                                                   >  1:460096/1‑150 (MQ=255)
                          tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  3:124887/150‑1 (MQ=255)
                          tCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGcc                                                                                                                  <  3:124907/150‑1 (MQ=255)
                            cAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCat                                                                                                                >  4:352110/1‑150 (MQ=255)
                             aGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCata                                                                                                               >  1:649564/1‑150 (MQ=255)
                               acacAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATAtt                                                                                                             <  1:663025/150‑1 (MQ=255)
                                cacaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTg                                                                                                            <  3:350515/150‑1 (MQ=255)
                                 acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  1:228280/138‑1 (MQ=255)
                                 acaGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGt                                                                                                           <  2:640565/150‑1 (MQ=255)
                                  caGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGtt                                                                                                          >  1:287127/1‑150 (MQ=255)
                                   agaagaCAAGTAGTAGACAACTGACTGATATCGAGAGATTTCTTCGAGGAATATTTTTTCCTTTAATTTCACAGAAGGCCGGATTTTCTGGTTTAAAATAAAAAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         <  2:739591/150‑1 (MQ=255)
                                   aGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTc                                                                                                         >  4:213349/1‑150 (MQ=255)
                                      aGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGGGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGTAGTGCCATATTGTTCCTc                                                                                                      >  2:894949/1‑150 (MQ=255)
                                       ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  1:231488/150‑1 (MQ=255)
                                       ggCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTGATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTca                                                                                                     <  1:520938/150‑1 (MQ=255)
                                        gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGGTTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCCTTTTGCAGTGCCATATTGTTCCTcac                                                                                                    >  4:384013/1‑150 (MQ=255)
                                        gCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTcac                                                                                                    <  2:152455/150‑1 (MQ=255)
                                         cAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTAttt                                                                                                                                    <  3:698896/117‑1 (MQ=255)
                                          aagaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACAt                                                                                                  >  3:319130/1‑150 (MQ=255)
                                           agaagTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATg                                                                                                 <  3:430095/150‑1 (MQ=255)
                                               gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                             <  2:782860/150‑1 (MQ=255)
                                               gTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcaca                                                                                             >  1:15755/1‑150 (MQ=255)
                                                tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacac                                                                                            >  1:702200/1‑150 (MQ=255)
                                                tAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGcacac                                                                                            >  2:87905/1‑150 (MQ=255)
                                                   aCAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACAtt                                                                                         >  2:1057483/1‑150 (MQ=255)
                                                    cAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGAGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTg                                                                                        <  2:316387/150‑1 (MQ=255)
                                                     aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       >  4:440458/1‑150 (MQ=255)
                                                     aaCTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       <  3:497051/150‑1 (MQ=255)
                                                       ctgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTa                                                                                     <  4:412992/150‑1 (MQ=255)
                                                        tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    >  4:522604/1‑150 (MQ=255)
                                                        tgactgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTaa                                                                                    <  3:470927/150‑1 (MQ=255)
                                                          actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATg                                                                                  <  2:707306/150‑1 (MQ=255)
                                                          actgaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATg                                                                                  >  1:635331/1‑150 (MQ=255)
                                                          actgaTATCCAGAGATTTCTTCGAGCAATATGTTCTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATg                                                                                  <  1:809093/150‑1 (MQ=255)
                                                             gaTATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaa                                                                               <  3:120727/150‑1 (MQ=255)
                                                                 tCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGaaaaaaa                                                                           >  4:456784/1‑150 (MQ=255)
                                                                    agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                        >  1:467797/1‑150 (MQ=255)
                                                                    agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                        >  1:907012/1‑150 (MQ=255)
                                                                    agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                        >  4:1509/1‑150 (MQ=255)
                                                                    agagATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGAc                                                                        >  2:170037/1‑150 (MQ=255)
                                                                      agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      >  3:36807/1‑150 (MQ=255)
                                                                      agaTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaa                                                                      >  3:382877/1‑150 (MQ=255)
                                                                        aTTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACaaaa                                                                    <  2:317074/150‑1 (MQ=255)
                                                                             ttcGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               <  4:140483/150‑1 (MQ=255)
                                                                                   cAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAg                                                         <  3:66896/150‑1 (MQ=255)
                                                                                    aaTATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAgg                                                        <  3:1509/150‑1 (MQ=255)
                                                                                           tttCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGAAAGGCGCAATa                                                 >  1:171010/1‑150 (MQ=255)
                                                                                             tCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                               <  2:635331/150‑1 (MQ=255)
                                                                                             tCCTTTATTATCCCAGAAGGCCGGATTTTCTGGTTTCCAGGCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTAGTAAAAAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGc                                               <  2:893135/150‑1 (MQ=255)
                                                                                             gcaatatttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGCCGCAATAGc                                               >  2:860840/8‑150 (MQ=255)
                                                                                              ccTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       <  3:428270/109‑1 (MQ=255)
                                                                                              ccTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       >  4:428270/1‑109 (MQ=255)
                                                                                                tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       <  3:427848/107‑1 (MQ=255)
                                                                                                tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTgg                                                                                       >  4:427848/1‑107 (MQ=255)
                                                                                                tttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            >  3:433354/1‑150 (MQ=255)
                                                                                                tttATTTTCCCAGAAGGCCAGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  3:547806/150‑1 (MQ=255)
                                                                                                 ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                           <  1:157608/150‑1 (MQ=255)
                                                                                                 ttATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGt                                           >  3:625470/1‑150 (MQ=255)
                                                                                                  tATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGtt                                          <  3:586253/150‑1 (MQ=255)
                                                                                                    ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            <  2:341550/146‑1 (MQ=255)
                                                                                                    ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAg                                            >  1:341550/1‑146 (MQ=255)
                                                                                                    ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat                                        >  4:558642/1‑150 (MQ=255)
                                                                                                    ttttCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttat                                        >  4:408198/1‑150 (MQ=255)
                                                                                                         ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                       >  2:274768/1‑130 (MQ=255)
                                                                                                         ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGgc                                                       <  1:274768/130‑1 (MQ=255)
                                                                                                         ccAGAAGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGttaatgaag                                   >  4:682952/1‑144 (MQ=255)
                                                                                                             aaGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTaa                               <  4:433354/150‑1 (MQ=255)
                                                                                                              aGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  3:352110/150‑1 (MQ=255)
                                                                                                              aGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  2:467797/150‑1 (MQ=255)
                                                                                                              aGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  2:649564/150‑1 (MQ=255)
                                                                                                              aGGCCGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAg                              <  1:845428/150‑1 (MQ=255)
                                                                                                                  cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               <  2:518978/113‑1 (MQ=255)
                                                                                                                  cGGATTTTCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGa                                                               >  1:518978/1‑113 (MQ=255)
                                                                                                                        ttCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    <  1:42716/150‑1 (MQ=255)
                                                                                                                        ttCTGGTTTCCAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTgg                    >  3:210479/1‑150 (MQ=255)
                                                                                                                         tCTGGTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTggg                   <  2:391636/150‑1 (MQ=255)
                                                                                                                           tGGTATCCAGTCGCGAA‑AAAAAATAACGACTATTTTTCTCATTTTGAACGGAAATATGGTTCCTCACAAACACACGTTGGTAATGAAAAAGATACAAGACACGGCGTAAGGCGCAATAACCAGTTATTAGAATTAAGGATGAATTGGGTg                 >  2:707441/1‑150 (MQ=255)
                                                                                                                             gTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       >  1:938198/1‑46 (MQ=255)
                                                                                                                             gTTTCCAGTCGCGAA‑AAAAGATAACGACTATTTTTCTCATTTTGCa                                                                                                                       <  2:938198/46‑1 (MQ=255)
                                                                                                                                  cAGTCGCG‑AAAAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCt          >  1:970963/1‑150 (MQ=255)
                                                                                                                                      tcgcgaaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGGTGAATTGGGTGAAGTGCTGAt       >  1:443445/8‑150 (MQ=255)
                                                                                                                                           aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:872618/148‑1 (MQ=255)
                                                                                                                                           aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  1:924600/148‑1 (MQ=255)
                                                                                                                                           aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  <  2:907012/148‑1 (MQ=255)
                                                                                                                                           aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  3:139349/3‑150 (MQ=255)
                                                                                                                                           aaaaaaGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGaaa  >  1:1046157/3‑150 (MQ=255)
                                                                                                                                           |||                                                                                                                                                   
ATTAACGGTAACTGCGCTAAATAGCATCCAGACACAGAAGGCAAGAAGTAGACAACTGACTGATATCCAGAGATTTCTTCGAGCAATATGTTTTCCTTTATTTTCCCAGAAGGCCGGATTTTCTGGTTTCCAGTCGCG‑C‑AAAAGATAACGACTATTTTTCTCATTTTGCAGTGCCATATTGTTCCTCACATGCACACATTGGTAATGAAAAAAAGACAAAACACGAGGTAAGGCGCAATAGCCAGTTATTAGAATTAAGGATGAATTGGGTGAAGTGCTGATTGAAA  >  NC_000913_3_pae_tpiA/1543884‑1544170

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