Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913_3_pae_tpiA 4,545,683 C→A 7.8% D224E (GAC→GAA)  fimD → fimbrial usher outer membrane porin protein; FimCD chaperone‑usher

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,545,6830CA 7.8% 159.7 / 4.4 63D224E (GAC→GAA) fimDfimbrial usher outer membrane porin protein; FimCD chaperone‑usher
Reads supporting (aligned to +/- strand):  ref base C (31/27);  new base A (3/2);  total (34/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.23e-01

TGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTC  >  NC_000913_3_pae_tpiA/4545536‑4545830
                                                                                                                                                   |                                                                                                                                                   
tGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACa                                                                                                                                                    >  1:436121/1‑149 (MQ=255)
    ttATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATac                                                                                                                                                >  2:448606/1‑149 (MQ=255)
            tCCCGGTATAAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGa                                                                                                                                        <  1:448443/149‑1 (MQ=255)
               cGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTGAACCTACTAACTTCGATAAATATTGGTGCGTGGCGTTTACGCGACACTACCACCTTGAGtt                                                                                                                                     >  1:180658/1‑149 (MQ=255)
               cGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTCCGTGTAGTTTACGCGACAATACCACCAGGAGtt                                                                                                                                     >  1:770599/1‑149 (MQ=255)
                 gTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTtat                                                                                                                                   <  2:208443/149‑1 (MQ=255)
                 gTATTAATGCCGGATTGCTCAATTATAATTCTAAAGGTGATTGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACACGGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAGTACCACCTGGGGTtat                                                                                                                                   >  1:469271/1‑149 (MQ=255)
                    ttAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGGTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAAc                                                                                                                                >  1:179516/1‑149 (MQ=255)
                       aTGCCGGATTGCTCAATTATAATTTCAACGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGAAAAAACCACCTGTAGTTATAACAGt                                                                                                                             >  2:758967/1‑149 (MQ=255)
                               ttGCTCAACTGTAACTGCAGCGGAAATAGTGTACAGAATCGTATTGGGGGTAACAGCCATTATGCATATGTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCTTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAg                                                                                                                     <  2:171613/149‑1 (MQ=255)
                                tGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAGACCTACAGAGTGGGTTAAATATTTGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCTACAgt                                                                                                                    <  1:514597/149‑2 (MQ=255)
                                 gCTCAATTAGAATTTCGGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGat                                                                                                                   <  2:838158/149‑1 (MQ=255)
                                   tCAATTATAATATCAGCGGAGATAGATTACAGAATCGGATTGGAGGTAACAACGATTATGCATATTTAAACCTACAGAGTGGGTTAAAAATTGGTGCGTGGCGTTGACGCGACAATACCACCTGTAGTTATAACAGTAGCGACAGatca                                                                                                                 >  1:540491/1‑149 (MQ=255)
                                     aaTTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTAGAGAGTGGGTTAAATATTGGTGCGTGGCGTTTATGCGGCAATAACACCTGAGATTATAACAGTAGCGACATatcaat                                                                                                               >  2:257202/1‑147 (MQ=255)
                                     aaTTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGGTATAACAGTAGAGACAGatcaca                                                                                                               >  2:257362/1‑147 (MQ=255)
                                                 aGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCaaaaa                                                                                                   <  2:794269/149‑1 (MQ=255)
                                                 aGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCaaaaa                                                                                                   <  1:688341/149‑1 (MQ=255)
                                                    ggAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAaaataa                                                                                                <  1:146976/149‑1 (MQ=255)
                                                      aaaTAGTGTACAGAAGCGGATTGGGGGTAACAGCTATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGTCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAaaataaat                                                                                              <  2:314843/149‑1 (MQ=255)
                                                      aaaGAGTGTACAGAATCGGAGTGGTGTTTACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAAGATTGCTGCGTGGCTTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAaaataaat                                                                                              <  2:239431/149‑1 (MQ=255)
                                                        aTAGTGTACAGAAACGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATgg                                                                                            <  2:654720/149‑1 (MQ=255)
                                                           gtgtACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGgcag                                                                                         <  1:475800/149‑1 (MQ=255)
                                                                 aGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATc                                                                                   >  1:675477/1‑149 (MQ=255)
                                                                  gAATCGGATTGGGCGTAACAGCCATTATGCATATTTAAACCCACAGAGGGCGTTAATTAATATTACAATGAAGTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAACAAATAGCAACATATCa                                                                                  >  2:441692/1‑149 (MQ=255)
                                                                    aTCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGGGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAAt                                                                                >  1:519808/1‑149 (MQ=255)
                                                                      cGGATTGGGGTTAACATCCATTATCCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATAc                                                                              <  2:15041/149‑1 (MQ=255)
                                                                          ttGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTgg                                                                          <  1:338954/149‑1 (MQ=255)
                                                                               ggTAAGCGCCATTATGCATATTTAAACCTGCGGAGGTGGTTAAATAGTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGa                                                                     <  1:174897/149‑1 (MQ=255)
                                                                                   aCAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCGTCCTGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCga                                                                 >  1:637879/1‑149 (MQ=255)
                                                                                   aCAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGTCTGGAGCga                                                                 <  2:783933/149‑1 (MQ=255)
                                                                                   aCAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCga                                                                 <  1:756205/149‑1 (MQ=255)
                                                                                   aCAGCCAGTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCga                                                                 >  1:670497/1‑149 (MQ=255)
                                                                                     aGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCgaga                                                               >  2:282814/1‑149 (MQ=255)
                                                                                     aGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGAAGCATATCAATACCTGGCTTGTGCggga                                                               >  2:414786/1‑149 (MQ=255)
                                                                                     aGCCATTATGCATATTTAAACCGACAGAGTCGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCgaga                                                               >  1:488401/1‑149 (MQ=255)
                                                                                      ggcATTAGGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGGAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGAc                                                              <  2:578507/147‑1 (MQ=255)
                                                                                      gCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGAAAATACCACCt                                                                                                                                           >  1:628570/1‑72 (MQ=255)
                                                                                      gCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGAAAATACCACCt                                                                                                                                           <  2:628731/72‑1 (MQ=255)
                                                                                      gCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGAAAATACCACCt                                                                                                                                           >  1:628731/1‑72 (MQ=255)
                                                                                      gCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGAAAATACCACCt                                                                                                                                           <  2:628570/72‑1 (MQ=255)
                                                                                          ttATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGAAAAAATAAATGGCAGCATATCAATACATGGCTTGAGCGAGACataa                                                          >  2:281578/1‑149 (MQ=255)
                                                                                            aTGCATATTTAAACCTACAGAGTGGGTTAGATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGATAGCACAAATAACTGGCAACATACCGATACCTGGCTTCAGCGAGACacaaca                                                        >  2:799343/1‑147 (MQ=255)
                                                                                             tGCATATTTAAACCTACAGAGTGGGGTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATAc                                                       <  1:436223/149‑1 (MQ=255)
                                                                                              ggctatTTGAACCTACAGAGTGGGTTAAATATTTGTGCGTGGCGTTCACGCGACAATACCACCTGAAGTTATAACAGTAGCGACAGATCATCAGGTAGCTAAAATAAATGGCAGTATATCAATACCTGGCTTGAGCGAGACataatctc                                                      <  2:671728/146‑4 (MQ=255)
                                                                                               cATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCg                                                     <  1:448606/149‑1 (MQ=255)
                                                                                                 tatTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGtt                                                   >  1:528068/1‑149 (MQ=255)
                                                                                                   tgtAAACCGACAGAGTTGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTAc                                                 <  2:893554/147‑1 (MQ=255)
                                                                                                           tACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCgg                                         >  1:348442/1‑149 (MQ=255)
                                                                                                           tACAGAGTGGGCTAAATATTGGTCCGCGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGCTAGCAACAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCgg                                         >  1:540173/1‑149 (MQ=255)
                                                                                                            aCAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTAAGCGAGACATAATACCGATACGTTCCCGGc                                        >  1:187241/1‑149 (MQ=255)
                                                                                                               gagTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGa                                     >  2:566981/1‑149 (MQ=255)
                                                                                                               gagTGGGTTAAATAGTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGa                                     >  1:565239/1‑149 (MQ=255)
                                                                                                                  cgggTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGc                                  <  1:347574/148‑1 (MQ=255)
                                                                                                                      ttAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTggg                              >  2:680050/1‑149 (MQ=255)
                                                                                                                      ttAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATAATTCACTATCGGTATTACGAACTAAAACCAGTACGTTCCCGGCTGACGCTggg                              >  1:914739/1‑149 (MQ=255)
                                                                                                                           tatTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGCCATAATACCGTTACGTTCCCGGCTGACGCGGCGTGATg                         >  2:33428/1‑149 (MQ=255)
                                                                                                                             ttGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCga                                                                 >  2:340213/1‑107 (MQ=255)
                                                                                                                             ttGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCga                                                                 <  1:340213/107‑1 (MQ=255)
                                                                                                                               ggTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTta                     <  1:891121/149‑1 (MQ=255)
                                                                                                                                     tGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCTAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCa               >  1:825375/1‑149 (MQ=255)
                                                                                                                                     tGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCa               <  1:441692/149‑1 (MQ=255)
                                                                                                                                               gcgACAATACCACGTGGAGTTATCACAGTAGCGACAGATGATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTAGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCGCACGCTGGGAGATGGTTATACTCAGGGCGgtagt     >  1:811483/1‑147 (MQ=255)
                                                                                                                                                  aCAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTc  <  2:917798/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
TGAGTTATGGGATCCCGGTATTAATGCCGGATTGCTCAATTATAATTTCAGCGGAAATAGTGTACAGAATCGGATTGGGGGTAACAGCCATTATGCATATTTAAACCTACAGAGTGGGTTAAATATTGGTGCGTGGCGTTTACGCGACAATACCACCTGGAGTTATAACAGTAGCGACAGATCATCAGGTAGCAAAAATAAATGGCAGCATATCAATACCTGGCTTGAGCGAGACATAATACCGTTACGTTCCCGGCTGACGCTGGGTGATGGTTATACTCAGGGCGATATTTTC  >  NC_000913_3_pae_tpiA/4545536‑4545830

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