Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 4,185,526 G→T K87N (AAG→AAT)  rpoC → RNA polymerase, beta prime subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA4,185,5260GT82.3% 256.2 / 61.1 113K87N (AAG→AAT) rpoCRNA polymerase, beta prime subunit
Reads supporting (aligned to +/- strand):  ref base G (7/13);  new base T (43/50);  total (50/63)
Fisher's exact test for biased strand distribution p-value = 4.59e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.31e-01

Only 100 of 113 total aligned reads displayed.

AAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCG  >  NC_000913_3_pae_tpiA/4185380‑4185671
                                                                                                                                                  |                                                                                                                                                 
aaaaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAGGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgc                                                                                                                                                >  1:748625/1‑150 (MQ=255)
aaaaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgc                                                                                                                                                >  3:489622/1‑150 (MQ=255)
aaaaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCGTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgc                                                                                                                                                <  4:674722/150‑1 (MQ=255)
  aaaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTGTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgcgg                                                                                                                                              <  1:427113/150‑1 (MQ=255)
   aaGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgcggc                                                                                                                                             >  1:1191486/1‑150 (MQ=255)
      ccGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGtt                                                                                                                                          <  1:1085684/150‑1 (MQ=255)
       cGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgcggc                                                                                                                                             <  2:179981/146‑1 (MQ=255)
       cGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgcggc                                                                                                                                             >  1:179981/1‑146 (MQ=255)
        ggAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGa                                                                                                                                        >  2:1251677/1‑150 (MQ=255)
         gAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGaa                                                                                                                                       >  2:92939/1‑150 (MQ=255)
         gAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGaa                                                                                                                                       <  2:434192/150‑1 (MQ=255)
         gAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGaa                                                                                                                                       <  1:311336/150‑1 (MQ=255)
         gAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGaa                                                                                                                                       <  1:1074266/150‑1 (MQ=255)
            aCCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTg                                                                                                                                    <  1:88367/150‑1 (MQ=255)
             ccATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGa                                                                                                                                   <  2:1191486/150‑1 (MQ=255)
              cATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGa                                                                                                                                        <  3:899380/144‑1 (MQ=255)
              cATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGAc                                                                                                                                  <  1:360666/150‑1 (MQ=255)
               aTCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGAcc                                                                                                                                 <  1:34401/150‑1 (MQ=255)
                   aCTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGa                                                                                                                             <  2:1068354/150‑1 (MQ=255)
                    cTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGAc                                                                                                                            <  3:94056/150‑1 (MQ=255)
                         gTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGa                                                                                                                             >  2:256665/1‑144 (MQ=255)
                         gTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAg                                                                                                                       >  4:385735/1‑150 (MQ=255)
                           tttcaggaaCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTa                                                                                                                     >  2:103556/8‑150 (MQ=255)
                             gTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCATATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACg                                                                                                                   <  2:1310358/150‑1 (MQ=255)
                                cAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCg                                                                                                                >  1:48581/1‑150 (MQ=255)
                                 aaaCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGt                                                                                                               >  1:467258/1‑150 (MQ=255)
                                   aCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGa                                                                                                             <  1:545179/150‑1 (MQ=255)
                                   aCCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGa                                                                                                             >  1:387772/1‑150 (MQ=255)
                                    ccAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAg                                                                                                            >  3:800522/1‑150 (MQ=255)
                                     cAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGc                                                                                                           >  3:615129/1‑150 (MQ=255)
                                     cAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGc                                                                                                           >  2:154282/1‑150 (MQ=255)
                                           gTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTTCGCCGTGAGCGTATgg                                                                                                     <  2:748625/150‑1 (MQ=255)
                                           gTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATgg                                                                                                     <  2:467258/150‑1 (MQ=255)
                                            tGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATggg                                                                                                    <  2:1160312/150‑1 (MQ=255)
                                             gACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGc                                                                                                   >  1:154125/1‑150 (MQ=255)
                                             gACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGc                                                                                                   >  2:521736/1‑150 (MQ=255)
                                              aCGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGcc                                                                                                  <  3:520760/150‑1 (MQ=255)
                                               cGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCca                                                                                                 >  1:1118194/1‑150 (MQ=255)
                                                ggCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcggc                                                                                                                                             <  3:557170/105‑1 (MQ=255)
                                                ggCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTgcggc                                                                                                                                             >  4:557170/1‑105 (MQ=255)
                                                  ccTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACAt                                                                                              >  4:67264/1‑150 (MQ=255)
                                                    ttttCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCg                                                                                            <  4:15344/150‑1 (MQ=255)
                                                    ttttCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCg                                                                                            <  2:703896/150‑1 (MQ=255)
                                                      ttCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGaa                                                                                          <  4:340618/150‑1 (MQ=255)
                                                          gcgcCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTgg                                                                                      <  1:311191/150‑1 (MQ=255)
                                                          gcgcCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTgg                                                                                      <  2:651693/150‑1 (MQ=255)
                                                              ccGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgcggc                                                                                                                                             >  2:411718/1‑91 (MQ=255)
                                                              ccGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTgcggc                                                                                                                                             <  1:411718/91‑1 (MQ=255)
                                                              ccGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTc                                                                                  >  2:820986/1‑150 (MQ=255)
                                                              ccGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTc                                                                                  >  3:180361/1‑150 (MQ=255)
                                                              ccGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTc                                                                                  >  1:1281580/1‑150 (MQ=255)
                                                                gTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTccc                                                                                >  1:1439207/1‑150 (MQ=255)
                                                                         ggCCGTTAAAAGATTACGCGTGCCTGTTCGGTAAGTACGAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCATACTAAAGTACGCCGGTATCGTATGGTCCACATCGAACTGGCTTCCCCGTCTgcgc                                                                       <  2:728317/150‑1 (MQ=255)
                                                                         ggCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTgcgc                                                                       >  4:204347/1‑150 (MQ=255)
                                                                          gCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACCACGTAGCGGCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGca                                                                      <  1:92939/150‑1 (MQ=255)
                                                                             ggTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGt                                                                                                                      <  4:303628/99‑1 (MQ=255)
                                                                             ggTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGt                                                                                                                      >  3:303628/1‑99 (MQ=255)
                                                                              gTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATc                                                                  <  1:1244540/150‑1 (MQ=255)
                                                                              gTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATc                                                                  <  4:800522/150‑1 (MQ=255)
                                                                               tAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCt                                                                 <  1:751451/150‑1 (MQ=255)
                                                                                 aaaGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTgg                                                               <  1:892252/150‑1 (MQ=255)
                                                                                   aGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGtt                                                             <  3:9746/150‑1 (MQ=255)
                                                                                   aGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGtt                                                             >  4:803498/1‑150 (MQ=255)
                                                                                    gATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTc                                                            >  4:114724/1‑150 (MQ=255)
                                                                                          gAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACGCCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGaaa                                                      <  4:474374/150‑1 (MQ=255)
                                                                                          gAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGaaa                                                      <  1:103556/150‑1 (MQ=255)
                                                                                            gTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATc                                                    <  1:103939/150‑1 (MQ=255)
                                                                                               ccTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATgg                                                                                                     <  1:602001/98‑1 (MQ=255)
                                                                                               ccTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATgg                                                                                                     >  2:602001/1‑98 (MQ=255)
                                                                                               ccTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCt                                                 <  3:67264/150‑1 (MQ=255)
                                                                                                cTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTg                                                >  4:590950/1‑150 (MQ=255)
                                                                                                     cGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAGAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTc                                           <  4:580694/150‑1 (MQ=255)
                                                                                                      ggTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTcc                                          <  2:48581/150‑1 (MQ=255)
                                                                                                       gTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTccc                                         <  1:880546/150‑1 (MQ=255)
                                                                                                         aaGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCgc                                       >  1:142532/1‑149 (MQ=255)
                                                                                                           gTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGt                                                                                                         <  2:1020624/82‑1 (MQ=255)
                                                                                                           gTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGt                                                                                                         >  1:1020624/1‑82 (MQ=255)
                                                                                                           gTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATg                                                                                                      <  1:1184369/85‑1 (MQ=255)
                                                                                                           gTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATg                                                                                                      >  2:1184369/1‑85 (MQ=255)
                                                                                                           gTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTAt                                     <  1:154282/150‑1 (MQ=255)
                                                                                                               aaGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAg                                                                                                                                      <  1:1411078/49‑1 (MQ=255)
                                                                                                               aaGCGCCTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAg                                                                                                                                      >  2:1411078/1‑49 (MQ=255)
                                                                                                                     cTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTctgctg                           >  2:348383/1‑150 (MQ=255)
                                                                                                                     cTGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTctgctg                           >  3:759092/1‑150 (MQ=255)
                                                                                                                      tGAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCtgctgc                          >  4:885251/1‑150 (MQ=255)
                                                                                                                       gAAACACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTgctgct                         <  2:1281580/150‑1 (MQ=255)
                                                                                                                        aaaCACCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCt                                                 <  2:1205749/125‑1 (MQ=255)
                                                                                                                          acacCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGa                      <  3:82455/150‑1 (MQ=255)
                                                                                                                            acCGTGGCGTCATCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGata                    <  4:615129/150‑1 (MQ=255)
                                                                                                                               gTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGc                 <  4:776106/150‑1 (MQ=255)
                                                                                                                                       aTCTGTGAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgcgc         <  3:275046/150‑1 (MQ=255)
                                                                                                                                       aTCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTgcgc         <  1:865740/150‑1 (MQ=255)
                                                                                                                                           gtgAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGc                 <  4:406593/138‑1 (MQ=255)
                                                                                                                                           gtgAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGc                 >  3:406593/1‑138 (MQ=255)
                                                                                                                                           gtgAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGata     <  1:732425/150‑1 (MQ=255)
                                                                                                                                            tgAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATc                                    <  1:959290/118‑1 (MQ=255)
                                                                                                                                            tgAGAATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATc                                    >  2:959290/1‑118 (MQ=255)
                                                                                                                                             gagaATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATc   <  1:452652/150‑1 (MQ=255)
                                                                                                                                             gagaATTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATc   <  3:25581/150‑1 (MQ=255)
                                                                                                                                              agaAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCg  >  2:899830/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                 
AAAAAGCCGGAAACCATCAACTACCGTACGTTCAAACCAGAACGTGACGGCCTTTTCTGCGCCCGTATCTTTGGGCCGGTAAAAGATTACGAGTGCCTGTGCGGTAAGTACAAGCGCCTGAAACACCGTGGCGTCATCTGTGAGAAGTGCGGCGTTGAAGTGACCCAGACTAAAGTACGCCGTGAGCGTATGGGCCACATCGAACTGGCTTCCCCGACTGCGCACATCTGGTTCCTGAAATCGCTGCCGTCCCGTATCGGTCTGCTGCTCGATATGCCGCTGCGCGATATCG  >  NC_000913_3_pae_tpiA/4185380‑4185671

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