Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913_3_pae_tpiA 482,230 G→A P725L (CCG→CTG)  acrB ← multidrug efflux system protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913_3_pae_tpiA482,2300GA98.9% 301.0 / NA 93P725L (CCG→CTG) acrBmultidrug efflux system protein
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base A (52/40);  total (53/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_pae_tpiA/482087‑482375
                                                                                                                                                |                                                                                                                                                 
catAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTa                                                                                                                                              >  2:1401590/1‑149 (MQ=255)
 atAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCACCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAACCTTCAGGGTAt                                                                                                                                             >  2:679608/1‑149 (MQ=255)
      cTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATGGATTCTCCCGGCGTCTGCGCTTTTTCCTTGTCGATATCAATCTTAAACTGCAGGGTATCTTcc                                                                                                                                        <  2:887585/149‑1 (MQ=255)
         tcttcACCCGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGGCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGa                                                                                                                                     <  1:600680/149‑1 (MQ=255)
          cttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAc                                                                                                                                    >  2:839458/1‑149 (MQ=255)
          ctgcacaCGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGCGTGGTGTTAATGTCGTTGATAGAAACCCCCAGCGCCAGCGCTGTTTCCTGGTCGATATCAGTATTAACCTGCGTGGTATCTGccggtc                                                                                                                                    >  2:838053/4‑145 (MQ=255)
           ttcACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCAGACCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGAcc                                                                                                                                   >  1:1290586/1‑149 (MQ=255)
              acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGtt                                                                                                                                >  1:922788/1‑149 (MQ=255)
              acacGACCGCGGTCGATAAAGTCGTTCACATAGCTGCAGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATTAATCTTAAACTGCAGGGTATCTTCCAGCCCGtt                                                                                                                                >  1:748181/1‑149 (MQ=255)
               cacGACCGCGGTCGATAAAGTAGTTGACATAGCTGCTGCCAAGAAAGCGTGTTGATGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGttt                                                                                                                               <  1:691325/149‑1 (MQ=255)
                      gcgGGCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTa                                                                                                                        <  2:360037/149‑1 (MQ=255)
                       cgGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGAGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTAc                                                                                                                       <  1:1395127/149‑1 (MQ=255)
                            gATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTgg                                                                                                                  >  1:364875/1‑149 (MQ=255)
                              tAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTc                                                                                                                >  2:85942/1‑149 (MQ=255)
                                 aGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTGGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAAc                                                                                                             >  2:1403945/1‑149 (MQ=255)
                                 aGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAAc                                                                                                             >  1:626910/1‑149 (MQ=255)
                                  gTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                            >  1:1110476/1‑149 (MQ=255)
                                  gTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACa                                                                                                            <  1:483762/149‑1 (MQ=255)
                                        cacaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAACCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAg                                                                                                      <  2:12454/149‑1 (MQ=255)
                                         acaTAGCTGCCGCCCCATGCAGCGCCCCGAGTGGTGTTAATGTCGTTGATAGAAGCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAATGGCAGCGTATATTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  1:209102/149‑1 (MQ=255)
                                         acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCTgg                                                                                                     >  2:953738/1‑149 (MQ=255)
                                         acaTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  1:742072/149‑1 (MQ=255)
                                          caTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAggg                                                                                                    <  1:69966/149‑1 (MQ=255)
                                           aTAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTCCGCTGGTCAACATATCAGAGt                                                                                                   >  2:1435156/1‑149 (MQ=255)
                                             aGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAg                                                                                                      >  2:207974/1‑144 (MQ=255)
                                             aGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAg                                                                                                      <  1:207974/144‑1 (MQ=255)
                                             aGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGc                                                                                                 >  1:1338808/1‑149 (MQ=255)
                                              gCTGCCGCCCCATGCAGCGCCCAGAGTGGGGTTAATGTCGTTGATAGAAAAACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAAACTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGgctgct                                                                                                >  1:1222369/1‑149 (MQ=255)
                                                   cgccCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCgct                                                                                           <  2:364875/149‑1 (MQ=255)
                                                        cATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTc                                                                                      <  2:1191116/149‑1 (MQ=255)
                                                             aGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaa                                                                                 <  1:749347/149‑1 (MQ=255)
                                                              gcgcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                                >  1:480462/1‑149 (MQ=255)
                                                               cgcCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGAGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCTGACCGTTTGGACGTACGCTGGTCAACATATCAGGGGGCTTCACTGTCTCTgcaagc                                                                               >  1:1094777/1‑149 (MQ=255)
                                                                 cccAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaagcaa                                                                             >  2:561892/1‑149 (MQ=255)
                                                                    agagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTg                                                                          <  1:85942/149‑1 (MQ=255)
                                                                    agagTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACCCAGATCCAAATATCCGCGTGCTTCGCTACTGCTTCAAGCAAcgt                                                                          >  2:280330/1‑147 (MQ=255)
                                                                           gtgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATAAGGGTGCTTCGCTGCTTCTGCAACCACCTGGTTAcg                                                                    >  2:1238917/1‑149 (MQ=255)
                                                                            tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   <  2:616656/149‑1 (MQ=255)
                                                                            tgtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTAcgc                                                                   >  2:392158/1‑149 (MQ=255)
                                                                             gtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     <  2:798034/114‑1 (MQ=255)
                                                                             gtTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAgg                                                                                                     >  1:798034/1‑114 (MQ=255)
                                                                              ttAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACgcgc                                                                 <  1:1401590/149‑1 (MQ=255)
                                                                                aaTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCGACTGGTTCCGCCCCt                                                               >  1:868805/1‑149 (MQ=255)
                                                                                 aTGTCGTTGATAGAAAGAAACAGCCCCTGCGCTTTTTCCAGGCGAACAGCAATCTTAAACTGAACGGTGACCAGCAGACCGATTGGACGTACACTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACCCGCCTg                                                              >  2:305300/1‑149 (MQ=255)
                                                                                 aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                                <  1:583033/131‑1 (MQ=255)
                                                                                 aTGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTgcaag                                                                                >  2:583033/1‑131 (MQ=255)
                                                                                  tGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                             <  1:1240400/149‑1 (MQ=255)
                                                                                  tGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGa                                                             >  1:1412681/1‑149 (MQ=255)
                                                                                    tCGTTGATAGAAACGCCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGTACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                           <  2:392914/149‑1 (MQ=255)
                                                                                    tCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGt                                                           >  1:1263197/1‑149 (MQ=255)
                                                                                    tCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTAGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGCACGTACGCTGGTCAACATATCAGGGTGCTTCGCTCCTTCTGCAAGCAACTGGTTCCGCGCCTGAGt                                                           >  2:70964/1‑149 (MQ=255)
                                                                                    tCGTTGATAGAAACACCCAACGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGCTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCACGCAACTGGTTACGCGCCTGact                                                           >  2:1221090/1‑147 (MQ=255)
                                                                                      ggtGATAGAAACACCCAGCGCCTGCGGTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                         <  2:1172009/147‑1 (MQ=255)
                                                                                      gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAACCAGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                         >  1:279801/1‑149 (MQ=255)
                                                                                      gTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCa                                                         <  1:489258/149‑1 (MQ=255)
                                                                                       ttGATAGAAACAACCAGCGCCTGCGCTATTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATATTCGAGACCGTTTGAACGTACGCTGGTCAACATATCAAGGAGCGTTGCTGGTTCTTCAAGCAACTGGTGTCGCGCCGGAGCCAg                                                        >  2:334589/1‑149 (MQ=255)
                                                                                         gATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGtt                                                      >  1:1306334/1‑149 (MQ=255)
                                                                                          aTAGATACACCCAGCGCCTGCGCTTTCTCCTGGTGGATATCAATCTGAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCGGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGttt                                                     <  2:255315/149‑1 (MQ=255)
                                                                                             gacacccccAGCGCCTGCGCTTTTTCCTGGTCGATATCAAGCTTAAACTGCAGGGGATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                  <  1:1411932/143‑1 (MQ=255)
                                                                                             gAAACACCCAGCGCCTGCGCTCTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTc                                                  <  1:679608/149‑1 (MQ=255)
                                                                                              aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTg                                                                                                                              <  1:528903/72‑1 (MQ=255)
                                                                                              aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTg                                                                                                                              >  2:528903/1‑72 (MQ=255)
                                                                                              aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                 >  1:291711/1‑149 (MQ=255)
                                                                                              aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                 >  2:1293550/1‑149 (MQ=255)
                                                                                              aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCg                                                 >  1:980277/1‑149 (MQ=255)
                                                                                              aaaCACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGACAGTTTTTCg                                                 >  1:865774/1‑149 (MQ=255)
                                                                                                acacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTg                                               <  1:387587/149‑1 (MQ=255)
                                                                                                 cacCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGa                                              <  2:75146/149‑1 (MQ=255)
                                                                                                  acCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTCCGCTGGTCAACATATCAAGGTGCTTCGCTGCTGCTGCAAGCAACTGGTTCCGCGCCTGAGTCAGTTTTTCATgca                                             >  1:481720/1‑147 (MQ=255)
                                                                                                     cAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGAGGTACGCTGTTCAACATAACAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCTTGAGCTAGTTTTTCGTGACCaa                                          >  1:1283390/1‑149 (MQ=255)
                                                                                                     cAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTAAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCATGAGTCAGTTTTTTGTGACCaa                                          >  1:299015/1‑149 (MQ=255)
                                                                                                       gcgcCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAgg                                        <  1:561892/149‑1 (MQ=255)
                                                                                                          ccgggggtttttCCTGGTCGATATCAATCTTAAATTGCAGGGTATCTTCCAGACCGTTTGGAGGTACGCTGGTCAACATATCAGGGTGCTACGCTGCTTCGGCAAGCACCTGGTTATGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCa                                     <  2:310922/142‑1 (MQ=255)
                                                                                                           cTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGCCCGTTTGGACGTACGATGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTGCGCGCCTGTGTCAGTTTTTCGTGACCAAGGCCAg                                    >  2:1210903/1‑149 (MQ=255)
                                                                                                             gcgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGcc                                  >  1:585565/1‑149 (MQ=255)
                                                                                                              cgcTTTTTCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCt                                 <  2:174567/149‑1 (MQ=255)
                                                                                                               gcTTGTTCCTGGTCGATATCAATCGTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTg                                <  2:118486/149‑1 (MQ=255)
                                                                                                                ttttttCCTGGTCGATATCAATCGTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTgg                               <  1:274/148‑1 (MQ=255)
                                                                                                                    ttCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                           <  1:839458/149‑1 (MQ=255)
                                                                                                                    ttCCTGGTCGATATCAATCTTAAACTGCAGGGGATCTTCCAGACCGTTTGGACGTACGCTGGTCAACCTATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCaa                           <  1:838053/149‑1 (MQ=255)
                                                                                                                     tCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAAt                          <  2:1051450/149‑1 (MQ=255)
                                                                                                                     tCCTGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCACTTGTTCGTGCCCAAGGCCAGCCTGGTCAAt                          <  1:1403945/149‑1 (MQ=255)
                                                                                                                        tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAg                       >  1:516469/1‑149 (MQ=255)
                                                                                                                        tGGTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGCGACCAAGGCCAGCCTGGTCAATCAg                       >  1:391677/1‑149 (MQ=255)
                                                                                                                          gTCGATATCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCt                     <  1:585732/149‑1 (MQ=255)
                                                                                                                                 tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTTCCCAAGGCCAGCCTGGTCAACCAGCtcaaagtc              >  1:953943/1‑149 (MQ=255)
                                                                                                                                 tCAATCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtc              >  2:801824/1‑149 (MQ=255)
                                                                                                                                   aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCtcaaagtcaa            >  1:421933/1‑149 (MQ=255)
                                                                                                                                   aaTCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTAAACCTATTAGGGTGTGTCTCTGCTTCTGCAAGCAACTGGTTACGCGCCTGGGTTGTATTTTCGTGTCCAAGTCCGGCCTGGTCAATCAGCttgaagtcaa            >  1:1208720/1‑149 (MQ=255)
                                                                                                                                     tCTTAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCGAGGCCAGCCTGGGCAATCAGCtcaaagtcaaag          >  1:1194743/1‑149 (MQ=255)
                                                                                                                                       ttAAACTGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGcc        >  1:446671/1‑149 (MQ=255)
                                                                                                                                        tacaCTGCAGGGTATCGTCCAGTCCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCg       <  2:293013/146‑1 (MQ=255)
                                                                                                                                             tGCAGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGc  >  1:807561/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
CATAAACTTTCTTCACACGACCGCGGTCGATAAAGTCGTTCACATAGCTGCCGCCCCATGCAGCGCCCAGAGTGGTGTTAATGTCGTTGATAGAAACACCCAGCGCCTGCGCTTTTTCCTGGTCGATATCAATCTTAAACTGCGGGGTATCTTCCAGACCGTTTGGACGTACGCTGGTCAACATATCAGGGTGCTTCGCTGCTTCTGCAAGCAACTGGTTACGCGCCTGAGTCAGTTTTTCGTGACCAAGGCCAGCCTGGTCAATCAGCTCAAAGTCAAAGCCGGTTGC  >  NC_000913_3_pae_tpiA/482087‑482375

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: