Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 2,876,388 C→A intergenic (‑59/‑193) cysD ← / → iap sulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA2,876,3880CA100.0% 269.0 / NA 86intergenic (‑59/‑193)cysD/iapsulfate adenylyltransferase, subunit 2/aminopeptidase in alkaline phosphatase isozyme conversion
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (44/42);  total (44/42)

TTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTT  >  NC_000913_3_pae_tpiA/2876244‑2876517
                                                                                                                                                |                                                                                                                                 
tttGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTaa                                                                                                                               >  2:919156/1‑149 (MQ=255)
  tGAGAATTCTGCCGCCCCCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGGCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAaga                                                                                                                             <  2:1094337/149‑1 (MQ=255)
  tGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAaga                                                                                                                             <  1:945276/149‑1 (MQ=255)
   gagaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGGGAGTAcgag                                                                                                                            >  2:106905/1‑149 (MQ=255)
   gagaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAagag                                                                                                                            <  1:1161703/149‑1 (MQ=255)
    agaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAgaga                                                                                                                           >  1:1107304/1‑149 (MQ=255)
    agaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAgaga                                                                                                                           >  2:216637/1‑149 (MQ=255)
     taaTTCCGCCGACACCCCGCTTATCATGTGGATGCTTTCCGCATCCAGTTGCCGCATGTTTGTAAGTCGTATTTGATCCATAACCGTTCCATTGCAATACCGCAATTTTCTTGCCATCAGATGTTTCGAATATAGAGAGAGTAAGAGaa                                                                                                                          <  2:624358/148‑1 (MQ=255)
     gaATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGaa                                                                                                                          <  2:1062506/149‑1 (MQ=255)
     gaATGCCGCCGCTACCTCGCGAGTAATGGGGATGCTTTGCGTCTCCAGTTGCCGCAGGTGAGTAAGGCGTATTTGATCCATAACCGTTCCTTTGCAATGCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGaa                                                                                                                          <  2:149697/149‑1 (MQ=255)
      aaTTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAAc                                                                                                                         <  1:761360/149‑1 (MQ=255)
             ccgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGaa                                                                                                                  <  1:1082009/149‑1 (MQ=255)
              cgccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGaaa                                                                                                                 <  1:228881/149‑1 (MQ=255)
               gccACCTCGCGAATAATGTGGATGCTTGCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAAt                                                                                                                <  2:654205/149‑1 (MQ=255)
                ccACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAAtt                                                                                                               <  2:174220/149‑1 (MQ=255)
                    cTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCa                                                                                                           <  2:1436085/149‑1 (MQ=255)
                     tCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCTTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCaa                                                                                                          >  1:522102/1‑149 (MQ=255)
                        cgAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTa                                                                                                       <  1:986394/149‑1 (MQ=255)
                         gAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCATTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAAAGAACGAATGAAATTACCAATTAg                                                                                                      >  2:1016774/1‑149 (MQ=255)
                          aataatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGa                                                                                                     >  2:1331371/1‑149 (MQ=255)
                           ataatGTGGAGGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGaa                                                                                                    <  1:1283935/149‑1 (MQ=255)
                            taatGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAAt                                                                                                   <  1:216637/149‑1 (MQ=255)
                              atGTTGATGCTTTCCTCCTCCAGTTTCCGCAGGTGAGTTAGTCGTATTTGATGCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGa                                                                                                 <  1:996567/149‑1 (MQ=255)
                                 tgGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGGGAGTAAGAGAAGGAATGAAATTACTAATTAGAATGagta                                                                                              >  2:346447/1‑149 (MQ=255)
                                   gATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGagtagt                                                                                            >  1:389592/1‑149 (MQ=255)
                                   gATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGagtagt                                                                                            >  2:302582/1‑149 (MQ=255)
                                    aTGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTTGACTATAGGGAGAGGAAGAGAACGAATGAAATTACCAATTAGAAGGAGTGGtt                                                                                           >  1:360981/1‑149 (MQ=255)
                                          tCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTaa                                                                                     <  1:857649/149‑1 (MQ=255)
                                            cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACg                                                                                   <  1:514527/149‑1 (MQ=255)
                                            cGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTACCATCAGATGTTTCGACTATAGCGAAAGTAAGAGAATGAATGAAATTACCAATGAAAATGAGTAGTTCCACAACg                                                                                   >  2:763146/1‑149 (MQ=255)
                                             gCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACgg                                                                                  >  1:305540/1‑149 (MQ=255)
                                              ccTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATATTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGa                                                                                 >  2:1246199/1‑149 (MQ=255)
                                               cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGaa                                                                                >  1:510428/1‑149 (MQ=255)
                                               cTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGaa                                                                                >  1:194203/1‑149 (MQ=255)
                                               cTCCAGTTGCCGCAGGTGAGTAAGGCGTATTTGATCCATTTCCGTTCCTTTGGAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGaa                                                                                <  1:659524/149‑1 (MQ=255)
                                                tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAAc                                                                                    >  2:852532/1‑144 (MQ=255)
                                                tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAAc                                                                                    <  1:852532/144‑1 (MQ=255)
                                                tCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAAt                                                                               >  1:86123/1‑149 (MQ=255)
                                                tCCAGTTGCCGCAGGTGAGTAAGTCGTAAGTTCGCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTTCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAAt                                                                               >  2:708659/1‑149 (MQ=255)
                                                 ccAGTTTGCACATGTGAGTAAGTCGTATTTGTTCCATAACCGGTCCTCTGTAATACCGCTATTTTCTGGCCATCAAATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATa                                                                              <  2:45021/149‑1 (MQ=255)
                                                  cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTCGTTCCTTAACGGACTaa                                                                             >  1:1109297/1‑149 (MQ=255)
                                                  cAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATaa                                                                             <  2:389592/149‑1 (MQ=255)
                                                         cGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTg                                                                      >  1:1413382/1‑149 (MQ=255)
                                                             ggTGAGTGAGTGGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                  <  1:1081692/149‑1 (MQ=255)
                                                             ggTGAGGAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACAGCTAACTGATTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaa                                                                  <  1:1129715/149‑1 (MQ=255)
                                                              gTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCaaa                                                                 <  2:1413382/149‑1 (MQ=255)
                                                                gAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCCGATGTTTCGACTATAGGGAGAGTAAGAGAAAGACTGAAATTACCAATTAGAATGAAGAGTTCCTTACAGGTATAACGATTTGGCAAAGc                                                               >  2:424727/1‑149 (MQ=255)
                                                                   tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCATTTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                            >  2:1397588/1‑149 (MQ=255)
                                                                   tAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTaa                                                            >  1:1075513/1‑149 (MQ=255)
                                                                   tAAGTCGTATTTGATCCATAACCGTACCTTTGCAATACCGCTATTTTCTTGCCATCAGATTGCCCGACTAGATGTTCACTAGTACAAAGAATGAAATTAGCAATGAGAATGACTAGTTCCTTACAGGCATCACGATTTGCAAAAGCTaa                                                            >  2:364430/1‑149 (MQ=255)
                                                                    aaGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAat                                                           <  2:194203/149‑1 (MQ=255)
                                                                     aGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAata                                                          <  1:302582/149‑1 (MQ=255)
                                                                       tCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCa                                                       >  2:1379104/1‑149 (MQ=255)
                                                                                 gccATTACCGTTCCTTTGCAATACCTCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCt                                              <  1:685996/148‑1 (MQ=255)
                                                                                   cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCtt                                                                                       <  2:221951/106‑1 (MQ=255)
                                                                                   cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCtt                                                                                       >  1:221951/1‑106 (MQ=255)
                                                                                   cATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGACTTGGCAAAGCTAATATCAAAAAGTGCTTa                                            >  2:140492/1‑149 (MQ=255)
                                                                                     tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                          <  1:1331371/149‑1 (MQ=255)
                                                                                     tAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATGAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAg                                          <  2:1051971/149‑1 (MQ=255)
                                                                                      aaCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                         >  2:1151971/1‑149 (MQ=255)
                                                                                      aaCCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAgg                                         >  1:1275302/1‑149 (MQ=255)
                                                                                         cGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATACGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAACTGCTAAACGCAc                                      >  2:900764/1‑149 (MQ=255)
                                                                                         cGTTCCTTTGCAATACCGCTATTTTCTTGCCATCACATGTTTCGACTATAGGGAGAGTAAGAGAACGAATCAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCAc                                      <  1:1334069/149‑1 (MQ=255)
                                                                                          gTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATAAGTAGTTCAGTCACGAAATAACAATTTGGCAAGGCTAATATCAAAAAGTGCTTAAGGCAcc                                     >  2:1364649/1‑149 (MQ=255)
                                                                                               tttGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGAtt                                >  2:16583/1‑149 (MQ=255)
                                                                                                 tGCAATACCGCTCACTTCTTGCCATCAGATGTTTCGACTCTAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGAATTGGCAAACCTAATATCAAAAAGTGCTTAAGGCACCGGATTTc                              <  1:608188/149‑1 (MQ=255)
                                                                                                   cAATACCGATAATTTCTTTCCATCAAATCATTTGACTATAAGGAGATTAAGAGAACTAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGTCACCGGATTTCgg                            <  2:888817/149‑1 (MQ=255)
                                                                                                    aagaCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCACAGCTAAGATCAAAAAGTGCTTAAGGCACCGGATTTCggg                           <  1:856318/146‑1 (MQ=255)
                                                                                                    aaTACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCggg                           <  1:1364649/149‑1 (MQ=255)
                                                                                                       aCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGc                                               >  1:1163461/1‑126 (MQ=255)
                                                                                                       aCCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGc                                               <  2:1163461/126‑1 (MQ=255)
                                                                                                         cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTAGGCACTAGAAGGAGTTGGCTCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGttt                      >  1:1421774/1‑149 (MQ=255)
                                                                                                         cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGttt                      <  1:968855/149‑1 (MQ=255)
                                                                                                         cGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCGTAAGGGATCGGATGTCGGGCGttt                      >  2:759145/1‑149 (MQ=255)
                                                                                                          gCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAGTATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTa                     >  1:894564/1‑149 (MQ=255)
                                                                                                           cTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGAGCGTTTAg                    >  2:261270/1‑149 (MQ=255)
                                                                                                            tATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAgg                   >  1:218223/1‑149 (MQ=255)
                                                                                                              gtgtCTTGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCAAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGaa                 >  2:301156/4‑149 (MQ=255)
                                                                                                               ttacatGCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAg                <  1:611934/144‑1 (MQ=255)
                                                                                                                ttCTTGCCATTAAATTTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGa               <  1:228308/149‑1 (MQ=255)
                                                                                                                     gCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGaa                                                                                                                      >  2:1312138/1‑41 (MQ=255)
                                                                                                                     gCCATCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGaa                                                                                                                      <  1:1312138/41‑1 (MQ=255)
                                                                                                                        aTCAGATGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTGCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAAt       >  1:1253858/1‑149 (MQ=255)
                                                                                                                             aTGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGttt  <  1:1151971/149‑1 (MQ=255)
                                                                                                                             aTGTTTCGACTATAGGGAGAGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGttt  >  1:988368/1‑149 (MQ=255)
                                                                                                                             aTGTATCCACTATAGGGAGAGTATGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTGACGGACTAACGATTTGGCAAAGCTAATATCAAAAAGTGCTCAAGGCACCGAATTTCGGACGTTTAGGAGGATTTGAGATTGttt  >  2:1408104/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                 
TTTGAGAATTCTGCCGCCACCTCGCGAATAATGTGGATGCTTTCCGCCTCCAGTTGCCGCAGGTGAGTAAGTCGTATTTGATCCATAACCGTTCCTTTGCAATACCGCTATTTTCTTGCCATCAGATGTTTCGACTATAGGGAGCGTAAGAGAACGAATGAAATTACCAATTAGAATGAGTAGTTCCTTAACGGAATAACGATTTGGCAAAGCTAATATCAAAAAGTGCTTAAGGCACCGGATTTCGGGCGTTTAGGAAGATTTGAAATTGTTT  >  NC_000913_3_pae_tpiA/2876244‑2876517

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